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ABSTRACT 


\ 

\ 

\ ^ 

A  synthesis  method  based  on  geometrical  optics  for  designing  a  dual 
reflector  antenna  system  with  an  arbitrary  phase  and  amplitude  distri¬ 
bution  in  the  aperture  of  the  second  reflector  is  presented.  The  first 
reflector  may  be  illuminated  by  a  pattern  with  an  arbitrarily  curved 
phase  front.  A  pair  of  first  order  ordinary  nonlinear  differential 
equations  of  the  form  dy/ dx  =  f(x,  y)  are  developed  for  the  system.' 
Questions  concerning  uniqueness,  existence,  and  bounds  for  the  solutions 
are  discussed.  Calculations  and  numerical  results  for  the  design  of 
a  uniform  amplitude  and  phase  dual-reflector  system  are  presented. 


The  diffraction  effects  of  the  small  reflector  are  analyzed  by  the 
methods  of  geometrical  diffraction.  Their  effects  upon  the  aperture 
distribution  of  the  larger  reflector  are  analyzed  in  detail.  Correction 
for  the  small  reflector  diffracted  field  is  obtained  partially  by  an 

jk,  h  .  h  i 


iterative  procedure  utilizing  the  above  synthesis  method 


.  i , 


/% 


eo^l^iwelys  ^''extensive  numerical  analysis  that  com¬ 
plete  correction  for  the  diffraction  effects  is  not  possible. 
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GENERAL  INTRODUCTION 

Conventional  dual- reflector  antenna  systems  have  been  based  largely  on 

the  Cassegrain  parabola-hyperbola  design  or  the  Gregorian  parabola-ellipse 

1  2 
design  .  Some  highly  specialized  exceptions  have  also  been  reported  . 

The  designs  are  all  based  on  the  principles  of  geometrical  optics  and  are 

limited  accordingly.  That  is,  the  reflectors  must  be  large  and  have  a  large 

radius  of  curvature  compared  to  the  wavelength  in  addition  to  other  restrictions. 

A  generalization  of  the  design  techniques  for  a  dual  reflector  antenna  that  is 
also  based  on  geometrical  optics  and  is  similarly  limited  in  application  is 
presented  in  Chapter  1.  It  has  been  found  that  an  arbitrary  phase  and  ampli¬ 
tude  distribution  can  be  developed  in  the  aperture  of  the  larger  reflector  of 
a  dual- reflector  system  with  an  arbitrarily  curved  phase  front  illuminating 
the  smaller  reflector.  The  design  procedure  is  sufficiently  general  so  that 
many  useful  variations  of  the  design  objectives  are  clearly  possible-  For 
example,  the  large  reflector  of  the  system  may  be  specified  to  have  a  circular 
cross-section  and  arbitrary  phase  distribution  (uniform, for  example)  and  the 
appropriate  sub- reflector  is  then  found  for  a  given  primary  radiation  pattern. 

The  methods  of  geometrical  optics  which  are  used  for  the  above  synthesis  are 
of  the  first  order  from  the  point  of  view  of  the  assymptotic  solution  of 
Maxwell's  equations,,  which  is  an  expansion  in  powers  of  the  wavelength. 

Second  and  higher  order  solutions  have  been  obtained  for  many  problems 
and  yield  the  diffraction  effects.  A  discussion  of  the  assymptotic  expansion 
and  the  use  of  higher  order  terms  in  the  geometrical  theory  of  diffraction  is 
made  in  Chapter  2,  This  t  heory  is  then  used  to  analyse  the  diffracted  field  of 
the  subreflector.  The  effect  of  this  diffracted  field  upon  the  aperture  distribution 


2 


of  the  larger  reflector  is  discussed  in  Chapter  3  in  detail,,  together  with  a 
presentation  of  numerical  results  illustrating  the  effects*  The  importance 
of  higher  order  terms  or  the  diffracted  field  due  to  the  larger  reflector  is 
generally  not  nearly  as  significant  as  is  the  diffraction  due  to  the  subreflector.. 
In  addition,  as  the  subreflector  is  made  smaller  in  wavelength,  diffraction 
becomes  more  significant  and  places  a  limitation  on  the  minimum  size  and 
consequently  the  minimum  blockage  the  subreflector  presents  to  the  main  or 
larger  reflector. 

An  attempt  at  including  the  diffracted  field  in  the  synthesis  procedure  and  hence 
correcting  for  it  exactly  proved  unsuccessful.  It  appears  that  in  the  context 
of  the  present  formulation  a  solution  for  this  problem  does  not  exist.  The 
complexity  of  the  problem  makes  a  proof  of  no  existence  difficult  although 
this  conclusion  is  strongly  implied  by  the  discussion  of  Chapter  4.  Since  an 
exact  solution  did  not  appear  possible,  an  attempt  at  an  Iterative  solution 
using  the  synthesis  method  developed  in  Chapter  1  was  made  and  is  discussed 
In  Chapter  5,  The  numerical  results  obtained  here  tended  to  support  the 
implications  made  in  Chapter  4  that  no  exact  solution  exists.  However, 
the  results  of  the  iterative  procedure  indicated  that  a  reasonable  compensation 
for  the  diffraction  effects  can  be  obtained.  Although  an  approximate  theory 
of  diffraction  was  used  in  this  iterative  procedure,  the  true  diffracted  field 
could,  in  principle,  also  be  used  in  this  method;  although  perhaps  with  con¬ 
siderable  inconvenience  in  proportion  to  the  increased  accuracy  obtainable. 


CHAPTER  1 


FIRST  ORDER  DUAL  REFLECTOR  SYNTHESIS 

Introduction 

In  this  chapter  a  method  for  the  synthesis  of  a  dual  reflector  antenna 
with  an  arbitrarily  specified  phase  and  amplitude  distribution  in  the 
aperture  of  the  main  reflector  is  developed.  An  arbitrary  primary 
source  may  be  used  to  illuminate  the  subreflector.  The  design  is 
based  ppon  geometrical  optics.  This  technique  for  synthesis  is  used 
later  to  obtain  an  iterative  solution  for  the  reflector  shapes  which 
accounts  for  the  diffraction  field  of  the  subreflector  to  some  extent. 


The  synthesis  method  utilizes  the  analytical  expression  of  the  geo¬ 
metric  optics  principles  together  with  the  geometry  for  the  reflectors 
illustrated  in  Figure  1  to  develop  a  pair  of  first-order  nonlinear 
ordinary  differential  equations  of  the  form 


dy 

dx 


f(x,  y) 


(1) 


which  leads  to  the  cross-sections  of  each  reflector  when  subject  to 
boundary  conditions  such  as 

y  (x  =  xmax)  =  0  ,  (2) 


The  above  differential  equation  can,  in  general,  be  solved  readily  by 
high-speed  machine  computations.  A  pair  of  such  solutions  for  a 
uniform  phase  and  amplitude  design  is  presented. 
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A  desirable  feature  about  the  form  of  the  differential  solution, 

Equation  (1),  is  that  frequently  much  information  can  be  predicted 
about  the  final  solution  before  a  machine  computation  is  attempted. 

The  questions  of  existence  and  upper  bounds  on  the  size  of  the 
reflector  (j  yjrnax)  can  be  answered  by  considering  little  more  than 
the  boundary  conditions  of  the  problem.  If  f(x,  y)  is  found  in  a  suitable 
form,  then  lower  bounds  for  the  radius  of  curvature  of  the  reflectors 
can  also  be  found,  The  procedures  for  finding  the  bounds,,  existence 
and  uniqueness  proofs  are  discussed  in  Part  III  and  can  also  be  found 
in  the  literature  on  non-linear  differential  equations  (see  Reference  8 
for  example)-  See  also  Reference  9  for  an  alternative  formulation 
of  the  problem. 

The  optical  principles  which  are  utilized  to  develop  equations  of  the 

form  given  in  (1)  for  the  system  are  the  following: 

a*.  Snell's  Law^,  Application  of  this  law  at  each  reflector  leads 

to  a  form  of  (1)  with  f  specified  in  terms  of  the  angular  variables 
shown  in  Figure  1, 

b.  Conservation  of  Energy  Flow^  along  the  ray  trajectories.  The 
requirement  that  the  energy  flow  be  conserved  in  any  solid  angle 
bounded  by  ray  trajectories  leads  to  one  equation  of  constraint 
for  the  system, 

c.  The  surfaces  of  Constant  Phase  form  Normal  Surfaces  ^  to  the 
ray  trajectories,  and  this  normal  congruence  is  maintained  after 
any  number  of  reflections  ^  (theorem  of  Malus),  This  principle 
leads  to  another  constraint  for  the  system. 


It  can  be  shown  that  all  these  principles  are  not  independent  (References 
2  -  6).  Hence  we  will  not  need  to  use  the  expression  for  Snell’s  law 


at  reflector  #2  (Equation  5),  The  optical  principles  together  with  the 
geometrical  bounds  of  the  system  and  the  primary  pattern  and  the  phase 
and  amplitude  distribution  in  the  large  reflector  define  a  unique  opti¬ 
cal  system  (sometimes  equivalent  alternative  systems  exist,  as  with 
the  Cassegrain  and  Gregorian  systems). 

The  design  method  is  described  for  surfaces  of  revolution,  in  particular 
although  the  antenna  design  method  can.be  readily  extended  to  include 
doubly  curved  surfaces  where  the  surface  coordinates  are  separable, 
and  possibly  also  under  other  special  conditions  (see  Reference  10). 
When  the  curvature  in  one  plane  is  specified  for  the  dual-reflector 
system,,  the  design  procedure  for  the  curvature  in  the  perpendicular 
plane  can  be  the  same  as  reported  herein  with  an  extension  similar  to 
that  made  by  Dunbar  for  a  single  doubly  curved  reflector. 


II.  Analytic  Expressions  of  Optical  Principles  and  Formal  Solution* 

The  principles  stated  earlier  can  be  expressed  as  the  following  basic 
equations  from  which  the  solution  is  developed. 

Snell's  law  at  reflector  #1  leads  to  (see  Figure  2) 


where 


dyi 


dx 


1 


=  tan 


e2  (xj  .ylr  x2,  y2) 
2 


=  arc  tan 


'  X2  "  X1  ^ 
«+?+yr  vi ) 


(3) 


(4) 


The  solution  presented  in  this  part  and  much  of  the  remainder  of  this 
chapter  is  published  in  Reference  11, 
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The  quantities  fl^(0j)  and  ^(Q\)  are  prescribed  such  that  e(  (0j)  +  ^(6^)  = 
given  constant  (see  Figure  1).  The  intersection  of  a  primary  field  ray 
with  the  yj  axis  determines  (0^)  and  f  (@l).  The  rays  are  found  from 
the  known  or  given  form  of  the  primary  field. 

Snell's  law  at  reflector  #2  leads  to  (see  Figure  3) 


where 


0j  =  arc  tan  | 


x2  ~  x2 
YZ 


(5) 


(6) 


The  conservation  of  energy  principle  is  exp?  ssed  differentially  as 

sin  61  ("aSEfj  =  X3p  (x2}  x2 

The  quantity  I j ( ® i )  is  the  power  density  of  the  primary  illumination. 
The  quantity  the  power  density  flow  normal  to  the  aperture 

of  reflector  #2  It  is  chosen  arbitrarily  except  for  the  required 

normalization: 


e 


V  1 

^  1  max 

0, 

1  mm 


lj(0j) sin  0 j  d0 j  = 


f 


2  max 


hplx2>  x2  dx2 


(8) 


2min 


The  phase  at  the  aperture  can  be  assigned  arbitrarily  as  C  (xL) 

P  L 

The  path  length  from  the  primary  illumination  is  given  by 


Cp(x2)  =  rl(er  yl}  +  p2(ei,yl*  XZ  y2)+P3(x2’  y2*x2)  +  Cpft(ei)  (9) 
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The  quantity  C  (0,  )  is  defined  by  the  primary  illumination  as  illus- 
PO  1 

trated  in  Figure  4.  An  arbitrary  primary  illumination 


t(0j)  pha 


se 


front  is  chosen  as  a  reference.  The  aperture  phase  C  (xl)  must  be 

P  c 

defined  with  respect  to  this  reference.  The  remaining  quantities  in 
Equation  (9)  are  easily  found  as 


*  M1  -  u 


sec  0 


r 


(10) 


and 


(U) 


(12) 


The  right-hand  side  of  Equation  (12)  is  found  by  expressing  the  theorem 
of  Malus  from  the  diagram  of  Figure  3  as 


dC 


dxl 


sin  0^ 


(13) 


Since  C  (xp  is  assumed  given,  then 


e3<*2> 


is  thus  determined 


directly. 


We  will  now  choose  as  a  single  independent  variable  the  quantity  x^  . 
The  remaining  dependent  variables  will  be  considered  functions  of 
X£.  We  will  then  derive  the  two  differential  equations 


d0, 


34"  =  fe,(ei-2> 


(14) 


and 
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j^l 

dx2 


yijxz>- 


(15) 


Upon  solving  these  equations  all  the  remaining  dependent  variables  can 
be  immediately  found  as  functions  of  Although  the  choice  of  xl, 

as  the  independent  variable  and  0^  and  for  the  differential  equations 
is  somewhat  arbitrary,  this  choice  does  allow  all  the  following  expres¬ 
sions  (f q  ,f^;etc.  )  to  be  found  in  ^uplicit  form.  An  equivalent  formula¬ 
tion  in  terms  of  dy^/dx^  =  ^(xpYj)  and  dy^/ dx^  =  can  a^SQ 

be  found  (see  Reference  9). 


The  quantity  dG^/dx^  is  found  directly  from  Equation  (7)  as 

d0 .  I-  (xUx!, 

-  tQ^yx2)  -  x^y  si n0—  • 


(16) 


This  equation  can  be  solved  immediately  (by  numerical  methods)  sub¬ 
ject  to  the  appropriate  boundary  conditions  (see  Figure  1  and  Reference  8 
for  example).  If  I^e^sinQj  and  (x^)  x2  are  Integra ble  then  we 
can  solve  (16)  very  easily  as 


6, 
r 1 


I1(e1)  sinGjdej  =  P^Qj)  = 


e 


I3p(x2)  x2  dx2  =  P3a(x2)*  (17) 


1  min 


2  min 


This  integral  solution  will  be  of  especial  value  if  we  can  find  ©^(x!,) 
explicitly  from  (17).  Otherwise  it  appears  just  as  well  to  solve  (16) 
since  the  numerical  solution  can  be  used  simultaneously  to  solve  (15). 
That  is,  the  solutions  of  (14)  and  (15)  can  be  set  up  on  the  same  computor 


program^ 
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Before  finding  dy^/dx^  explicitly  we  will  find  dx^/dx^  as  a  function 

of  (0j,  yj,  dyj/dx^Jx^K  i *  Yl ) *  X2^1^1’X2^’  and  y2^1’yl’X2^  explicitly. 
Upon  finding  dyj/dx^  we  will  then  have  obtained  a  complete  solution  for 

the  system. 


From  the  diagrams  of  reflector  #1  it  is  seen  directly  that 

X1  =  W  yl>  =  [^^9l)"yll  tan0i* 

To  find  dxj/dx^  ^9 y^,  dy^/dx^;  we  use 


dx1  Qxj  d9j  9xj  dyj 

=  n5_0^  dx^  +  dY1  dx2 


with 


0x 


-  =  fie  tan  0X  +  (^-y^  sec2  9X 


whe  r  e 


=  <^/d01(  and 

9  x. 


dY  i 


=  -  tan  9 


1' 


We  can  find  y?(0l '  y  l  :x;9  and  ^(O^.y^x^)  explicitly  by  solving 
Equations  (9)  and  (13)  simultaneously.  This  gives 


y2  =  y2(e 


M  [(xi  -xz)2  +  (a+^-yi )2  “ 
i*yrx2J=  '  - 


(CL^fi-  yr) 


(18) 


(19) 


(20) 


(21) 


(22) 


and 
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with 


(23) 


E  =  Cp  -  sec  0j  -  C  {-Ojj. 

Now  from  Equations(3)  and  (4)  we  obtain 


(24) 


dy  dx 

^T-  =  ^  (0ryrx2)  W 


(25) 


whe  re 


fjtOj.  yx;x^)  =  tan 


/  x2  "  X1 

9.  -  arc  tan  ^  m - 

1  V  0L+/^+  y2  - 


(26) 


With  Equations  (18)  -  (24)  substituted  into  Equations  (25)  and  (26)  we 
obtain  explicitly 


dyl 

H5dT  =fy1(0l’yl;x2) 


|50itan0 


+  (fi-Yi)  se°Z ei 


f0l  ^ 


( 1  +  tan 0 ^  f . ) 


(27) 


Upon  solving  (27)  for  y^(x!>)  and  either  (16)  or  (17)  for  0j(x!,), 
then  all  the  remaining  dependent  variables,  Qy,  xi>  anc^  x2  can 

be  found  also  as  functions  of  x!>  by  application  of  Equations  (4), 
(18),  (22),  and  (23)  respectively. 
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III.  Existence,  Bounds  and  Uniqueness 


For  the  purpose  of  determining  existence  and  uniqueness  of  the  solution 
and  obtaining  bounds  for  and  ,  it  is  convenient  to  use  the  form 
of  the  equations  given  in  (3)  and  (5)  (see  also  Reference  9).  For  any 
specific  synthesis  problem,  it  is  frequently  possible  to  estimate  bounds 
for  the  solutions  y^(x^)  anc^  y^x^)  *n  ac^ition  to  predicting  uniqueness 
merely  by  using  the  forms  given  in  (3)  and  (5).  The  methods  for  employ¬ 
ing  this  useful  property  are  easily  seen  from  a  study  of  the  theory  of 
these  equations  such  as  given  in  Reference  8.  Some  discussion  of  this 
is  given  in  Appendix  I. 


IV.  Uniform  Phase  and  Amplitude  Examples 


The  solution  for  the  dual-reflector  system  which  will  produce  a  uniform 
phase  and  amplitude  distribution  in  the  aperture  of  reflector  #2  is  found  for 
a  primary  illumination  1^(0^)  =  cosn0^  It  is  assumed  that  the  primary 
phase  front  is  circular  with  center  location  given  by  the  constants.  cC  and 
.  In  this  case  we  will  take  the  phase  center  of  primary  radiation  as  the 
phase  reference  for  the  system  so  that  C  (0j)  =  0,  The  aperture  ampli¬ 


tude  distribution  will  be  given  by  IjpC*-?)  “  constant  found  from  (8).  The 
aperture  phase  distribution  C^(x^)  is  also  given  as  a  constant  which  is 
determined  from  (9)  and  the  boundary  conditions.  We  may  note  immediate! 
that  in  this  case 


=  0 


and 


(28) 
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The  boundary  conditions  will  be  chosen  so  that 


0  —  x  ’ 

1  min  2  min 


=  x 


2  min 


=  0. 


Four  more  independent  boundary  conditions  are  available.  We  will 
choose 


c,  .x,  ,  0.  » 

2  max  1  max  1  max 


and  CC . 


The  remaining  conditions  are  found  f  rom  Figure  1  as 

S 


'i 


max 


tan  0 


1  max 


< 


6,  » arc  tan 

2  max 


fa 


C  (xL)  =(  - H- 

p  '  Z  \  cos  6, 
r  V  1  max 


x2max  ~  X1  max 

aTT~ jS 

_ \  +  (  au+  fa 


cos  0 


2  max 


Since  the  overall  scale  of  the  structure  is  arbitrary,  the  value  of 


x-  will  be  chosen  as  either 

2  max 


(29) 


a. 


c0  ^  =  +  1,  which  leads  to  a  convex  type  of  solution  for 

£max  *  * 


reflector  #1 


or 


x2max  =  "  wkich  l^ads  to  a  concave  type  of  solution  for 
reflector  #1. 


In  this  example  we  find  that  Equation  (17)  is  integrable  and  that 
0  j(x^)  =  ® i ( x 2 )  can  found  explicitly  as 
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0j  =  arc  cos 


/n+1  I 

W  1  -  M  x, 

\  c  ) 


with 


w  v  n+1  _ 

M  =  1  -  cos  0 


1  max 


(30) 


(31) 


From  Equation  (22),  with  dC  /dx'  =  sin  0-  =  0  and  C  (0.)  =  0,  we 

P  c  *  Pq  1 


obtain 


y2 =  y2(9i*  yi  ;x2^  = 


where  rj  is  given  by  (10). 


(Cp"  rl)  "  (x2  "  xl)  ~  yl) 

2«C+^-  Yj  +  c  -  rj) 


It  is  now  possible  to  construct  the  function  (Equation  27) 


dYi 


=  fyi(ei,y>sx2)  =  fy1(yl;x2) 


(33) 


explicitly  and  solve  for  and  x^  as  functions  of  x 


It  should  be  mentioned  that  the  alternative  solution  discussed  previously 
in  terms  of  dy1/dx1  =  f^Xj,  y^  and  dy2/dx2  =  f2(x2>  y  )  leads  in  the 
present  example  to  explicit  equations  for  fj(Xj,  y^)  and  f2(x2>  Y2) 

(see  Reference  9). 

Numerical  computations  were  made  for  the  following  special  case: 


n  =  16,  x,  =  0.  1  ,flC=  0,  and  ft,  =  30' 

i  max  1  max 


The  values  chosen  are  reasonably  practical; 
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The  results  for  the  convex  (Cassegrain  type)  solution  are  depicted  in 
Figure  5.  A  parabola,  matched  at  two  points,  is  drawn  for  comparison 
with  reflector  #2.  The  results  for  the  concave  (Gregorian  type; 
solution  are  depicted  in  Figure  6.  In  both  cases,  a  numerical  and 
graphical  check  of  the  results  indicate  that  the  solutions  conform  to 
the  optical  constaints  as  required.  The  smoothness  of  the  contours 
(large  radii  of  curvature)  indicate  that  the  optical  design  will  be  a  useful 
one  in  this  case. 
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CHAPTER  2 

THE  SYNTHESIS  PROBLEM  AND  THE  L  UN  EB  ERG  -  KLIN  E  EXPANSION 

I.  Introduction 

The  synthesis  method  developed  in  Chapter  1  was  based  on  the  principles 
of  geometrical  optics.  Although  these  principles  were  discovered 
independent  of  electromagnetic  theory,  they  may  be  derived  from 
Maxwell's  equations.  For  the  design  of  optical  instruments  which 
operate  at  wavelengths  measured  in  angstroms,  the  connection  with 
Maxwell's  equations  if  often  unimportant.  At  radio  and  microwave 
frequencies  where  the  wavelength  becomes  significantly  large,  the 
connection  of  geometrical  optics  to  Maxwell's  equations  is  important, 
for  it  allows  us  to  determine  how  accurate  the  geometrical  optics 
approximation  is  and  how,  if  possible,  geometrical  optics  designs  may 
be  improved. 

For  example,  the  geometrical  optics  field  scattered  from  the  subreflector 
of  our  dual  reflector  system  will  have,  in  general,  a  sharp  cutoff  at 
the  shadow  boundary  as  is  illustrated  below  (a).  For  subreflectors 
many  wavelenths  in  diameter,  the  scattered  field  observed  appears  to 
have  just  such  a  sharp  discontinuity.  When  the  subreflector  is  of  the  order 
of  7  or  8  wavelengths  and  less,  significant  changes  in  the  true  scattered 
field  from  the  geometrical  optics  field  may  be  noticed.  Diffraction 
effects  become  significant  and  one  observes  a  scattered  field  such  as  is 
illustrated  below  in  (b).  In  addition  to  the  ripple  that  is  noticeable,  a 
large  change  in  the  amplitude  is  observable  at  the  shadow  boundary. 

A  detailed  discussion  and  analysis  of  this  diffraction  effect  will  be  made 
in  later  chapters. 


Another  consideration  that  is  necessary  at  larger  wavelengths  is  the 
effect  of  a  small  radius  of  curvature  for  a  reflecting  surface.  That 
is,  we  would  like  to  know  what  radius  of  curvature,  in  wavelengths, 
is  necessary  for  geometrical  optics  to  be  a  valid  approximation  to 
electromagnetic  theory.  Although  this  problem  cannot  be  solved 
exactly  for  most  practical  cases,,  considerable  insight  into  the 
effects  of  diffraction  and  small  radii  of  curvature  can  be  obtained 
by  deriving  the  principles  of  geometrical  optics  from  Maxwell's 
equations.  There  have  been  many  such  derivations  ’  ,,  but  the  most 
satisfactory  appears  to  be  that  developed  by  Luneberg  and  Kline^ 
in  the  form  of  an  asymptotic  solution  to  Maxwell's  equations.  This 
expansion  in  inverse  powers  of  k(2Tr/x  where  >«=  wavelength)  is  shown 
to  have  as  its  first  term  the  geometrical  optics  field.  Hence 
geometrical  optics  may  be  referred  to  as  a  first  order  solution  to 
Maxwell's  equations. 

Strictly  speaking,  under  certain  conditions,  the  higher  order  terms  of 
the  Luneberg-Kline  expansion  satisfy  the  geometrical  optics  principles 
also.  The  important  condition  is  the  usual  far  field  assumption  that 
the  fields  are  considered  only  at  a  great  distance  from  their  sources. 
Hence  the  primary  source  field  in  the  synthesis  problem  may  be 
represented  exactly  by  geometrical  optics  if  the  far  field  is  con¬ 
sidered  or  if  the  subreflector  lies  in  the  far  zone  of  the  primary 
source . 

An  additional  reason  for  discussing  the  Luneberg-Kline  expansion  is  its 
similarity  to  the  asymptotic  expansion  of  the  Sommerfeld  solution  to 
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the  half-plane  problem.  The  higher  order  terms  of  the  asymptotic 
solution  to  the  half -plane  problem  describe  the  sharp  edge  diffracted 
field.  The  second  term  provides  an  accurate  description  of  the  far 
zone  diffracted  field.  One  distinction  between  the  Luneberg-Kline 
expansion  and  the  asymptotic  expansion  of  the  half-plane  solution  is 

A 

that  the  latter  is  in  powers  of  k  1  . 

In  general  it  is  not  possible  to  obtain  an  asymptotic  solution  to  a 

diffraction  problem  in  powers  of  k  ^  where  ^  is  an  integer  unless 

3 

certain  conditions  are  satisfied.  The  original  paper  by  Kline  derives 
and  explains  these  conditions  quantitatively.  In  essence,  these  condi¬ 
tions  require  the  diffracting  obstacle  to  have  a  convex  curvature  and 
to  cast  no  optical  shadows  when  illuminated  by  the  incident  field.  A 
surface  that  is  concave  at  any  point  may  cause  reflected  rays  to  cross  at 
a  caustic.  The  asymptotic  expansions  do  not  exist  on  the  caustics  of  a 

field  defined  by  a  ray  system,  A  particular  expansion  may  however  be 

1 4 

valid  away  from  a  caustic.  It  can  be  shown  that  more  terms  of  the 
expansion  are  generally  necessary,  however,  to  evaluate  the  field  nearer 
to  the  caustic.  Away  from  the  caustics,  then,  the  field  diffracted  from 
a  surface  which  casts  no  shadow  boundaries  may  be  expressed  as  an 
asymptotic  expansion  in  k~  ..  Such  surfaces  may  be  infinite  paraboloids 
and  hyperboloids 


When  the  diffracting  surface  casts  a  shadow  boundary  it  is  still  usually 
possible  to  express  the  solution  in  terms  of  an  asymptotic  expansion 

with  the  distinction  that  non-integer  powers  of  k  must  frequently  be 
used  As  was  mentioned  earlier,  the  exact  half-plane  solution 
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$ 

may  be  expanded  in  powers  k  where  is  an  integer.  Unfortunately, 
in  most  of  these  cases  the  expansions  are  known  only  after  the  exact 
solution  is  found..  Nevertheless,,  it  should  be  noted  that  all  the  expan¬ 
sions  have  the  same  first  term,  the  geometrical  optics  solution.  Hence 
the  field  diffracted  from  the  reflectors  considered  in  this  work  will 
always  consist  in  part  of  the  surface  reflected  field  given  by  geometrical 
optic  s. 


The  Sommerfeld  half-plane  solution  suggests  that  a  substantial  part  of 
the  diffracted  field  in  our  problem  is  due  to  edge  currents  near  the 
knife  edge  of  the  reflector.  These  currents,  as  will  be  shown  in  the 
next  chapter,  give  rise  to  a  field  which  is  of  higher  order  in  powers  of 
1/k  than  the  geometrical  optics  field.  In  problems  for  which  the 
solution  is  not  known  apriori  it  has  been  possible  to  obtain  excellent 


results  for  the  diffracted  field  by  clever  use  of  the  edge  current 
791011 


concept 


Application  of  this  app.  eximation  will  be  made 


for  the  subreflector  diffracted  field  in  the  next  chapter. 


When  the  reflecting  surface  is  curved  two  additional  considerations  are 
important.  The  first  concerns  the  possible  excitation  of  a  surface 

5 

diffracted  wave  by  rays  incident  tangentially  to  the  surface  .  In  the 
subreflector  problem  this  is  possible  due  to  secondary  edge  diffracted 
rays  enamating  from  the  edge  current.  Their  effect  is  estimated  and 
shown  to  be  small  in  the  next  chapter.  An  additional  consideration  is 
the  accuracy  of  the  geometrical  optics  field  when  the  reflecting  surface 
has  a  small  radius  of  curvature. 
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This  last  question  is  too  difficult  to  answer  in  a  precise  and  conclusive 
manner.  For  example,  it  has  been  shown  that  the  geometrical  optics 
solution  is  an  exact  solution  to  the  vector  electromagnetic  diffraction 
problem  of  a  plane  wave  incident  on  an  infinite  paraboloid  of  revolu¬ 
tion^,,  This  solution  is  exact  and  independent  of  the  paraboloid  curva¬ 
ture  at  the  apex.  This  result  is  unusual,  however,  and  in  fact,  pecu¬ 
liarly  enough,  does  not  apply  to  the  equivalent  acoustic  problem.  In 
general,  when  the  radius  of  curvature  of  the  reflecting  surface  is 
larger,  then  fewer  terms  of  the  expansions  in  powers  of  k  are 

needed.  This  has  been  shown  in  a  number  of  cases  investigated  by 

5 

Keller,  Lewis,  and  Seckler  ..  They  apply  the  Luneberg -Kline  expansion 
in  analyses  of  parabolas  and  hyperbolas  of  revolution  and  also  of 
cylindrical  shape.  The  solutions  for  these  cases  are  generally  difficult 
to  interpret  since  the  geometrical  parameters  describing  these  shapes 
(radius  of  curvature,  focal  length,  etc.)  do  not  enter  into  the  expansions 
in  a  reasonably  simple  way.  The  application  of  the  expansion  to  reflec¬ 
tion  from  a  circular  cylinder  allows  a  direct  interpretation  of  the  effect 
of  surface  curvature  on  the  convergence  of  the  expansion.  As  discussed 
earlier,  the  application  to  this  problem  is  approximate  and  does  neglect 
diffraction  due  to  shadowing,  etc.  This  example  is  presented  after  the 
theory  is  developed  in  a  more  quantitative  manner  next. 
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II, 


The  Luneberg -Kline  Expansion  and  Geometrical  Optics 


We  will  assume  a  time  harmonic  0 


■  iaJt 


dependence  and  suppress 


5  6 

this  factor.  Our  notation  will  follow  those  of  Keller  and  Schensted 


which  are  the  same  if  we  set 


n 


n 


E  .  We  have  three 
n 


sets  of  waves;  the  primary  wave  incident  on  reflector  #1,  the 


—N. 


secondary  wave  incident  on  reflector  #2,  and  the  tertiary  wave 

_ \ 

passing  outward  through  the  aperture.  We  synthesize  for  a  given 

tangential  in  the  aperture  with  specified  arbitrarily.  The  vector 

-iWt 


U  is  the  electric  vector  with  an  0 


time  factor  suppressed. 


The  Luneberg-Kline  asymptotic  expansions  for  the  fields  are 


U, 


£  ^/ik)n 


£  tffcn/ (ik)n 


ikC 


P1 


£ikC 


c  .  -A.  +  c 

pl  1  po 


c  t  —  +  r* v  +  c  I 

pZ  2  1  po  (2.) 


(1) 


tf3  =  £/^/j/(ik)n  ^  P3  j  Cp3  +^2  +t*j  +  CpJ(3) 


The  notation  is  the  same  as  that  used  previously  except  that^.  is  the 
distance  from  the  Y-axis  for  j  =  1,  the  distance  from  reflector  #1 
for  j  =  2,  and  the  distance  from  reflector  #2  for  j  =  3,  in  each  case 
along  the  ray  which  is  normal  to  the  optical  phase  fronts  defined  by 


C  .  = 
PJ 


constant  • 


(4) 


In  order  to  have  (1),  (2),  and  (3)  satisfy  Maxwell's  equations  we  first 
insert  them  into  the  homogeneous  vector  wave  equation 

(V2  +  kz)lj.  =0 

and  set  to  zero  the  coefficients  of  like  powers  of  k”  ^  .  This  leads 
to  a  recurrence  relation 


'  4"  V^n-l 


(5) 


where 


The  divergence  condition 


0)  gives 


(6) 


The  boundary  condition  on  a  perfectly  conducting  reflector  surface 
defined  by  the  unit  normal  fl  gives 


Y\  x  /%  =  o 


Finally,  from  the  wave  equation  we  also  obtain  the  eikonal  equation  of 
geometrical  optics 


(The  subscripts  j  have  been  dropped  from  (5)  through  (8)).  Hence  (1) 
through  (3)  satisfy 
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V2  +  k2 
V-U 


u 


0 


(9) 


and  H  *  U  =  0 


if  (5)  through  (8)  are  satisfied. 


Equation  (8)  implies  that  there  is  a  system  of  rays  which  form  orthogonal 
trajectories  to  the  wavefronts  =  constant.  This  wavefront  is  a  true 

wavefront  for  the  first  term,  AC  ,  of  the  expansion.  It  is  generally 
not  a  true  phase  or  wave  front  for  the  field  when  more  terms  are  used. 


Equation  (5)  is  a  first  order  ordinary  differential  eqnvtion  in  the  orthog¬ 
onal  coordinate  system  defined  by  the  ray  trajectories  and  phase  fronts. 
This  equation  permits  us  to  relate  the  field  at  one  point  of  a  ray, 

to  the  field  at  another  point,  ,  when  the  behavior 

between  and  is  known  for  the  term.  This  equation  describes 

Z  a 

the  flow  of  energy  along  a  tube  of  rays  .  The  Tfc-1  term  is  a 
measure  of  the  diffusion  of  energy  between  adjacent  tubes  of  rays. 

Energy  flow  is  confined  within  a  tube  of  rays  when  = 

This  is  true  for  the  /l  =  0  or  geometrical  optics  term. 

Finally,  equation  (6)  is  of  very  great  significance  insofar  as  it  helps  to 
describe  the  vector  nature  of  the  field.  We  may  write  (6)  as 


VV  <  --  AK.-  -y-AK 


•i . 


That  is,  the  longitudal  component  of  along  /i  ,  is  given  by 


-V  *  .  • 

fir*  / 


Hence  /££  is  always  normal  to  the  ray  since  ~  O'. 


At  first  order,  f\  -  0,  equations  (5)  through  (7)  become 


“34:  +  (7  cp )  -  ° 

(10) 

Vcp  -  0. 

(ll) 

ft  x7\r0  -  0. 

(12) 

Equations  (8),  (10),  (11),  and  (12)  are  the  usual  geometrical  optics 


laws . 
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III.  Example  of  the  Effect  of  Surface  Curvature  on  the  Accuracy  of 
Ge  om  e t r i c aT  Optic  s  Refl e ction  " 

Geometrical  optics  principles  clearly  do  not  account  for  the  diffraction 
effects  at  a  knife  edge.  One  may  consider  the  radius  of  curvature  at 
such  an  edge  as  being  zero  (or  the  Gaussian  curvature  as  being  infinite). 
It  is  well  known  that  geometrical  optics  does  give  good  results  for  the 
fields  scattered  from  large  flat  or  nearly  flat  reflectors.  Hence  one 
may  state  that  if  the  radius  of  curvature  of  the  reflecting  surface  is 
large  then  geometrical  optics  or  the  first  term  of  the  asymptotic 
expansion  of  the  field  gives  an  accurate  result  for  the  field.  In  this 
section  we  wish  to  give  some  insight  into  the  problem  of  just  how  large 
in  wavelengths  the  radius  of  curvature  of  the  reflector  surface  must 
be  in  order  for  geometrical  optics  to  give  accurate  results. 


In  the  introduction  to  this  chapter,  it  was  mentioned  that  a  specific  answer 

to  this  question  cannot  be  given  since,  for  example,  the  geometrical 

5  6 

optics  solution  is  exact  for  the  infinite  paraboloid  9  .  a  number  of 

5 

different  problems  have  been  analyzed  by  Keller,  Lewis,  and  Seckler 

by  use  of  the  Luneb erg -Kline  expansion  wherein  the  convergence  properties 

-  77. 

of  the  k  series  is  illustrated.  In  many  of  their  problems  the  k 
expansion  is  in  fact  a  true  asymptotic  solution  since  no  shadow  boundaries, 
caustics,  etc.  ,  exist  in  the  problem.  The  convergence  dependence  on 
curvature  is  best  illustrated  however  by  applying  the  Luneber g-Kline 
expansion  to  scattering  from  a  metal  circular  cylinder.  The  application 
is  approximate  since  diffraction  effects  exist  here  and  are  neglected; 
nevertheless,  the  problem  is  illustrative  since  the  radius  of  curvature 


of  the  scatterer  is  a  constant  and  the  convergence  dependence  on 
this  radius  is  easily  discernible.  Furthermore,  away  from  the 
shadow  boundary  and  shadow  region  the  expansion  is  correct. 
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For  the  sake  of  brevity,  the  complete  details  of  the  derivation  of  this 
solution  are  not  presented  here.  A  complete  derivation  is  given  by 
Keller  et  al  .  In  essence,  the  geometrical  optics  solution  is  found 
by  direct  application  of  the  stated  principles  and  the  higher  order 
terms  are  then  found  from  the  geometrical  optics  solution  by  applica¬ 
tion  of  the  recurrence  relation  discussed  earlier. 


A  plane  wave,  U.  js  incident  upon  a  perfectly  conducting  right 

circular  cylinder  from  the  left  as  illustrated  below  (c).  The  cylinder 

radius  is  bQ  «  The  electric  field  is  oriented  parallel  to  the  cylinder 

axis  so  that  the  total  field  E  =  U.  =  0  on  the  cylinder.  The 

z  t 

incident  field  rays  are  parallel  to  the  x-axis.  The  reflected  wavefront 
is  not  circular  (except,  as  always,  it  is  very  nearly  so  in  the  far  zone). 
The  reflected  rays  and  phase  fronts  form  an  orthogonal  coordinate  system. 
The  rays  are  tangential  to  the  caustic  curve  which  is  defined  by 


sm 


(  (Si  2  ) 


(13) 


where 


is  the  distance  from  the  caustic  to  the  reflector  surface. 


The  solution  for  the  scattered  field  is  obtained  in  a  non -orthogonal  coor¬ 


dinate  system  wherein  a  point  is  defined  by 


(a4  ).  The  distance 
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i“  that  taken  along  a  reflected  ray  from  the  caustic  to  the  point  of 
observation  P.  The  angle  is  that  made  by  the  reflected  ray  and 

positive  X-axis.  The  solution  for  the  reflected  or  scattered  field  is 
the  Luneberg -Kline  expansion 


(14) 


The  expansion  is  found  by  application  of  the  recursion  formula  to  the 
geometrical  optics  solution.  The  coefficients  are  independent 

of  b  ,  >4  ,  and  /$>  and  are  given  by  the  recursion  relations 


*  j'1  [<2i  +  «  +  -  4t  -  2"  - 

+  (2j  -  4t  +  2n  +  5)(2j  -  4t  -  2n  +  3 t-l.n-1 
+  [  24(j -l)(j -2t-n)  -  6^.2,  t>n., 

(15) 

+  12(1  -  j ) ( 2 j  -  4t  -  2n  +3)  t  t_i ,  n.i 

*  9(2j  -5)0  -V)a}.3,t,n.l 


+  9Uj-5)(2J- 


and 
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and 


£^otn  "  “  Z,  0^  jtn 

j-' 


fl  -  -2 

\Aooo 


(16) 


(17) 


First  of  all,  it  should  be  noticed  that  the  asymptotic  solution  (14) 
blows  up  at  the  shadow  boundary  ^2  =  0  as  may  be  expected.  Away 

from  the  shadow  boundary,  £  >  o,  we  wish  to  investigate  the  conver¬ 
gence  properties  of  the  expansion  with  respect  to  the  radius  of  curvature 
bQ.  Except  for  the  factor,  bQ  occurs  only  together  with  k 

as  a  factor  (l/bQk)n  ,  In  general  we  are  interested  in  the  field  far 
away  from  the  reflecting  surface  v-here  b^ .  In  that  circumstance 

we  may  consider  only  the  far  field  solution  j  =  0  and  make  the  approxi¬ 
mation 

e ik  u ' 1  sin^  *  It  a 


n-o  -^o 


‘otn 


(18) 


>in^  )  (sin  $ ) 


-2t  . 


(16ibQk  sin'c  )  (sin  /2T  ) 


In  this  important  case  we  note  that  the  geometrical  optics  solution  in 
the  far  field, 


U  /v/ 
o 


3 


sin/^  ^ 


ik(/4,-  *2*  sin  c~ ) 


■  & 


(19) 


is  accurate  uniformly  in  <0  0>^O)  so  long  as  (k  b  )  is  large.  The 

'  r  r  °  y*  n  o-i 


percentage  error  in  this  specific  case  is 


100/(1 6b 


oksin'4>| 
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when  the  effect  of  only  the  second  term  is  considered.  Note  that  the 
error  increases  nearer  to  the  shadow  boundary  /=0. 

This  example  illustrates,  quantitatively  for  this  case  only,  that  the 
geometrical  optics  solution  is  more  accurate  for  reflecting  obstacles 
with  larger  radii  of  curvature.  As  a  crude  measure  of  the  accuracy 
of  the  geometrical  optics  solution  to  an  arbitrary  reflection  problem 
we  may  use  the  factor  (1/kb^)  as  the  approximate  order  of  magnitude 
of  the  second  term  in  the  asumptotic  expansion.  We  may  regard  b 
as  the  local  radius  of  curva  ture  of  the  reflecting  surface  and  the 
accuracy  thus  obtained  will,  roughly,  apply  to  the  field  reflected 
from  that  portion  of  the  surface.  This  confirms  what  is  known  generally 
through  experience.  It  is  expected  that  k  and  b^  should  occur  together 
since  any  other  result  would  contradict  the  fact  that  it  is  the  size  of 
bQ  as  measured  in  units  of  wavelength  only  that  is  significant.  That 
result  may  be  expected  generally. 
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CHAPTER  3 

THE  EDGE  DIFFRACTED  FIELD 


Introduction 

A  consideration  of  great  importance  in  any  dual  reflector  system 
is.  the  extent  and  effect  of  the  diffracted  field  scattered  from  the 
smaller  or  subreflector.  The  geometrical  optics  field  scattered 
from  the  si^breflector  is  the  field  that  is  considered  in  the  synthesis 
method  developed  in  Chapter  1.  From  the  discussion  in  Chapter  2, 
it  is  clear  that  this  field  is  first  order  with  respect  to  the  Luneberg- 
Kline  series  field  description,  The  remaining  terms  or  the  diffracted 
field  scattered  from  this  reflector  may  or  may  not  be  of  great  con¬ 
sequence  depending  upon  the  radius  of  curvature  and  overall  size 
of  this  reflector  in  wavelengths.  A  smaller  s ub- reflector  will  have  a 
scattered  field  which  is  more  different  from  the  geometrical  optics 
field.  This  is  expected  from  experience  and  is  borne  out  in  the  numer¬ 
ical  results  presented  later. 

In  general  it  is  important  to  have  as  small  a  sub-reflector  as  possible. 
Besides  the  obvious  mechanical  reasons  of  weight  support,  etc., 

(which  are  of  very  great  practical  significance)  the  blockage  of  the 
main  reflector  by  the  smaller  one  reduces  the  overall  gain  and  raises 
the  sidelobes  of  the  tertiary  pattern.  The  diffracted  field  of  the  larger 
main  reflector  is  generally  negligible  in  the  vicinity  of  the  aperture  of 
this  reflector.  'Ihat  is,  geometrical  optics  is  very  accurate  for  a  very 
large  reflector  The  far  field  of  the  main  reflector  aperture  distribu¬ 
tion  is  usually  determined  by  a  Kirchhoff  *  integration  over  the  aperture 
which  accounts  for  diffraction  effects  with  reasonable  accuracy.  For 
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a  uniform  aperture  phase  distribution  the  far  zone  field  is  a  purely 
diffraction  determined  field. 

We  will  use  the  term  uKirchho££  integration"  rather  frequently  in 
what  follows  so  that  a  brief  definition  of  this  term  is  in  order.  The 
definition  given  by  Baker  and  Copson^  is  interpreted  broadly  here. 

Their  definition  applies  more  directly  to  diffraction  through  a  hole 
in  a  screen  rather  than  from  a  metal  scatterer.  If  E^,  is  the  field 
incident  upon  the  aperture  and  generated  by  the  sources  in  free  space 
(no  screen  or  scatLerers  present),  then  a  Kirchhoff  integration  would 
assume  that  E^,  H-  is  the  total  field  in  the  aperture,  and  the  tangential 
field  on  the  screen  is  zero.  Strictly  speaking,  neither  assumption  is 
true.  For  scattering  from  a  metal  scatterer,  we  will  interpret  a 
Kirchhoff  integration  as  follows: 

If  E^,  H-  is  the  field  incident  upon  the  scatterer,  we  assume  that  the 
current  distribution  on  the  scatterer  surface  is  given  by  2  ft  x  H- 
where  n  is  normal  to  the  surface.  Those  portions  of  the  surface 
which  are  shadowed  by  the  scatterer  are  assumed  to  have  zero  current. 
The  assumptions  made  here  are  similar  to  those  made  for  the  hole  in 
a  screen  problem  and  are  similarly  approximate.  An  integration  of 
these  approximate  surface  currents  to  obtain  the  scattered  field  will 
be  termed  a  Kirchhoff  integration 

The  purpose  of  this  chapter  will  be  first  to  determine  the  diffraction 
field  of  the  subreflector.  In  addition  we  will  determine  the  actual  error 
in  aperture  phase  and  amplitude  caused  by  this  second  order  effect..  In 
subsequent  chapters,  we  will  discuss  various  efforts  to  correct  the 
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diffraction  field  effects  in  the  main  reflector  aperture  and  numerical 
results  will  be  presented. 

An  exact  determination  of  the  field  scattered  from  reflector  shapes 
that  are  used  in  practice  is  not  possible  in  general  except  for  the  flat 
circular  disk^.  A  Kirchhoff  integration  over  the  subreflector  surface 
such  as  that  made  by  Rusch  has  several  disadvantages.  The  accuracy 
for  very  small  reflectors  (less  than  4  A  in  diameter)  is  questionable. 

The  analysis  requires  laborious  numerical  computations.  The  entire 
reflector  shape  is  required.  This  is  so  since  this  analysis  yields  the 
first  order  geometrical  optics  field  and  the  diffracted  field  in  the  same 
computation. 

The  diffraction  field  is  primarily  an  edge  phenomenon  caused  and  deter¬ 
mined  by  the  abrupt  discontinuity  of  the  reflector  surface  and  the  actual 
edge  shape  (knife  edge,  rolled  edge,  wedge  or  square,  etc.).  Hence 
many  very  successful  diffraction  theories  have  been  developed  which 

depend  upon  only  the  shape,  position,  and  slope  of  the  edge  of  the 

4  5 

scatterer.  Although  such  edge  shapes  as  the  wedge  *  or  cylindr ic ally 
tipped  edge^  could  be  treated  we  will  confine  ourselves  to  the  most  sam¬ 
ple  subreflector  edge,  the  knife  edge.  In  particular,  we  will  apply  the 

7 

geometrical  theory  of  diffraction  as  developed  by  Keller  .  As  discuss¬ 
ed  in  Chapter  Z,  this  theory  is  based  on  the  Sommerfeld  half-plane 
solution^.  Other  methods^*  ^  which  are  similar  have  been  based 
on  the  same  solution  and  have  yielded  very  satisfactory  results  as 
determined  experimentally  1  ® 9  The  geometrical  theory  of  diffraction' 

is  somewhat  simpler  to  apply,  and  yields  very  accurate  results  for 
reflectors  of  the  order  of  one  wavelength  or  greater  in  diameter. 
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It  should  be  emphasized  that  the  exactness  of  the  solutions  for  the 
diffracted  field  found  in  this  chapter  is  of  no  great  importance  to  the 
general  conclusions  which  may  be  drawn  from  the  numerical  results 
presented  later  in  this  chapter,,  The  edge  diffracted  fields  found  in 
this  chapter  .will,  however,  be  very  accurate  in  general  If  the  primary 
source  amplitude  and  phase  does  not  vary  rapidly  near  the  edge.  The 
edge  diffraction  theory  is  based  on  an  assumption  that  locally,  very 
close  to  the  edge,  the  physical  situation  resembles  the  half-plane 
problem  of  Sommerfeld, 

The  edge  diffracted  field  has  a  non-circular  phase  front.  Opposite 
edges  of  the  subreflector  as  shown  in  Figure  1  are  illuminated  sym¬ 
metrically  by  the  primary  source.  According  to  the  theory  of  geomet¬ 
rical  diffraction  each  edge  contributes  to  the  total  field  only  In  the  plane 
normal  to  the  edges  and  containing  the  axis  of  the  system.  Since  each 

edge  is  a  directional  source  oriented  in  different  directions,  there  can 

13 

be  no  fixed  phase  center  for  the  combination  ,  An  equivalent  geomet¬ 
rical  optics  description  of  this  field  will  therefore  require  the  generality 
used  in  the  description  of  the  primary  source  in  Chapter  1, 

The  descriptive  geometrical  parameters  are  illustrated  in  Figure  Z. 

Since  it  will  be  convenient  later  to  add  the  primary  and  edge  fields  at 

the  aperture,  the  particular  edge  field  ray  intersecting  the  primary 

field  ray  at  the  aperture  is  illustrated  in  the  Figure,  Je(  ^ is  the 

angular  power  density  of  the  diffracted  field  in  the  same  units  as  1^(0^)* 

I~  (Xi)  is  the  aperture  power  flow  density  normal  to  the  aperture,,  The 
o  e  z 

edge  field  ray  intersects  the  y  axes  at  a  distance  de  (  (V^)  from  the  _ 
aperture.  Not  illustrated  in  Figure  2  is  the  reference  wavefront 
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Ce  (  $9)  which  serves  the  same  purpose  for  this  field  description 
o  L 

that  Cpo(0,)  served  for  the  description  of  the  primary  source  field, 

The  quantity  Ce(  ^ is  defined  as  the  aperture  phase  given  by 

+  Ce  .  Our  objective  in  this  chapter  is  to  find  the  values 
o 

of  these  parameters  for  various  frequencies  and  reflector  systems. 


II.  Edge  Diffracted  Fields  and  Sources 

We  propose  to  define  precisely  what  is  meant  by  an  edge  diffracted 
field  or  an  edge  source  in  this  section.  This  will  be  done  by  describ¬ 
ing  the  Sommerfeld  solution  to  the  half-plane  problem  in  detail. 

Consider  the  diagram  on  the  following  page  which  is  a  cross-section 
of  the  semi-infinite  or  half-plane  located  at  4>  ~  0.  We  assume  a 
plane  wave  is  incident  on  the  half-plane  at  an  angle  .  The  incidence 
is  normal  in  the  sense  that  there  is  no  z -dependence,  the  plane  wave 
travels  in  the  plane  of  the  diagram.,  The  polarization  is  either  normal 
or  parallel  to  the  edge. 

The  diagram  is  divided  into  three  regions,  I.  II,  and  III  by  boundaries 
G~  and  G;  and  the  half -plane.  The  boundaries  G  and  G.  are  the  geo- 
metrical  optics  shadow  boundaries.  Region  I  contains  the  incident  and 
reflected  geometrical  optics  fields,  Region  II  contains  only  the  incident 
geometrical  optics  and  Region  III  is  a  shadow  region,  being  optically 

dark.  If  we  subtract  these  fields,  the  geometrical  optics  fields,  from 
the  total  field  solution,  there  will  remain  what  is  referred  to  as  the 
edge  diffracted  field.  If  we  examine  this  edge  diffracted  field  near  the 
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edge  and  on  the  surface  we  find  an  associated  current  distribution  on 
the  half -plane  that  decays  rapidly  away  from  the  edge.  In  fact,  in  the 
far  zone  of  the  edge  diffracted  field,  when  1/  /  kr  r  (r  -  distance 

from  the  edge),  the  energy  all  appears  to  flow  from  the  edge  itself. 

That  is,  the  edge  is  the  center  of  phase  of  the  far  edge  diffracted 
field.  Hence  Keller,  in  his  geometrical  theory  of  diffraction  is  able 
to  represent  this  diffraction  by  rays  emanating  from  the  edge  itself. 

We  must  note  that  the  geometrical  optics  representation  of  the  far 
zone  edge  diffracted  field  is  not  a  geometrical  optics  field.  The  geo¬ 
metrical  optics  field  is  that  first  term  of  the  asymptotic  expansion  in 
inverse  powers  of  k  where  k  is  allowed  to  approach  an  infinite  value. 

As  we  will  see  shortly,  the  far  field  term  pf  the  edge  field  is  a 
i//T7  term  (k  and  r  always  occur  together  in  this  dimensionless 
problem),  Hence  the  edge  field  becomes  zero  for  k-?oo,  except  at 
the  shadow  boundaries  where  the  asymptotic  expansion  does  not  con¬ 
verge,  However,  the  1  / 1/"~ k r  term  in  this  problem  is  the  first 

correction  to  the  geometrical  optics  field,  It  is  the  second  term  of 

n 

an  asymptotic  expansion  in  k  ^  of  the  exact  field,  the  first  term 
being  the  geometrical  optics  field. 

Letting  UQ  be  either  the  electric  or  magnetic  field  depending  upon 
whether  E  or  H  is  oriented  parallel  to  the  edge,  we  express  the 
incident  plane  wave  by 

-  ikrcos(cj)-  )  -  i  w  t 

Uq  O  { 2L-  s  uppre  ssed).  (;  1 ) 

8 

If  U  is  the  scattered  E  or  H  field,  as  appropriate,  then  the  solution 


can  be  expressed  asymptotically  (first  two  terms  only)  inside 
Regions  I,  II,  or  III,  as 
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It  is  clear  from  the  above  equations  that  the  asymptotic  solution  blows 

I 

up  at  the  shadow  boundaries  given  by  <J>  =  tt  ±  cjC-  .  The  negative  and 
plus  signs  are  used  accordingly  as 


when  U  -  E  (parallel  to  edge) 

Z 


when  U  =  H  (parallel  to  edge). 

z  ^ 


(5) 


If  we  subtract  the  geometrical  optics  field  from  the  above  equations, 
we  are  left  with  the  edge  diffracted  field  which  is  expressed  identically 
in  all  regions,  excluding  the  shadow  boundaries,  as 


Ue  ^ 


(1  +  i) 
4yfnkP 
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ikr 


sec 


(JU2£.)  *  sec  (i±f^-) 


(6) 


These  expressions,  (2)  through  (6),  are  valid  to  second  order  in  the 
asymptotic  expansion  of  the  exact  solution  for  the  diffracted  field 
and  for 


0  2  <J>  2it  exclusive  of  <J>  -  G  ,  .  (7) 

The  accuracy  of  these  equations,  for  a  given  (kr),  is  poorer  near 
4>  -  Gr  and  =  G^,  since,  of  course,  the  convergence  of  the 

asymptotic  expansion  is  more  slow  closer  to  Gr  and  G-.  Near  the 
shadow  boundaries  Gr  and  G{  we  must  use  the  exact  expression  for 
U.  We  will  give  this  expression  only  for  the  field  near  G  and  in 
Region  I,  since  this  is  the  region  which  concerns  us  in  the  synthesis 
problem.  With  8  defined  as  the  angle  away  from  cj>  -  Gr  in  a  negative 
<J>  direction,  we  have 


ue  (9<  Gr)  =  -  (  ±  1/2  e 


i  k  r  cos 


(8) 


i  u  C  f 2  sin  (  £/2) 

1  -  l)  _,ikr  coso  l  V  n  \  o  / 


_i)  ^ 

±  7  w 


.  ;  TT  s-t-'  £* 

>0  e  2  ~  dt 


uiii  eLkr 

4>/ir  k  r  '  sin  (<yC  +  cT/2) 


where  0  =  tt  -  -  4*  and  the  integral  in  (8)  is  the  complex 

14 

Fresnel  integral  which  is  well  ta1  ?lated. 


We  may  note  in  (6)  that  the  phase  center  of  this  edge  diffracted  field 


49 


is  at  r  a  0  (Le.  ,  the  phase  of  this  function  is  a  constant  independent 
of  9},  and  the  energy  radiated  appears  to  flow  from  this  point.. 

Now  in  order  to  obtain  an  approximate  solution  to  problems  more 
complex  than  that  of  the  half -plane,  many  other  analysts  have  used 
the  above  Sommerfeld  solution  as  the  basis  for  their  approximations. 

In  all  cases  they  assume  that  the  knife  edge  of  their  diffracting  obstacle 
is  locally  straight  and  that  the  conductor  is  locally  flat.  Braunbek  . 
is  analyzing  the  flat  circular  disk,  assumed  that  the  currents  near 
the  edge  of  the  disk  were  the  same  as  those  for  the  half-plane  (decaying 
as  (r)  away  from  the  edge)  and  integrated  the  contributions 

of  such  currents  along  the  edge  contour  to  obtain  the  diffracted  field. 
Millar^  and  Clemmow^  analyzed  the  diffraction  from  an  artibrarily 
shaped  aperture  in  a  plane  conducting  screen  when  a  plane  wave  is 
normally  incident  upon  the  screen  from  one  side.  They  assumed  that 
a  fictitious  filamentary  edge  current  exists  along  the  aperture  edge. 

This  current  is  assumed  to  radiate  the  same  field  pattern  as  the 
extended  current  that  Braunbek  assumed,  when  the  Braunbek  field  is 
evaluated  in  the  far  zone  |~ (k  r)  “  ^  I ^  <^I  .  Despite  this  latter  ass  ump¬ 

tion,  Millar  and  Clemmow  obtained  excellent  results  even  up  to  the 
aperture.  In  fact  their  results  are  excellent  even  for  apertures  on 
the  order  of  a  wavelength  in  diameter  (see  the  experimental  results 
of  Welch12). 

7 

Keller,  in  his  theory  of  geometrical  diffraction  from  knife  edges, 
also  uses  the  Sommerfeld  edge  diffraction  Solution.  However,  be  does 
not  sum  the  contributions  from  all  parts  of  the  edge.  The  edge  current 
at  a  given  point  on  the  edge  contributes  to  the  radiated  field  only  in 


the  plane  normal  to  the  edge  at  that  point.  This  method  assumes 
that  the  same  fictitious  current  filament  exists  along  the  edge  that 
Millar  and  Clemmow  assume. 

In  the  formulation  for  the  edge  source  which  will  follow,  we  will  use 

i. 

ihe  method  of  Keller.  His  far  field  results  for  the  slit  in  a  plane 
(analogous  to  our  cylindrical  problem)  and  also  for  the  circular 
aperture  in  a  conducting  plane  (analogous  to  our  axially  symmetric 
problem)  are  excellent  even  for  slit  widths  and  aperture  radii  of  the 
order  of  one  wavelength.  In  the  slit  problem,  th  5  edge  is,  of  course, 
straight  (as  in  our  cylindrical  problem).  In  Keller's  proble  ns,  as  in 
those  of  Millar  and  Welch,  the  edge  is  part  of  a  flat  conducii?  screen. 
Braunbek^,  however,  also  analyzed  the  problem  of  diffraction  by  a 
plane  wave  incident  upon  a  semi-infinite  truncated  funnel.  He  also 
used  the  Sommerfeld  half-plane  solution  as  the  basis  for  his  approxima 
tion  in  this  problem.,  In  this  case,  the  edge  is  part  of  a  curved  surface 
The  principle  effect  of  the  curved  surface  is  not  in  altering  the  basic 
nature  of  the  edge  currents,  but  in  altering  their  orientation  in  space 
and  in  altering  the  surface  reflected  geometrical  optics  field.  These 
changes  are  considered  in  the  analysis  made  here. 

An  additional  assumption  made  in  the  analysis  is  that  the  field  incident 
on  the  edge  is  locally  a  plane  wave.  This  assumption  is,  or  course, 
consistent  with  the  other  assumptions  that  the  edge  and  screen  are 
locally  straight  and  flat. 
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III.  The  Edge  Diffracted  Field 

Under  the  assumptions  of  symmetry  which  we  have  made,  the  present 
formulation  will  in  general  apply  for  either  the  cylindrical  or  rotation- 
ally  symmetric  system.  The  differences  are  minor,  so  we  will  dis¬ 
cuss  only  the  somewhat  more  practical  rotationally  or  axially  sym¬ 
metric  case, 

A  Combining  the  Edge  Fields  from  Opposite  Edges 

The  plane  of  incidence  of  the  field  incident  on  the  edge  is  normal 
to  the  plane  containing  the  tangent  to  the  edge.  In  this  plane 
we  have,  according  to  Keller,  radially  directed  rays  emanating 
from  the  point  caustic  at  the  edge  as  illustrated  on  the  following 
page.  The  ^-dependence  of  this  distribution  is  given  by  equa¬ 
tions  (6)  or  (8).  The  angle  q/C  can  be  determined,  since  the 
slope  of  the  reflector  at  the  edge  and  the  angle  of  incidence  of 
the  primary  source  are  given  boundary  conditions  of  the  synthesis 

problem:  0,  ,  and  ( )  =  tan  1/2  [~0  -  0  \ 

dxj  max  max  L  max 

There  are  two  edge  sources  radiating  fields  that  must  be  summed 

in  order  to  determine  J  (  /K),  Ce  (/fi),  and  de  These 

e  Ct  q  c 

sources  lie  at  opposite  edges  of  the  reflector  in  the  plane  contain¬ 
ing  the  axis  of  the  axially  symmetric  system.  Each  source  con¬ 
tributes  to  the  field  at  the  point  P  as  illustrated  in  the  diagram 
following  (c). 

The  field  radiated  from  edge  No.  1  will  be  designated  /\Te^  (r ^ 
This  field  is  constructed  of  three  factors  The  first 
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is  given  by  equations  (6),  (7)  and  (8)  except  that  for  the  axially 

-1/2 

symmetric  case  the  curvature  (kr)  in  formulas  (6)  through 

(8)  becomes  the  Gaussian  curvature  when  P  is  in  the  far  zone 

A 

field  of  the  edges  .  We  designate  this  factor  as  Ue  ^  (rj,  4>j). 

In  this  expression  we  will  write  -  C^-i  since  it  is  clear 

that  the  angle  of  incidence  is  the  same  for  both  edges.  The  field 
%  is  the  edge  field  for  an  incident  wave  of  unit  amplitude  and 
zero  phase  at  the  edge.  The  amplitude  and  phase  of  the  incident 
wave  at  both  edges  is  the  same  (polarization  neglected  for  the 
moment)  and  given  by  the  boundary  conditions  and  the  given  pri¬ 
mary  source  field.  We  find  the  amplitude  and  phase  of  the  wave 
incident  at  the  edges  as  follows: 


The  primary  source  electric  field  evaluated  on  reflector  No.  1 
is  given  by  (for  the  preferred  polarization  component) 

Ej  =  Gj  eik  [ri  +  cpo(0l)] 

where  is  a  factor  proportional  to  the  Gaussian  curvature. 
Consider  the  diagram  following  (d). 

At  the  reflector  edge  we  have 

Jlm  -  (0lm>  =  m 

The  curvature  of  the  wavefront  at  the  reflector  depends  on 


See  Appendix  II  or  Reference  7.  See  also  equation  (24b). 
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R0j  (0l)  and  rj  (0j).  From  the  diagram  we  find  that 

R0,  =  -(  —  —  -)  sin  9 ^ ,  and  RQ  =  -  (  -^-2f  )  sin  9,  .  (9a) 

1  d  Gj  lm  d  0i  m  rn 

The  radius  of  curvature  in  the  plane  of  the  above  diagram  is 
then 

R1  =  R0l  +  rj *  and  Rlm  =  Rolm  +  rlm  at  the  ed«e*  <9b> 

Hence  we  find  for  Gi  at  the  edge 

-  1  * 

Glm  *  ■J  <kRlm)  <k  nj 

The  phase  of  the  incident  primary  field  at  the  edges  is  given  by 

/oik  (ji  +  C  (0i  =  ,0 

'  max  po  *  'max'-*  d  (10) 

Hence  the  edge  diffracted  field  at  edge  No.  1  is  given  by 


^*1  <*,■  V  =d/imGim  %  <rr  <h>  e‘ 


m 


(11) 


and  similarly  the  edge  diffracted  field  at  edge  No.  2  is  given  by 

ik 


<lTe2  (r^,  <j>2 


r? *  ~ 


m 


Glm  Ue2  (r2  ’  ^2)  S 


m  m 


(.12) 


Note  that  is  not  the  complete  formula  for  the  Gaussian  curvature  for 
the  primary  field.  The  Gaussian  curvature  is  given  by 

GAl  =  |jk  Rj)  (  k  r2  sin  .  See  Appendix  II  for  GAe,  the 

Gaussian  curvature  for  the  edge  field,  which  is  obtained  similarly. 
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w 


The  edge  diffracted  field  at  P  is  (ri  and  r2  2  Xi  ) 
°  A  i  max 


Afe  s  ATej  (rj.  <t>i )  ±  (Ue2  (r2> $2) • 


(13) 


Whether  the  plus  or  minus  sign  in  (13)  is  used  depends  on  the 
vector  symmetry  of  the  problem.  In  the  acoustic  case  we 
always  use  the  plus  sign..  In  the  cylindrical  case  we  generally 
would  use  the  plus  sign  also.  In  the  axially  symmetric  case, 
when  the  incident  vector  field  on  the  edges  is  polarized  identically 
with  respect  to  the  edge  for  all  rotation  angles  (vector  symmetry) 
then  we  must  use  the  minus  sign.  When  the  incident  field  is  lin¬ 
early  polarized  in  the  axial  case  (E-field  in  same  absolute  dir¬ 
ection  everywhere)  we  must  perform  an  averaging  of  the  fields 

3 

around  the  edges  such  as  was  performed  by  Rusch  .  How  this 
may  be  done  will  be  discussed  subsequently.  For  the  axially 
symmetric  case  near  the  axis  we  cannot  use  (13)  accurately, 
since  the  axis  is  a  caustic  of  the  system.  A  caustic  approxima¬ 
tion  will  be  presented  shortly  and  is  discussed  in  Appendix  II. 


Our  next  steps  should  be  to  first  let  P  become  far  from  the 
edges,  i.e.  ,  let  (krj)‘^>>l  and  (kr^)’'  ^*>1  (this  is  actually 
implied  already  when  we  assume  r,  and  r^  2  Xi  ).  This 
is  consistent  with  geometrical  optics  approximations.  To  deter¬ 
mine  Je  <<&>•  C«o  (<&■  and  de  (  )  we  must  first  deter¬ 

mine  the  phase  contours  of  the  far  field.  This  will  complete 
the  determination  of  the  edge  field.  Before  proceeding  to  det¬ 
ermine  these  quantities  we  should,  however,  clarify  three 
important  details  about  the  expressions  for  Ue^  and  Ue^.  The 
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first  concerns  the  precise  form  of  Ue  when  the  indicent. 
plane  wave  polarization  is  neither  normal  nor  parallel  to  the 
plane  of  incidence.  The  second  question  concerns  the  precise 
value  of  <j>.  near  G^,  the  shadow  boundary,  at  which  point 
equation  (6)  must  be  replaced  by  (8).  The  third  question  con¬ 
cerns  the  manner  in  which  the  caustic  axis  for  the  axial  system 
will  be  accounted  for. 


B.  Polarization  at  an  Angle  0  with  Respect  to  the  Edge 

In  order  to  resolve  the  question  concerning  the  incident  polar¬ 
ization,  let  us  define 


A  =  -  IL+iL  GA  e 


i  k  r 


sec  (  -k-yr-) 


and  (14) 


B  =  - 


(1  +  i)  /oikr 

4VV 


oAp  e1 


i  _ 


.  cj>  +  c ^  \ 

sec  (  — - * - ) 


(15/ 


when  cj)  is  not  near  Gr,  and 


B  .  -  lll±  ga  eikr 

4  vV 
3tt 


\[2?ga^1  e"ikr(i+cos[<i>+c7C]j 

(15a) 


e  4 


.  ir 
-1  — 
4 


n  b  GAe  c°s 


exp  (i  y  qf)  d  %  -1 


Co 


when  *  is  near  G^. 
T  r 
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The  quantity  Ga6  is  the  appropriate  "curvature"  for  the 
system  as  explained  in  Appendix  II,  where 


G  a6  -  1  /  (k  r  s  in 


1/2 


f""t  I  )  (^n  the  far  zone  away 


( 1 4a) 


from 


rc 


=  0  > 


The  angle  /  is  the  angle  from  the  horizontal  defined  in 

equations  (26)  and  (27).  See  also  Figure  1  of  this  chapter. 

The  validity  of  (15  a)  when  =  l/(kr  sin^^iP  I  ) 

^e  I  2  1 

g 

is  justified  by  the  asymptotic  form  for  (15a)  (see  Sommerfeld  ) 
in  a  three  dimensional  system.  Equations  (14),  (15),  and  (15a) 
are  identical  with  (6)  and  (8)  when  G.  is  accounted  for. 


When  the  incident  electric  vector  is  polarized  parallel  to  the  edge 
we  have 

Ue  =  UeH  =  A  -  B  (16) 

and  when  the  electric  vector  is  polarized  normal  to  the  edge  we 
have 

ue  *  Uei  =  A  +  B  (17) 

Now  let  the  incident  electric  vector  be  polarized  at  an  angle 
with  respect  to  the  edge.  Then  we  find 


Ue  j  j  =  (cos^)(A+B)  and 

(18) 

UeJ_  =  (sindHA  +  B). 

(19) 

* 


See  Appendix  II. 
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Hence  we  see  that  the  edge  diffracted  field  polarization  does  not 
remain  fixed  for  all  cj>.  In  this  event  we  will  consider  the  field 
component  that  is  normal  to  the  (^-direction  along  a  ray  as  a 
cross -polarized  component  which  mast  be  accounted  for.  The 
■desired  component  in  the  $ -direction  is  found  from 


U0  cos 
ell 


=  u« 


1 1  # 


cos 


¥  (A  -  B) 


and  Ue  sin  ^  =  Ue  ^  ^ 


so  that 


=  sin 


2  T  (  A  +  B) 


U, 


Uell\  + 


UriY  =  A  ’  B  (cos  2 


(20) 


When  the  polarization  is  either  normal  or  parallel  to  the  edge, 
there  is  no  cross -polar ized  edge  field.  It  is  clear  that  in  case 
the  preferred  polarization  is,  for  example,  right  circularly  pol¬ 
arized,  then  a  "cross-polarized"  left  circular  component  will 
exist  in  the  edge  diffracted  fiends.  (Note:  If  a  right  circular 
wave  is  incident  on  the  edge,  then  the  principle  polarization  ref¬ 
lected  is  left  circular  as  it  is  for  reflection  from  any  conducting 
surface.  This  fact  was  taken  into  account  in  the  previous  statement.) 

Equation  (20)  is  suggestive  of  an  "averaging"  method  which  can 
be  used  when  the  incident  polarization  is  not  rotationally  sym^ 
metric  or,  in  other  words,  those  cases  where  ^  =  ^(9)  and  0 
is  the  angle  of  rotation  about  the  axis  of  the  axially  symmetric 
system.  For  example,  if  the  incident  wave  on  the  subreflector 


See  equation  (20a)  which  follows  (24), which  follows  . 


is  linearly  polarized,  that  is  u-  0,  then  we  may  use  as  an 
average  edge  field 
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Ue  =  A, 

for  on  the  average  we  do  have  ch  tt/4  and  y-  3 tt/ 4  for 
the  top  and  bottom  halves  of  the  reflector.  When  using  Ue  -  A 
we  would  also  use 

ATe  =  /Vej  +  /lTe2 

from  equation  (13),  and  hence  would  not  get  a  zero  on  the  axis 
or  caustic  from  this  equation. 

C.  Field  Near  the  Shadow  Boundary  <j>  -  Gr 

We  will  determine  a  value  of  cj>  ^  Gr  at  which  we  will  exchange 
(15)  for  (15a)  in  evaluating  the  field.  This  value  is  somewhat 
arbitrary  and  is  chosen  for  analytical  convenience  and  for  reason¬ 
able  approximation.  Before  considering  this  question  further, 
it  will  be  necessary  to  decide  upon  a  value  of  r  at  which  to  eval¬ 
uate  B  when  cf>  =  tt  . 


At  =  tt  „  the  value  of  B  at  the  aperture  is  given  for 

r  =  =  f2m*  Since  the  main  reflector  is  in  the  far  zone 

V  max  { 

of  the  subreflector,  the  r  dependence  of  equation  ( 1  5a)  is  given 
essentially  by  the  factor 


(1  +  i) 


ga  G 


i  k  r 


4  /7" 

when  cj>  &L)  -  Gr  (see  the  asymptotic  expansion  of 

Sommert’eld^) .  Hence  for  our  purposes  we  may  use 
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B  =  -  GAe  S 


i  k  r 


4  >nr 


-i 


'B?  °Ae  I  fij  <£ 


k  Pirn  [i  +  cos(9  +  ^)]3ir 


-1 


.^re' 


IT 

1  j 


C=  °Ae  (p2m)cos(^-) 


-  1 


o  exp  (i  7  'T ')  d  ^ 


e  4 

(15  b.) 


where 


-1 


,  Z  ^  2  1 

p2m>  =  \ 

k  [2m  Sm  1  2m 

and 


(15c) 


r 

0 

1  2m 

2m 

(See  Appendix  II). 


(15d) 


Our  problem  now  is  to  connect  equations  (15)  and  (15b).  The  best 
way  to  connect  these  equations  would  probably  be  to  find  two  values 
of  4>  for  which  they  are  equal  in  amplitude  and  in  phase.  In  this 
way  we  could  maintain  the  continuity  of  both  the  amplitude  and 
phase  of  the  edge  field.  Connecting  the  amplitude  as  indicated 
in  the  diagram  iollowing  (e)  would  in  principle  be  simple,  although 
connecting  the  phase  may  lead  to  some  difficulties.  Whereas  the 
range  4>  ^  4>  is  small  for  the  amplitude  connection,  it  may 

be  large  for  the  phase  connection.  For  the  sake  of  simplicity  in 
analysis  and  because  the  range  pr  4>  —  Gr  will  be  small,  we 
will  use  instead  a  somewhat  simpler  connection.  For  the  phase 
of  B  in  the  region  <|>  ^.Cp  we  will  use  the  phase  of  equation 
(15).  This  phase  is  constant  over  the  entire  range  0^<j>*=Gr. 

This  is  justified  because  the  region  <|>  <^<t>s&Gr  is  very  small;. 
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usually  less  than  1°  or  2°.  Hence 


phase  of  B 


(_il  +  kr)  for  ("0  —  ^  Gr ) . 


(21) 


We  will  connect  the  amplitude  by  assuming  that  the  amplitude 
is  constant  in  the  region  4^  —  4*  ^  Gr  and  is  given  by  equation 
(15b)  evaluated  at  <j>  =  G^.  The  value  of  4>r  is  determined 
by  finding  the  intersection  as  shown  in  the  diagram  on  the  follow¬ 
ing  page  (f)4  If  we  evaluate  the  value  of  (15b)  at  cj)  ~  Gr  (the 
integral  term  vanishes  at  <j>  ~  when  r  is  finite)  and  set 
the  amplitude  of  this  value  equal  to  the  amplitude  of  equation 
( 1 5)  we  find 


cj>r  ~  +  2  arc  sec 


.  1/^  , 

?m  sm  N> 


2  m 


\TZtt1  k 

and  the  magnitude  of  B  in  the  range  ^Gr  given  by 

( 1+i )  „  _ikr 


(22) 


B 


4/tP 


GAe-e 


{fZi?  GAe  ( 


(23) 

for  (4>r^:(t)^Gr) . 


In  summary  then,  we  may  state  that  Ue  is  given  by  (20)  where 
A  is  given  by  (14)  and  B  is  given  by 


(A+i)  G  Ae<2 


i  k  r 
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VL 


4/^r1 
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4vTrP 
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sec  (1 ±sk) 
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for  (  0*t(|><r4>  ) 


(24) 


^  Ae  ^ 
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'Jt*1  GAe  (Am)  i  for  (ct»rz*4>ZGr) 


where  is  given  by  (22)* 
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D. 


Before  proceeding  to  determine  Je  (  J^),  de  (  (fe),  and  Ceo{(j?2), 
two  additional  points  concerning  the  range  of  <j>  for  equation 
(20)  (and  the  above  equations  for  A  and  B)  should  be  mentioned. 
As  was  indicated  in  equation  (7),  we  require 

0^-^2tt  (7) 

Because  of  the  fact  that  <j>  is  divided  by  two  in  the  equations,  the 
solutions  are  not  periodic  in  2it ,  but  rather  over  4tt.  Hence 
we  shall  define 


Ue  =  0  unless  (  0^4>^2tt).  (20a) 

Taking  account  of  (20a)  may  become  important,  for  example,  if 
the  subreflector  is  very  curved  at  the  edge  and  the  ray  from 
edge  No.  2  in,  the  direction  of  Q^max  (toward  X2m  re^ector 
No.  2)  is  blocked  or  shadowed.  In  this  event  the  value  of  <j>  for 
this  ray  from  edge  No.  2  will  be  negative  and  (20a)  would  indicate 
that  Ue  =  0  as  desired. 


Value  of  G 


Ae 


Equation  (24)  specifies  two  separate  ranges  of  values  over  which 

two  different  equations  for  B  are  valid.  There  are  two  possible 

expressions  for  G^e  depending  upon  whether  -0  or  not 

near  the  caustic  at  P?  =  0)  and  whether  we  are 

considering  the  top  edge  (LL)  or  bottom  edge  (tLJ. 

m  e  2 


We  have  in  general  that  G^e  ^  ^ 


(kr)/(kr  sin  |  ir  k  Xjm) 


where  the  +  sign  is  for  U  and  the  -  sign  is  for  U  . 

.1  *'  e  2 
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For 


sm 


r. 


>  > 


Xlm/r 


-1  1/2 

we  have  G^e  —  kr  sin 


r2 


as  the  far  zone  approximation  off  the  caustic.  Near  the  caustic 
for  any  finite  r  we  should  use 


GAe  =  kr 


sin  P,  ±  (Xlm/ 


In  the  vicinity  of  «  0  we  have  as  an  approximation  for 
r  =  rc  ^(ckC+^3)  +^3*  An  estimate  for  rc  may  be  taken  as 


r  =  A  . 
c  Zm 


(24a) 


More  suitable  estimates  may  be  made  for  specific  problems. 


When  =  0,  GAe  has  th'e  same  value  for  the  top  and  bottom 
edges  (U  and  Ue^).  Away  from  P^  -  0  when 

J  s  in  r2 1  |.Xim/  rc|  ,  the  value  of  Ga^  is  again 

independent  of  the  edge  designation.  As  an  approximation  in  the 
small  angular  sector  near  the  caustic,  we  will  take 


for  both  edges  and  thus  simplify  our  computations  considerably. 
For  the  case  of  vector  symmetry,  this  will  give  f\Jfe  ~  0  (see 

equation  (13)  ^  for  =  0  which  is  correct  so  the  error  will  be 
over  a  small  angular  region  comprising  a  very  small  solid  angle. 
This  is  explained  further  in  Appendix  II.  We  choose  the  +  sign 
in  equation  (24b)  to  avoid  the  infinity  which  would  occur  with  the 
-  sign  when  the  field  of  edge  No.  2  is  evaluated  near  =  0, 
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E.  The  Total  Edge  Field 

In  order  to  compute  the  actual  edge  diffracted  field  in  the  far 
zone  of  the  two  edges  combined,  we  let  and  hence 

rays  r^  and  to  point  P  from  each  edge  become  parallel. 

We  designate  their  common  angle  with  the  L axis  as  Hj. 

In  accordance  with  Figure  1  and  with  reference  to  the  distance  R 
from  the  origin,  we  have 


rl  =  R  -  xlm  sin  l~2 

and 

r2  =  R  +  Xlm  sin  fl. 


(25) 

(25a) 


Equations  (25)  and  (25a)  are  used  to  replace  r^  and  r^  only 

i  k.  r 

in  the  phase  factor  indicated  in  the  expressions  for 

Uej  and  Ue2  which  are  obtained  from  (14)  and  (24).  The  r  term 
in  the  equations  for  Gy^e,  equations  (24a)  and  (24b),  are  replaced 
by 


R 


(25b) 


in  accordance  with  the  customary  far  field  assumptions. 


From  Figure  1  we  also  readily  find 

4>1  =  -f  +  r2  t(-e^p6;m) 

<t>2  -  —  -  f2  +  (itm  -  e2m, 


TT 

T 


(  9lm  +  e2m  j 


and 


(26) 


(27) 


C^'  = 


2 


(28) 
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With  (25)  through  (28)  substituted  into  the  expressions  for 
A  and  B  in  (14)  and  (24)  we  obtain 


Ai  = 


i  k  r 


(i+  DGAee 

4  \|*tP 


■  i  k  X  i  Yn  ®  in.  R 


e. 


sec(-S..t  9lm.) 


A2  = 


B 1  = 


i  k  r 


(l  +  i)GAee 

L  4^ 


+ikXlm  sin  Q 


i  k  r 


( 1  +  i)GAe^ 


L  4vTfP 


e 


ilkXlmsin  6 


sec  (  A  ~  9!m  ) 


-  esc  ( — A— 1  6 2m  ) 


(30) 


(31) 


for  (  0  2<|>i  <4>r) 


B2  = 


V. 


and 


ik  r 

sin  r2 

— 

— T 

(i+i)GA,e 

+ikX. 
p  1  m 

/ 

.?  +  ®2rp  ) 

4^~iT7 

CSC  (  — 

(32) 


for  (  0<cct>2<:<t>r) 


B,  = 


(l  +  i)GAee 

4  n 


i  k  r  |  . ,  v.  P 

-lkXim  Sln  '2 

le 


-i 


^  G  fp2m) 


(31a) 


for  (  ?G  ) 


V®  2 


i  k  r 

( 1  +  i)G Ae  e 

+  ikxlm  sin  (1 

e 

^  G'1  (  pZm) 

4  T7 

{ 

(32a) 


for  (  4>r  ^Gr) 


From  equations  (11),  (12)  and  (20)  we  get 

Ue  ^  =  A^  “  Bj  cos  2  $ \  and  Ug^  =  A^  ~  ^2  2  or 


/u 


ei 


(j+lAn  Gim  GAe^ 


m 


4fP 


-ikXimsin  ?z 


'H±i 


lxa) 


(33) 


sec  (■ - ^ 

+  esc  (  —  2m )  cos  2 
for  (Oai(qZ(jir) 
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rOi 


n  iMK  +  tl 

n 

-  ( 1  +  i)<^l  1  m  Gl  xn  GAe 

+  ikXim  sin  \z  : 

b 

4^ 

sec( 


^2  ~  9-j 


mi 


(34) 


a r 


e  i 


-( l  +  i)i?im  Glm 


ik(R+^mn 


4  fP 


e 


m  sinQ 


-  esc  2mj  cos  2<Y 


for  (0 ^<\>z  ■£.$■?) 
sec  9lpri)  -  fin  (33a) 

G"‘<P2m»  cos  2<^ 

for  (4>r^4>iS:Gr) 


ATe2  = 


-(1+i  >3  lm  Gl 

m  pAe 


ik(R  +Tm 


4J1P 


+ikxlm  sin  ^2 


sec  (  ^  ~,6  lm) 


(34a) 


-  \f2rr  G-  1  (pZvn)  cos  2  / 
for  (4>r^.  ct)z^Gr). 


For  the  sake  of  convenience  in  future  manipulations,  we  define 
from  equations  (33)  and  (34), 


z 

z 


1 


2 


sec( 


J?  1 9.1™.)  +  esc  ( JhlQ Zm.  ■)  cos  2</ 


sec(—  -  ilSL)  -  esc  (-^  '-.r^mx  cos  2  0 
_  2  2  1 


-Bt^Zxnx  cos  z  (¥ 


sec(—  +2  9^m)  GAe  (  Am)  cos  2/ 


for  (0  z<(>r)  (35) 


for  (0  <<)? 2  s4>r)  (36) 


for  (4>r  <  4>  ].  ^cGr)  (35a) 


sec( 


(1  -lf2?GA1e  (  ^2m)  cos  2(P 


z,. 


for  (4>r  ^<j)2  2"Gr)..  (36a) 


Note  that  Zj.j  is  the  field  contribution  from  edge  No.  1  aird  that 
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f 

— *2  *s  t^ie  field  contribution  from  edge  No.  2.  A  common  factor 
from  each  of  these  contributions  is  omitted  as  may  be  seen  in 
equation  (37)  below. 


From  equation  (13)  we  get  for  the  edge  diffracted  field  subject 
to  the  conditions  of  symmetry  specified  earlier,  the  result 


where  the  appropriate  values  of  Z_q  and  Z_i2  are  to  be  used 
from  equations  (35),  (35a),  (36),  and  (36a)  depending  on  the 
values  of  <j>i  and  <J>2- 


The  edge  diffracted  field  given  by  (37)  is  complex  and  no  fixed 
phase  center  exists  for  this  field.  Hence  the  phase  contours  are 
not  circles.  Therefore  P2  £  From  (37)  we  can  derive 

n2  (  ^e  an<^  Ce Q ( (^2 ) •  The  value  of  ^  e  (  or  Je  (  ) 

can  be  found  directly  from  (37).  When  P>  (  <j^)  is  found  sub¬ 
sequently,  then  we  can  obtain  (  e((j^)  or  Je  (<j^).  As  will  be  made 
more  evident  shortly,  we  cannot  obtain  J e  (d^)  from  Je  (f^)  simply 
by  substituting  /^>  (<J^)  for  P^ •  The  transformation  is  more 
complex. 


72 


f.  j;  ( p2) 

i  2 

We  will  obtain  Je  (  f?)  =  jj  e  (  fl)  directly  from  (37).  There 
is  some  arbitrariness  in  this  determination  depending  on  our 
definition  of  1^(0-^)  =  ^  (0^)  from  the  given  primary  field.  We 
have  assumed,  however,  that  the  primary  field  is  obtained  by: 


T  1  pik(rl  +  GPo)  a 

E1(61)  =  Gi  exP  <ik  Lri  +  CPo  (0lU  )ilj  <9l>  #1  <38> 

1  jkRi  kr^ 

With  this  definition  of  ^  (0j)  we  factor  1/kR  and  obtain  from 
equation  (37), 


Je  (  r2)=  <^)q?«imGL  [*2l  ®«2  (k  X^sin^+I-Zj  ±\) 


(39) 


sin^  (k X  sin 


in  r2)J 


where 


1 /  (sin 


H>  +  X,  /r,). 


2 | 1  ~lm 

As  we  expect,  this  power  flow  is  of  a  higher  order  in  1 /k  than 
the  primary  field. 


(39a) 


G.  Reference  Wavefront  of  the  Edge  Field  -  R  (  f^)_ 

We  desire  to  find  the  reference  phase  contour  curve  given  by 
R(  Q).  As  pointed  out  earlier,  the  curves  of  constant  phase  are 
not  circles  for  the  edge  field.  We  will  choose  there,  a  particu¬ 
lar  curve  given  by  the  phase  being  set  equal  to  zero  for  equation 
(37).  Although  any  phase  value  is  adequate,  the  value  zero  is 
the  reference  contour  for  which  the  distance  from  the  curve  to 
the  axis,  along  a  normal  to  the  curve,  is  equal  to  "CeQ(  Jf^)> 
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the  reference  phase. 


To  simplify  the  notation,  the  following  substitutions  are  made; 


?  =  kXlm  SinPZ 
y  =  k  (R  +  ^)  +  rr/4 

Yl  =  Z1  ±Z2’  YZ=  *  Z1  ±Z2 ‘ 
Equation  (37)  may  then  be  written  as 


(40) 


(40a) 


(40  b) 


-  Sim  G 1  m  OAe 
4/71 


[y^  cosJ^cos^  -  sinjfsin^ 


r 

-|  Y  j  sin  1  c os  ^ 

Up. 

+  Y^  sinj  cos  74 

J  ^ 

Note  from  the  above  equation  that  the  phase  distribution  of  (XT' 
is  determined  in  large  measure  by  the  relative  magnitude  of 
Yj  and  Y^J  the  sum  and  the  difference  of  the  edge  field  con¬ 
tributions  from  each  edge. 


The  phase  of  P[^  is  equal  to  zero'  if 


Yj  sin^fcos^+Y^  s  in  ^  c  o  s  =  0 


Hence  we  obtain 


y-  - 


arc  tan  ( 


Y  z 


Y 1 


tan  £  ) 


or 


* 


Actually  zero  or  180°.  The  following;  additional  condition  is  required 
for  zero  degrees:  Fy^  cosj0  cos  f  -  Y2  sinTsin£]>0- 
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R  =  R(  ^2.)  =  arctan  (  tan  £>, 


(41) 


(A  convenient  branch  for  should  be  chosen,  although  this 

can  be  seen  to  be  arbitrary.  Once  is  chosen  it  is  generally 
preferable  to  maintain  R ( f^)  continuous  by  the  appropriate 
branch  choice  in  (41).  It  is  important,  however,  that  the 
branch  choice  for  j  be  chosen  so  that  R  is  far  from  the 
caustic  of  the  field;  near  would  be  a  good  choice  for 

example). 


H .  The  Transformation  Function  P?  (  ^ 

The  transformation  function  P>  (  <^>)  can  be  obtained  from 
the  phase  contour  R(  f^)  given  by  (41).  We  need  in 

order  to  obtain  de(^,),  and  Ce„(d/,). 


evu2/ 


From  the  diagram  following  (g),  we  note  that 


P2  =  -J-  -  A2  and  tan  A2  =  ctn 

7=  & 


(42) 


;2-  (42  a) 

If  y-j  (x^)  is  the  rectangular  equation  for  the  curve  R(  f^),  then 

dy-j  =  dR  sin  +  R  cos  A^  d  A 

dx2  =  d  R  cos  Ag  -  R  sin  A^  d  A z>  and 


(42b) 


dy 


dx 


-t-  =  tan  (it  -  ^  )  =  -tan 


(42c) 


Equations  (42)  through  (42c)  yield  for  (  P2) 


which  gives,  numerically,  the  result 


r2  =  r2  ( <& 


(43a) 


In  order  to  evaluate  (43)  we  find 


dR 


-1 


dP  "  kd+  tV/YlT*“5> 


-  Y  j  Y2  tan  5 


(44) 


Y; 


sec 


+  h  tan^  +  Xl 

Y 1  Y  i 


where  the  prime  denotes  differentiation  with  respect  to 


The  Reference  Phase  Functions  de  and  Cen 

The  reference  quantities  de  (  anc*  CeQ  (  can  also  be 

obtained  from  a  knowledge  of  R  (  f^),  equation  (41).  If  we  refer 
to  the  figure  in  part  H  above  and  apply  the  law  of  sines  to  the 
diagram  we  obtain  first 


_ R 

sin  (tt  -  (¥2) 


Jo^  +  ^S)  -  dp  _ 

sin  tt-(tt-J^)-  1 


Note  that  ^3)  is  a  given  constant  for  the  reflector  system. 

From  the  above  equation  we  obtain 


=  (<^+Z3) 


S  in  ( (pz  +  r2) 


sin 


R 


(45) 


where  R  is  given  by  (41). 


If  we  apply  the  law  of  sines  a  second  time  to  the  same  figure 
we  obtain 


sin  (ir  -  (K) 


-C 


eo 


sin  r , 


which  gives 


sin(  I2) 

_  sin  ( 

R 


(46) 


In  both  equations  (45)  and  (46)  for  de  and  CeQ  we  note  that 
in  order  to  obtain  de  (  (j^)  and  Ce.0  (  (j^) .  we  must  first  solve 
equation  (43)  numerically  for  (  [p. 


(  £.2) 

Having  found  (  P^),  equation  (39).  and  (jf(  p^),  equation  (43), 
we  are  prepared  to  determine  Je  (  <^K) . 

First  we  will  explain  why  Je  (  J^,)  ~  Jg  (  The  9uan,;i~ 

ties  J  and  J  represent  power  flows  per  unit  solid  angle.  At 
a  point  on  the  wavefront  given  by  the  values  anc*  R  ( 

Jg  (evaluated  at  this  value  of  P^)  gives  the  value  of  the  power 
flow  through  a  solid  angle 

s  in  d  d  0 

Here  we  have  d9  as  the  angle  of  rotation  about  the  axis  of  the 
system.  At  the  same  point,  the  same  power  flows  through  a 
different  solid  angle  given  by 
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Je  designates  the  power  flow  per  unit  solid  angle  referring  to 
this  latter  angle.  This  is  evident  from  the  diagram  (h)  on  the 
following  page , 

Since  the  total  power  flow  through  both  solid  angles  is  the  same 
we  have 


Jg  (  P2)  sin  r2  d  r2  *  Je  (  <£>  sin  d  ^ 


With  considerable  manipulation, 


can  be  evaluated  directly  from  (43). 


(47) 


(48) 


With  d0(  (jf2),  Ce  ('^),  and  Jg(  determined  from  equations 
(45),  (46),  and  (47)  together  with  the  auxiliary  equation  (43)  for 
(b  (  P2)  we  have  obtained  the  complete  edge  field  in  a  form  which 
will  allow  us  to  compute  the  edge  field  distribution  in  the  aperture 
of  the  main  reflector  after  reflection  from  this  reflector. 
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IV.  Computed  Edge  Diffracted  Field  Results 

In  this  section  the  results  of  numerical  computations  for  the  edge 
diffracted  field  will  be  presented  and  briefly  discussed.  The  compu¬ 
tations  presented  here  will  not  include  the  edge  field  scattered  from 
the  main  reflector  into  the  aperture.  The  presentation  of  those  results 
is  deferred  until  the  next  section. 

The  results  presented  here  include  the  diffracted  fields  for  subreflec¬ 
tors  of  various  sizes  in  wavelengths  and  for  various  reflector  shapes. 
They  are  plotted  and  illustrated  as  a  function  of  from  0°  to  180°. 
Between  about  170°  to  180°  the  results  are  in  error.  The  oscillations 
are  180°  out  of  phase  with  the  true  oscillations  of  the  tield.  This  is 
a  computational  error  not  due  to  the  edge  diffraction  theory.  The  aver¬ 
age  field  in  this  area  (averaged  over  the  oscillations)  is,  however,  cor¬ 
rect,  so  this  portion  of  the  plots  was  included  in  any  case.  Hence  we 
will  observe  the  back  radiated  field  away  from  the  main  reflector 
in  addition  to  the  field  incident  on  the  main  reflector. 

It  should  be  stressed  that  the  geometrical  optics  field  has  been  sub¬ 
tracted  from  the  total  field  scattered  from  the  reflector.  Only  the 
diffracted  or  second  order  field  is  considered. 

A  large  number  of  input  parameters  or  boundary  conditions  are  nec¬ 
essary  for  each  computation  of  a  curve..  Most  of  these  parameters 
remain  fixed  for  most  of  the  curves  computed.  It  will  be  far  simpler 
and  clutter  the  graphs  much  less  if  a  set  of  standard  values  for  these 
parameters  is  adopted  and  only  deviations  from  this  set  of  values  is 
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mentioned  on  the  graph  and  in  the  ensuing  discussion.  A  similar  pro 
cedure  is  followed  in  the  next  section  and  later  chapters,  so  we  will 
refer  to  the  following  parameters  as  Standard  Set  No.  1,  and  present 
these  values  in  the  following  table: 


STANDARD  SET  NO.  1 


SYSTEM  =  Axially  symmetric 


lm 


Xl 


m 


30° 


+  1.00 


^  =  <^  +  (3  =  fP 


*2 


m 


+  10.  0 


=  constant  =  0 

Ix  (Qj)  =  cos16  d1  so  ^  (0lm)  *=  1  / 1 0 

SIGN EG  =  +1  for  respectively  cancellation  or  addition  of 

the  edge  field  from  each  edge  at  the  center  (  =  0). 

It  will  be  chosen  as  +1  for  addition  as  the  standard 


j' 


edge. 

polarization  angle  with  respect  to  the  ege.  It  is  chosen 
as  0°  for  the  standard  value. 


A,  Diffraction  of  Flat  Discs  and  Strips 

The  edge  diffracted  field  of  a  plane  wave  incident  on  a  flat  con¬ 
ducting  strip  or  disk  is  of  no  direct  interest  to  the  designer  of 
dual  reflector  systems,  but  does  serve  the  purpose  here  of  check¬ 
ing  our  computational  procedure  against  previously  computed 


7 

work  and  also  checking  the  theoretical  accuracy  of  this  method 
of  computing  the  diffracted  field. 

7 

Keller  computed  the  diffracted  field  of  a  plane  wave  incident 
on  a  conducting  screen  with  a  strip  aperture  and  circular  aper¬ 
ture.  The  results  depicted  in  Figures  3,  4  and  5  are  for  the 

diffracted  field  of  a  plane  wave  on  a  conducting  strip  and  circu- 

15 

lar  disk,  the  complementary  screens  ..  By  Babinet's  principle 
the  results  are  equivalent  with  the  principle  exception  being  that 
the  backscattered  field  in  one  system  is.  equal  to  the  forward  scat¬ 
tered  field  in  the  other.  (The  polarizations  must  also  be  altered, 
but  this  makes  little  difference  for  these  cases).,  Hence  Keller^ 
presents  his  results  for  the  forward  scattered  field  only  and 
these  compare  very  well  with  our  backscattered  fields  (  = 

0°  to  90°)  depicted  in  Figures  3,  4  and  5.  In  addition  the  results 
depicted  in  these  figures  extend  the  published  results  of  Keller 
to  the  region  90°  ^  180°. 

For  the  case  of  the  metalic  strip  with  k  =  8,  Figure  3,  the 

depicted  results  compare  very  accurately  with  the  theoretically 

17 

exact  results  published  by  Karp  and  Russek  .  This  is  signifi¬ 
cant  since  the  Kirchhoff  theory^  predicts  substantially  different 
results.  The  Kirchhoff  theory  predicts  zeros  in  the  field  whereas 
the  Keller  theory  predicts  non-zero  minima  in  approximately  the 
same  angular  regions.  The  Kirchhoff  theory  predicts  a  fixed 
phase  center  at  the  center  of  the  aperture,  whereas  the  edge  dif¬ 
fracted  field  as  mentioned  earlier  does  not  have  a  fixed  center  of 
curvature  for  the  phase  front.  This  example  is  important  for 


83 


oar  purposes  because  it  demonstrates  the  breakdown  of  the 
Kirchhoff  theory  for  small  apertures  or  reflectors  (  k  =  8 
implies  a  diameter  or  slit  width  of  approximately  i.6  wave¬ 
lengths)  and  at  the  same  time  illustrates  the  success  of  the 
edge  diffraction  theory.  Since  we  wish  to  examine  the  diffrac¬ 
tion  effects  for  small  reflectors  of  this  order  (two  or  more  wave¬ 
lengths)  in  diameter  the  distinction  between  the  theories  is 
critic  al. 

For  larger  diameter  reflectors  and  for  curved  reflectors  the 
two  theories,  the  Kirchhoff  and  edge  diffraction  theories,  agree 
well  as  will  be  demonstrated  in  the  next  part  (B),  However,  it 
must  be  remembered  that  the  edge  diffraction  theory  begins  to 
fail  badly  when  the  amplitude  of  the  incident  wave  varies  very 
rapidly  near  the  edge.  For  example,  if  the  incident  wave  has 
a  zero  at  the  edge,  then  the  edge  diffraction  theory  gives  a  zero 
diffracted  field  which  is,  of  course,  wrong.  The  Kirchhoff  theory 
gives  essentially  the  same  answer  for  the  diffracted  field  whether 
the  field  varies  rapidly  or  slowly  near  the  edge.  That  is,  its 
solution  is  not  ciitically  dependent  on  the  edge  condition- 

With  regard  to  Figure  5,  those  readers  who  compare  that  figure 
with  the  same  case  illustrated  in  Reference  7  will  note  some 
small  differences.  These  may  be  attributed  to  the  difficulty  in 
setting  up  the  computer  program  for  a  plane  wave  incident  on  a 
flat  reflector  in  the  axially  symmetric  case.  The  program  was 
designed  for  the  input  parameters  of  more  practical  dual  reflec¬ 
tor  designs  such  as  those  parameters  in  Standard  Set  No,  1  and 
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large  changes  in  the  program  were  not  warranted  to  merely 
correct  small  deficiencies  in  one  curve. 

B.  Edge  Diffracted  Fields  from  a  Hyperboloid 

(for  comparison  with  the  Kirchhoff  theory  and  experiment) 

The  edge  diffracted  field  of  a  conducting  hyperboloid  of  revolution 
is  presented  in  Figures  6,  7,  8  and  9  for  hyperboloids  of  diameter 
8  A  and  16  A  •  The  8  A  results  in  Figure  6  are  for  3  0 
(H  -plane  cut),  and  the  16  A  results  in  Figures  7,  8  and  9  are 
for  tt/2  (E-plane  cut).  The  three  Figures  7,  8  and  9  pre¬ 

sent  varying  degrees  of  detail  for  the  16  A  results. 

These  results  are  intended  for  comparison  with  the  results  that 

Rusch  obtained  by  Kirchhoff  integration  and  also  with  the  exper- 

1 8 

imental  results  that  Rusch  has  obtained  .  (Rusch  presented 
similar  results  in  Reference  3).  The  agreement  between  the 
results  of  the  two  theories  is  generally  good,  both  theories  pre¬ 
dicting  the  same  major  characteristics  of  ripple  and  high  back 
lobe  (  Q  >-120°)  as  do  the  experimental  results  also.  For  those 
readers  making  an  actual  comparison  between  the  results,  it 
should  be  noted  that  Rusch  presents  his  results  with  the  geomet¬ 
rical  optics  field  included  as  is  characteristic  for  results  obtained 
by  Kirchhoff  integration. 

It  might  be  added  that  Figure  6  through  9  (and  later  results)  are 
the  first  results  obtained  by  edge  diffraction  theory  for  curved 
reflectors  known  to  the  author  and*  hence,  these  results  do  extend 
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the  known  usefulness  of  the  theory. 

One  important  result  seen  in. the  figures  is  that  the  forward  dif- 
fracted  field  (  Q  ^90°)  is  as  large  as  or  greater  than  the  field 
diffracted  toward  the  main  reflector  (  ^  90°).  The  forward 

diffracted  field  is,  however,  spread  out  over  the  angular  regions 
depicted  in  the  graphs,  whereas  dual  reflector  systems  designed 
for  high  gain  (easily  better  than  99%  of  all  such  antennas)  will 
focus  the  energy  diffracted  into  the  main  reflector  and  hence 
increase  its  significance  in  the  main  lobe  area  of  the  tertiary 
pattern.  The  forward  scattered  field  contributes  to  the  overall 
system  sidelobes. 

C.  Edge  Diffracted  Fields  as  a  Function  of  Frequency 

Consideration  of  the  edge  diffracted  field  as  a  function  of  frequency 
or  k  is  of  great  importance  to  this  general  study.  One  of  the 
objectives  as  explained  earlier  was  to  examine  the  possibility 
of  using  a  smaller  than  usual  subreflector  for  the  dual  reflector 
system  in  order  to  reduce  blockage  by  the  subreflector.  In  addi¬ 
tion,  a  general  understanding  of  diffraction  effects  for  any  size 
reflector  is  of  great  importance. 

With  the  above  considerations  in  mind  a  study  of  the  edga.-ii.eld 
as  a  function  of  k  was  made  and  is  presented  in  Figures  10 
through  24,  In  addition  to  plotting  Jg  (  p^)  for  from  0°  to 

180p,  several  auxiliary  studies  of  other  parameters  were  also 
made  and  presented  among  these  figures  Except  where  it  is 


(x!>  yx)  origin  (  Q).  For  R o>^>  it  is  apparent  that 
In  other  words,  as  explained  earlier,  the  wavefront  appears  more 
spherical.  We  are  concerned,  however,  about  the  non-circularity 
of  the  wavefront  near  the  main  reflector*  and  Figures  13  and  14 
show  how  the  ray  directions  discussed  above  compare  in  this 
region  for  k  =  2tt  and  k  =  6 tt.  In  both  cases,,  we  find  that 
' —  ^2  over  the  maj°r  portion  of  the  180°  range.  Later  we 
will  find  it  extremely  valuable  to  make  the  approximation  that 
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rz  for  computational  purposes.  This  approximation 


will  be  most  necessary  near  ^80°  where  xi  , 

3  2,  £  2m 

since  this  is  the  general  area  where  the  edge  field,  has  the  largest 
effect  in  the  main  reflector  aperture.  It  is  seen  from  Figures  13 


and  14  that  the  approximation 
ion. 


<k  ~  C,  is  excellent  in  this  reg- 


We  now  turn  to  a  consideration  of  J'  (  fl)  as  a  function  of  fre- 

e  6 

quency  k.  It  is  clear  from  the  previous  discussion  that  j'  (  T-,) 

e  Lt 

will  differ  little  in  any  important  respect  from  J_i(  The 

results  of  this  study  are  plotted  in  Figures  15  through  19  for 
k  =  2tt,  4tt,  6tt,  8tt,  and  20tt,  respectively.  For  our  purposes 
in  this  section,  the  two  characteristics  to  take  note  of  are: 

1)  The  average  amplitude  away  from  the  shadow 
boundaries  diminishes  as  1/k  with  k„ 

2)  The  interference  effect  between  opposite  edges 
becomes  larger  in  amplitude  and  more  rapidly 
oscillatory  as  k  increases. 


From  the  point  of  view  that  we  would  like  to  eventually  deform 
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the  reflector  shapes  to  compensate  for  the  edge  field,  the 
first  fact  aids  us  and  the  second  fact  helps  to  defeat  this  purpose. 
It  is  impractical  to  attempt  to  correct  a  rapidly  oscillating  field 
by  optical  techniques.  However,  as  we  will  see  later,  the  field 
is  usually,  if  not  always,  negligible  in  those  regions  where  it 
oscillates. 

D.  Effect  of  Reflector  Shape  and  Field  Symmetry 

Before  proceeding  to  compute  the  actual  aperture  edge  field,  it 
will  be  interesting  to  study  the  effect  on  the  edge  field  distribu¬ 
tion  of  varying  the  primary  field  symmetry  and  the  subreflector 
shape . 

W  e  will  va.ry  the  prima.ry  field  symmetry  by  varying  the  relative 
phase  of  the  top  and  bottom  edge  sources.  That  is,  we  will  let 
SIGNEG  =  -1  instead  of  +1,  as  in  Standard  Set  Noc  1.  This 
type  of  primary  field  symmetry  occurs  quite  often almost 
enough  to  justify  a  complete  set  of  numerical  data  for  this  case.. 
For  the  most  significant  example,  we  refer  to  the  circularly 
polarized  primary  distributions.  In  addition,  monopulse  cr 
tracking  feeds  have  this  symmetry.  As  indicated  in  Figure  20, 
however,  the  only  significant  difference  between  the  field  of  this 
symmetry  and  the  former  is  that  a  zero  exists  at  P7  =  0°.  The 
remaining  general  characteristics  are  the  same  except  that  the 
oscillations  are  essentially  180°  out  of  phase. 

We  have  already  observed  the  diffe  reuctus  between  the  edge  field 
of  aflat  reflector  and  a  curved  reflector.,  The  d'lft'e rences  were 
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V. 


significant.  In  Figures  21  through  24  we  have  plotted  the  edge 
diffracted  field  as  we  varied  x2  ,  keeping  all  other  parameters 
fixed.  This  change  of  in  effect  changes  the  "tilt"  of  the 

subreflector  edges.  Referring  to  Figure  1,  this  change  is  mani¬ 
fested  by  a  change  in  As  x£m  increases,  ^  decreases. 

Therefore,  the  angular  position  of  the  shadow  boundary  changes. 
Of  course,  the  first  order  optics  field  also  changes.  On  Figures 
21  through  24,  the  position  of  the  shadow  boundary  is  indicated 
by  an  arrow.  The  change  in  position  of  the  shadow  boundary  and 
the  edge  field  is  very  small  over  a  range  of  x~  from  x~  = 
10  to  =  40.  This  range  covers  a  majority  of  practical 

cases  and  allows  us  to  draw  the  conclusion  that  the  choice  of 

x-  /x,  =  10  for  Standard  Set  No.  1  is  sufficiently  general 

2m  lm  7  & 

for  practical  purposes. 


The  Edge  Field  Aperture  Effects 


In  this  section  we  will  investigate  the  direct  effect  of  the  edge  diffracted 
field  in  the  main  reflector  aperture.  The  edge  field  is  reflected  by 
the  main  reflector  into  the  aperture  wherein  its  amplitude  per  unit 
area  will  be  determined.  To  assess  the  importance  of  the  aperture 
edge  field,  we  must,  of  course,  compare  it  to  <0  ^  or  1^  ,  the  aper¬ 
ture  field  generated  by  the  primary  source. 


We  will  find  that  the  edge  field  is  quite  significant  for  subreflector 
diameters  of  less  than  six  wavelengths  and  of  diminishing  significance 
for  greater  diameters .  These  results  are  in  substantial  agreement 
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with  those  found  from  general  engineering  practice  over  the  years. 
Without  presenting  any  statistics  to  document  this  statement,  we  will 
nevertheless  state  that  the  vast  majority  of  dual  reflector  systems  have 
subreflectors  greater  than,  usually  much  greater  than,  five  wave¬ 
lengths  in  diameter.  The  results  of  this  section  it  is  hoped  will  clar¬ 
ify  in  a  quantitative  way  the  reason  for  the  five  wavelength  diameter 
minimum  that  appears  in  standard  engineering  practice. 


In  order  to  calculate  the  aperture  edge  field  and1  tertiary  field  ( i) 


3p 


or  )  in  both  amplitude  and  phase,  we  must  choose  a  particular 
reflector  system.  As  was  pointed  out  earlier,  the  particular  system 
chosen  will  have  little  effect  on  the  substance  of  the  results.  We  will 
choose,  therefore,  a  Cassegrain  system  with  the  Standard  Set  No.  1 
parameters  as  the  boundary  conditions  for  most  calculations.  A 
parabola-hyperbola  system  is  much  easier  to  handle  analytically  than 
the  dual  reflector  designs  found  in  Chapter  1,  since  we  can  obtain 
closed  form  expressions  for  the  surface  contours  of  the  parabola- 
hyperbola  system.  Furthermore,  to  date,  this  system  is  the  most 
common  system  used  in  practice. 


In  this  section,  specifically,  we  will  first  define  the  parabola-hyperbola 
system  and  find  and  C^  (C  =  constant  for  a  Cassegrain  or 

Gregorian  system;  see  Chapter  1  for  example).  We  will  then  determine 
the  aperture  edge  field  and  present  extensive  numerical  calculations. 


A.  The  Cassegrain  System  and  the  Aperture  Field  Due  to 
the  Primary  Source 


The  main  reflector  of  the  Cassegrain  system  is  depicted  in  the 
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sketch  following  (i)  and  defined  by  the  quadratic  equation 


2  „ ,  .  2 

x2  =  -  4f2  ^2  +  X2m 

^2  =  (  X2m/  4f2)  -  (1/4f2>X2  =  a  +  bx2‘ 


(49) 


The  quantity  is  the  focal  length.  The  focus  depicted  below 
will  be  a  focal  point,  of  course,  for  the  hyperbola  subreflector 
also.  The  slope  at  any  point  on  the  main  reflector  is  given  by 


dy2  j  dx^  =  2b  .  (50) 

The  focal  length  is  found  from  defined  parameters  as 

f2  -  1/2  (\)  +  p1  +  \[  (  )) '  +  P j ) ^  +  xZm  )  (51) 

where  we  define 

])  =  cL*  +  ^3  -  ^3  ~  constant  (  ch  £■  0).  (51a) 

(See  equation  (58)  for  the  definition  of  p^  =  f^-k^.) 


The  hyperbola  geometry  is  depicted  in  the  sketch  (j)  following 
sketch  (i).  With  f ^ ,  the  distance  between  hyperbola  focii,  deter¬ 
mined  from  the  given  parameters  by 

q  =  9/2  ( - ^£0 - )  (52) 

x2m  "  x lm 

and  the  distance  in  the  sketch  found  as 

kj  =  1>-  fr  (53) 


we  define 
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kk  *  <k5  +  -L  +  f i  >  A 


■  \ 


-J 


k2  -  k2  f2 

\  l  l 


(54) 


and 


,2  ,2  2 

d  =  f1  -  c 


(55) 


and  obtain  the  equation  for  the  hyperbola  surface  as 


*1  s  ±  M 


or 


y,  =  k,  ±  \|  c  + 


— - 1 

d  ,  ,  v2  ,2 

(Yl  -  kl>  -d 

■  2 


2  ,  c2'  2 


j2  1 


(56) 


(57) 


The  sketch  (k)  given  on  the  following  page  of  the  parabola 
and  hyperbola  together  illustrates  the  composite  system.  Note 
that  0  for  this  system.  This  system  is  precisely  graphed 

in  Figure  25.  This  system  gives  a  focal  length  over  diameter 
ratio  for  the  parabola  of  a  little  greater  than  0.3,  a  common 
ratio  found  in  practice. 

With  kj  defined  for  the  hyperbola,  we  can  define  and  compute 


Pi  =  fi  -  ki .  (5S) 

The  maximum  value  of  y^  in  this  system  is  easily  found  as 

y  =  k,  -  c.  (59 ) 

7mx  1  '  '  ' 


Before  proceeding  to  determine  the  tertiary  field,  it  will  be 
interesting  to  find  the  secondary  field  or  the  geometrical  optics 
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field  reflected  by  the  hyperbola.  We  will  compute  this  numer¬ 
ically  for  the  standard  primary  field  1^(0^)  -  cos  ^  0^.  The 
actual  numerical  computations  are  made  by  assigning  a  value  to 
yj  and  computing  x^  from  (56).  We  then  find 


=  arc  tan  (- 


X_L 


l)+  Yi 


) 


(60) 


and 


=  arc tan  ( 


xi 


Yl  +  Pi 


(61) 


We  now  need 


d  0 


in  order  to  evaluate  (0-,)  from 

d  02  P  Z 


1^(0^)  by  application  of  the  energy  conservation  principle..  This 
derivative  is  obtained  by  first  finding 


d  X] 
dy, 


{72]  (  x 


Yl  -  M 


(62) 


d  0 


dy 


i-  = 


cos^  0 


dx_I  +  _ 


L(  j)  -  Yj>  dyj  (  y-YiT  J 


(63) 


and 


diz 

dy  1 


2  Q 

cos  0 . 


1  d  xi 

(y !  +  px)  ^y~\ 


(Y!  +-P1!)2  J 


(64) 


so  that 


d0!  =  dQu  / 

dyi  /  dyi 


d  0 


(65) 


With  (7  2)  we  apply  the  conservation  of  energy  principle  and  find 


i,  (e>)  = 


(  sln  0i. 

sin  02 


(dil 

d02 


(0,) 


(66) 


The  numerical  results  for  Standard  Set  No.  1  parameters  are 
depicted  in  Figure  26. 
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We  now  wish  to  obtain  the  tertiary  field  1^  .  This  can  be 
obtained  from  which  is  already  determined  by  (66)  if 

we  find  the  derivative  ('d  02/  dxj)  and  apply  the  conservation 
of  energy  rule  between  the  two  fields.  In  the  parabola-hyperbola 
system,  we  should  note  that 


x2  =  xj  .  (67) 

We  will  allow  to  be  the  independent  variable  and  assign  its 

value  first  and  then  determine  y^  ,  With  y^  determined,  we 
may  find  1^  as  previously  shown. 


Given  we  find  from  (49)  and  (50)  the  values  of  y^  and 

dy^/dx^  respectively.  Then 

tan  9 £  ~  (  x2/(yz  +  V  +  Pp)  (68) 

is  found.  We  solve  (56)  and  (57)  simultaneously  to  obtain 
Yj  ('tan.62)  given  by 


y^  =  (k^  d2  +  Pj  c2  tan2  02)/(d2  -  c2  tan2  ©2).  (70) 

and 

.,2,2  2  2  2  2,  2q  .  /, 2  2+  2  /7  n 

y  =  (d  kj  -  c  d  -  Pj  c  tan  02y(d  -c  tan  02)-  (71) 


In  (69)  we  use  the  +  sign  if  (yqn/yqd)  ^0  and  the  -  sign  other¬ 
wise,  where 

.2.2'  2  2 
y  =  k.  d  +  p,c  tan  0, 
y qn  1  r 1  2 
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and 


yqd  =  d  -  c“ tan  9z 

If  y  ^  is  approximately  zero,  then  use 


yl  =  0/2>  (yqin  /  V 


where 

AZ  .  2  2,2  2  2  2ft 

y  =■  d:  k,  -  c  d  -  P  t  c  c  an  0?  . 

qin  1  1  * 

Having  found  y^  we  obtain  as  discussed  already.  We 

still  need  dx^dB^  =  d  >^2,/ ^  ^ Z  to  °ktain  the  tertiary  field. 
From  (7  6)  we  obtain  this  value  as 


de . 


cos  0: 


dxz  (y2  +  +  Pi)  L 


1  -  tan  9  2 


(72) 


By  applying  the  conservation  of  energy  principle 


I_  (0_)  sin  0.  d0.  =  H  (x0)  x0  dx. 

2p  2'  2  2  3p  2'  2  2 


(73) 


we  obtain 


X3d  (x2}  =  X2d  <02>  2 L) 

3P  *  2p  2  X2  dx2 


(74) 


The  primary  fie]d  aperture  phase,  C  ,  is  constant  and  may  be 

P 

obtained,  as  a  check,  by 


C  =  r,  + 
P  1 


f2  *  f: 


-2  +  C 

3  Po 


(75) 


where 


r  ^  =  x  j  /s  in  0  j  , 


(76) 


and 


0  by  definition. 


C 

Po 

Numerical  results  are  depicted  together  with  results  for  the 
aperture  edge  field  to  be  discussed  next. 

Edge  Field  in  the  Aperture 

We  now  wish  to  determine  the  edge  field  amplitude  and  phase 

(I^e  or  i)^e  and  Ce)  in  the  main  reflector  aperture.  We 

will  determine  and  C  at  xi  *=  x-,,  the  same  point 

3e  e  Z  Z 

for  which  I0  and  C  were  determined.  This  will  facilitate 
3p  p 

the  complex  addition  of  the  two  fields  in  the  aperture  which  will 
be  made  in  Chapter  5. 

It  will  simplify  the  computations  considerably  if  we  assume 
that  r2  d2  or  that  the  true  ray  normal  to  the  reference 
phase  front  of  the  edge  diffracted  field  is  in  the  same  direction 
as  the  ray  through  the  (x^,  y^)  origin.  The  justification  for  this 
approximation  in  terms  of  its  accuracy  was  presented  earlier 
in  Figures  13  and  14.  It  should  be  noted  that  we  will  use  the 
true  phase  of  the  edge  field  at  the  given  point  on  the  wavefront. 

It  is  only  the  direction  of  the  ray  through  the  reference  wave- 
front  which  is  approximated. 

The  projection  of  the  approximate  edge  field  ray  as  it  is  reflected 
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from  the  main  reflector  into  the  aperture  at  the  point  x*  is 
illustrated  in  the  diagram  (1)  on  the  following  page.  Note  the 
difference  between  this  diagram  and  Figure  2,  where  the  ray 
is  in  the  (j^  direction.  The  ray  reflected  from  the  main  reflec¬ 
tor  is  in  the  ^  direction.  The  point  of  incidence  of  the  edge 
field  ray  on  the  main  reflector  is  designated  (x^,  y^)  in  the 
same  reference  system  as  (x^,  y^).  A  detailed  drawing  of 
this  ray  and  reflector  structure,  to  scale,  is  shown  in  Figure  27. 
The  equations  for  the  main  reflector  shape.  and  slope  in  terms 
of  (x^>  y^)-  is  the  same,  of  course,  as  that  given  in  equations 
(49)  and  (50)  for  the  (x^»  y^>)  system, 


x0  =  + 


( 


4f2  ?3  +  X2m 


^3 


=  (~^m  -)  -  (— — >  X, 


4f2 


4  £2 


and 


d  x^ 


tnr*  X3 


(7  9) 

(80) 


(81) 


Given  the  value  of  the  independent  variable  x7  =  x.-,  we  would 
like  to  find  P,  directly.  It  is  difficult  to  do  this  directly  so 
we  will  find  when  given  f 2>  From  the  geometry  we  have 


tan  P2  =  x3/(y3  +  l)  ) 

from  which,  when  solved  simultaneously  with  (80),  we  obtain 


(82) 


The  values  of 


-  2f^  ctn  r2  + 
y2  and  dy^/dx 


»4  ctn  P2  f2  *  4£2  x2m  *  (83) 

3  are  then  found  from  (88)  and  (81). 
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With  Snell's  law  at  the  point  (x3>  y^)  expressed  as 

dy^/dx^  =  -tan  (  — — — 7-...,u3  ) 
and  (84)  and  (81)  solved  simultaneously,  we  obtain 

-  2  arc  tan  (x^/zf^). 

From  (85)  and  the  geometry  we  then  readily  obtain 


(84) 


(85) 


x 


2 


x3  +  y3  tan 


(86) 


Now  given  x^  =  x^  we  may  solve  for  by  solving  the 

transcendental  equations  for  (81)  through  (86).  This  is  not 

very  difficult,  since  f. ^  (x^)  a  monotonic  function. 


We  now  need  dx^d  in  order  to  find  I3e  by  application 
of  the  conservation  of  energy  principle  to  Jg  or  (note  that 

the  assumption  that  (j^  ^  P>  implies  also  that  J  J^). 

This  value  is  found  by  differentiation  of  previous  results  as 


dx2 


dx 


d  r2 


?■  +  tan  T) 


dx3 ' '  d  r2 


(87) 


where 


d  <K 
d  fS 


■  *  -  ( 


COS 


+  y3  tea2  y  (7^ 


(88) 


and 
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2f- 


c  sc 


ft 


x3/J4ctn  r2f2 


4f2^+x2m 


•  (89) 


Applying  the  conservation  of  energy  principle,  we  obtain  fin¬ 
ally  that 


I 


3e 


sin 

x2 


(90) 


In  order  to  determine  the  edge  field  phase,  Cg,  at  the  same 
point  in  the  aperture  we  find 


p4  =  x3  csc 

p5  =  (x2  -  x3)  csc  (  <^), 

and  finally 

Ce  *  f>4  +  ^5  *  R' 

Note  that  C  ^-R  according  to  the  approximation 

K  ~  r2.e° 


(91) 

(92) 


(93) 


In  Chapter  5  we  will  use  the  synthesis  method  of  Chapter  1  in 
an  attempt  to  correct  the  aperture  edge  field  effect.  In  order 
tc  carry  out  the  synthesis,  it  will  be  remembered  that  the  der¬ 
ivative  of  the  desired  aperture  phase  distribution  must  be  found. 
For  the  work  in  Chapter  5.  we  will  need,  therefore,  to  find 

d  Cp  _  d  Cp  By  the  theorum  of  Malus,  we  find  directly 
d  x£  ~  d  x2 " 


that 
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d  Ce^dx7  =  sin 


where  (T  is  given  by  (85). 


(94) 


C*  Aperture  Edge  Field  Computations 

Aperture  edge  field  computations  of  and  G  were  made 

for  a  wide  range  of  frequencies,  for  normal  and  parallel  polar- 
izations  (  dP =  tt / 2,  and  dK =  0°),  and  for  even  and  odd  sym¬ 
metry  of  the  primary  field  (SIGNEG  =  +1  and  -l)  incident  on 
the  subreflector.  These  results  will  be  presented  graphically 
and  briefly  discussed  in  this  part.  Unless  it  is  specifically 
stated  otherwise,  the  Standard  Set  No.  1  parameters  introduced 
earlier  may  be  assumed, 


Figures  28  through  29  illustrate  the  edge  field  amplitude, 
contribution  to  the  aperture  field  in  comparison  with  the  primary 
source  contribution,  *0^*  for  frequencies  k  =  2tt,  4t r,  6tt, 

8t r,  IOtt,  16tt,  and  20tt.  The  primary  source  contribution  does 
not  change  with  frequency,,  of  course,.  The  important  general 
characteristics  of  the  edge  field  contribution  are  the  same  for 
all  frequencies,  and  as  will  be  seen  subsequently,  for  different 
polarizations  ($)  and  symmetries  (SIGNEG).  The  character¬ 
istics  are: 


1)  A  large  lobe  exists  near  the  outer  periphery  of  the 

main  reflector  *  Although  we  would  expect  this  lobe 

to  be  centered  around  x,  we  find  that  it  is  shifted 

2  m 

towards  a  value  of  slightly  less  than  x2m'  This 

shift  may  be  accounted  for  in  part  by  the  approximation 
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1*2  (see  Figures  13  and  14),  Some  of  the 

shift  is  due,  however,  to  a  change  in  the  amplitude 
distribution  caused  by  reflection  from  the  main  ref¬ 
lector,  For  a  smaller  value  of  k  or  a  wider  edge 
diffraction  lobe,  this  shift  is  less  important,  since 
it  is  approximately  the  same  for  all  values  of  k 
considered.  The  amplitude  in  this  diffraction  lobe 
region  has  been  calculated  with  the  shadow  boundary 
correction  discussed  earlier  in  this  chapter  and 
should  be  reasonably  reliable.  This  contribution 
is  of  major  importance  and  will  be  discussed  further 
later. 

2)  Away  from  the  shadow  boundary  lobe;,  the  field 

decays  rapidly,  the  rate  of  decay  increasing  with 
frequency.  In  the  shadow  boundary  lobe  region  and 
in  this  decay  region,  there  is  no  interference  effect 
between  opposite  edges  of  the  subreflector..  The  top 
edge  field  dominates  in  this  region.  It  is  possible, 
however,  to  include  some  contribution  from  the  bot¬ 
tom  edge  despite  the  optical  blocking  of  the  curved 
surface  of  the  subrefle^tor .  The  inclusion  of  sur- 
face  diffracted  rays  is  discussed  by  Keller  .  How¬ 
ever,  for  this  polarization  (  $  =  0),  the  transition  to 
the  region  where  interference  between  edges  takes 
place  is  so  smooth  that  it  is  clear  that  the  present 
theory  is  adequate.  For  the  other  polarization 
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(  d  ~  tt/2)  we  will  see  shortly  that  a  discontinuity 
does  exist  in  this  transition,  but  the  discontinuity 
is  of  no  qualitative  significance, 

3)  In  the  region  where  the  fields  of  the  two  edges  inter¬ 
fere  and  oscillations  do  take  place  x*  <:  6)  the 
field  amplitude  levels  off,  but  at  a  rather  low  level 
-20  db.)  with  respect  to  the  shadow  boundary  lobe 
and  This  is  significant  since  there  is  very 

little  possiblity  that  the  rapid  oscillations  in  this 
region  could  be  corrected  by  optical  techniques 
(i.  e.  ,  reshaping  the  reflector  surfaces).  However, 
the  low  amplitude  for  all  frequencies  in  this  region 
makes  a  correction  of  this  effect  unnecessary. 


From  the  figures,  then,  it  is  clear  that  the  only  significant  con- 
tribution  of  the  edge  diffracted  field  is  the  shadow  boundary  lobe 
where  the  amplitude  is  of  the  same  order  as  the  primary  source 
contribution.  Actually,  this  contribution  is  even  more  signifi¬ 
cant  than  is  apparent  from  the  figures.  First  of.  all,  the  area 
subtended  by  the  small  radial  segment  of  the  shadow  boundary 
lobe  is  quite  large,  since  this  effect  is  at  the  outer  edge  of 
the  main  reflector-  Secondly  and  even  more  significant  is  the 
fact  that  this  contribution  to  the  total  aperture  distribution  is 
usually  ^180°  out  of  phase  with  the  rest  of  the  aperture  field. 
This  effect,  which  is  demonstrated  quite  clearly  in  Cha.pter  5,  is 
particularly  important  with  regard  to  the  gain  of  the  antenna. 
Furthermore,  it  is  not  possible  to  simply  change  the  value  of  S 


1Q7 


in  a  new  parabola-hyperbola  antenna  system  in  order  to  cause 
the  edge  field  to  add  in  phase  instead  of  out  of  phase  at  the  aper¬ 
ture.  It  is  reasonably  clear  from  Figure  2  that  changing  the 
boundary  conditions  does  not  alter  the  relative  phase  between 
these  two  fields  (note  that  only  the  top  edge  makes  a  contribution 
to  the  shadow  boundary  lobe).  What  is  actually  required  to  make 
the  edge  field  contribution  such  that  the  resultant  aperture  phase 
distribution  is  approximately  uniform  is  to  alter,  fundamentally, 
the  reflector  contours.  This  is  done  in  Chapter  5  with  some 
success . 

Let  us  examine  a  little  more  closely  the  origin  and  character 
of  the  shadow  boundary  lobe.  In  the  brief  succession  of  diagrams 
(m)  is  depicted  the  degree  of  approximation  obtainable  with  larger 
and  larger  subreflectors  of  the  desired  geometrical  optics  pattern 
(the  dashed  line).  This  is  similar  to  using  more  and  more  terms 
of  a  Fourier  series  to  approximate  a  square  wave.  The  difference 
between  the  two  curves  at  the  vertical  edge  of  the  dashed  curve 
is  the  shadow  boundary  lobe.  In  general  the  amplitude  at  this 
point  does  not  vary,  but  the  steepness  of  the  approximating  curve, 
does  ,  and  hence  the  width  of  the  shadow  boundary  lobe  also  does 
vary  in  the  same  manner  with  increasing  k.  This  is  the  calculated 
result  illustrated  in  the  figures. 

For  k  =  2tt  we  find  that  the  shadow  boundary  lobe  is  significant 
for  approximately  9  <  x^  <  10  and  hence  10%  of  the  radius  or 
approximately  20%  of  the  area  of  the  aperture  since  the  lobe  is 
near  the  outer  periphery.  Since  this  lobe  is  out  of  phase  with 
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the  rest  of  the  aperture,  the  net  reduction  of  anticipated  gain 
for  the  antenna  may  be  as  great  as  40%,  surely  an  intolerable 
figure  in  practice.  For  the  case  of  k  =  6tt,  we  obtain  a  similar 
figure  of  about  15%  net  reduction  in  possible  gain.  Even  if  these 
figures  were  30%  and  10%,  qualitatively,  the  significance  is  appar¬ 
ent.  A  subreflector  of  about  6  A  in  diameter  (k  =  6tt)  is  the 
minimum  tolerable  size  that  can  be  used  in  practice.  This,  as 
was  mentioned  earlier,  agrees  with  generally  accepted  standards 
found  in  practice.  For  much  larger  subreflectors,  the  effect 
becomes  progressively  less  significant  as  may  be  expected. 

In  Figures  31  through  36  are  illustrated  the  edge  field  phase 
distributions  in  the  aperture,  Ce>  for  k  =  2t t,  4tt,  6tt,  and  8rr. 
Figures  31  and  32  illustrate  the  same  distribution,  k  ~  2tt:(  but 
in  Figure  3  1  the  various  braches  of  are  connected  to  illustrate 

the  continuous  distribution.  The  various  branches  are  presented 
segmented  in  most  of  the  other  illustrations.  The  change  in 
branch  for  arises  from  the  branch  choice  of  the  arctangent 

in  equation  (41)  for  R.  Since 

Ce  '  f4  +  P5  -  R 

this  branch  choice  is  reflected  in  the  results  for  C  also, 

e 

Each  branch  in  the  phase  curves  for  represent  another  lobe 

of  oscillation  in  the  amplitude  „  Hence,,  for  higher  values 

j  e 

of  k  there  are  more  oscillations  or  lobes  and  more  branch 
segments  in  the  curves  for  If  the  phase  front  of  the  edge 


u< 


field  were  perfectly  circular  with  the  center  of  curvature  at  the 

focal  point  of  the  main  reflector  parabola,  then  the  phase  curves 

for  ,  when  they  were  connected  continuously  as  illustrated 

in  Figure  3  1?  would  appear  as  shown  in  the  diagram  (n).  They 

wohld  be  perfectly  formed  squared  steps  of  t r  radians  or  X/Z 

height.  The  step  jumps  would  occur  at  zeros  of  For  a 

slightly  non-circular  edge  diffracted  field  from  the  subreflector 

or  one  with  the  center  of  curvature  not  at  the  focal  point  of  the 

parabola,  the  resulting  distribution  might  appear  as  shown  in  the 

diagram  (o).  In  this  case  the  amplitude  S)  does  not  have 

zeros  and  the  steps  of  are  not  clearly  delineated  as  before. 

An  excellent  example  of  C  for  this  case  is  illustrated  in 

Figure  48.  The  less  circular  the  edge  field  diffracted  wavefront 

becomes  the  smoother  the  curves  for  C  become. 

e 


In  Figures  37  through  39  are  illustrated  the  aperture  edge  fields,. 

,  when  -the  polarization  incident  on  the  edge  of  the  subreflector 
is  normal  to  the  subreflector  edge  (  tt/2).  This  may  be  inter¬ 

preted  alternatively  as  an  E-plane  cut  whereas  when  -  0,  the 
fields  then  are  in  the  H -plane  (planes  parallel,  to  electric,  E,  or 


magnetic,  H,  fields).  The  oscillations  in  the  region  where  the 
top  and  bottom  edge  fields  interfere  (x^^6)  are  considerably 
greater  than  for  the  0  case.,  While  they  are  still  small  com¬ 

pared  to  che  sb?do\v  boundary  lobe,  they  are  not  entirely  insignifi¬ 
cant,  especially  for  the  lower  frequencies,  k  =  Ztr  and  4rr.  At 
around  -  7.  7,  we  find  a  discontinuity  in  the  field  that  is 

due  to  the  sudden  introduction  of  edge  source  No.  2  (the  bottom 
edge)  into  the  computations*  At  this  point,,  the  surface  of  the 
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subreflector  no  longer  optically  blocks  the  rays  from  edge  No.  2 

and  the  interference  between  edges  commences  for  x^  less 

than  this  value.,  As  .mentioned  earlier >  the  creeping  or  surface 

7 

wave  that  can  be  excited:  by  the  edge  diffracted  field  and  "seen" 
around  the  curved  surface  of  the  reflector  is  much  more  signifi¬ 
cant  for  this  polarization.  While  such  a  correction  to  the  com¬ 
puted  values  would  give  a  more  continuous  transition  at  this  point 
(x2l  ^7.7),  this  effect  is  not  considered  very  important  for  the 
purposes  of  this  work.  The  effect  is  interesting,  however,  and 
indicative  of  the  general  validity  of  Keller's  theory"^  when  properly 
interpreted. 


In  Figures  40  through  42,  the  edge  field  phase  in  the  aperture, 

C  ,  is  illustrated  for  k  =  2tt,  4t t,  and  6ir  with  (J4  ,/2.  Note 
e 

that  this  phase  distribution  is  more  sharply  stepped  than  the 
equivalent  case.  As  discussed  earlier,,  this  result  is  con- 


sistant  with  the  lobe  structure  found  for 

0 


when 


/e 


as  compared  to  that  for 


and  implies  that  a  more  spherical 


wavefront  for  the  edge  diffracted  field  exists  here. 


Figures  43  through  48  illustrate  the  amplitude  and  phase  for  the 
case  when  odd  symmetry  (SIGNEG  ^  -1)  exists  in  the  primary 
source  distribution.  Both  polarizations,  0  and  rr/2,  were 

used  in  computations  for  k  =  6rr.  There  is  no  difference  of  any 
consequence  for  the  SIGNEG  =  - 1  case  relative  to  the  SIGNEG  ~  +1 
case.  A  zero  for  <J^)  exists  at  x}  »  0  for  SIGNEG  =  -1  whereas 

jC  iw 

a  maximum  exists  at  this  point  for  SIGNEG  =  +  l,  but  this  expected 
result  of  the  symmetry  is  of  no  important  consequence.  The 


same  differences  between  the 


<y=  0  and 


=  tt/2  cases  exist 
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for  SIGNEG  =  - 1  as  do  for  SIGNEG  =  +1.  This  is  true  for  both 
the  phase  and  amplitude  results  as  is  clearly  evident  from  the 
figures.  No  deviation  of  these  results  need  be  expected  for  other 
values  of  k. 
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CHAPTER  4 


ON  THE  POSSIBILITY  OF  EXTENDING  THE  SYNTHESIS 
METHOD  TO  INCLUDE  THE  EDGE  DIFFRACTED  FIELD 

Introduction 

The  discussion  and  computations  presented  in  Chapter  3  illustrated 
the  importance  of  the  effect  of  the  edge  diffracted  field  upon  the  aper¬ 
ture  distribution.  The  synthesis  method  developed  in  Chapter  1 
showed  how  a  specified  aperture  amplitude  and.  phase  distribution 
may  be  obtained  by  shaping  the  reflectors'  c  ro  ss  -  sec  tions  when  the 
diffracted  field  is  not  included,,  It  would  be  ideal  if  we  could  extend 
the  synthesis  method  of  Chapter  1  so  as  to  include  the  edge  diffracted 
field's  contribution  to  the  aperture  distribution.  In  this  chapter  we 
will  discuss  an  attempt  at  finding  such  a  solution  which  proved  unsuccess¬ 
ful,  The  method  for  obtaining  a  solution  which  was  attempted  appears 
to  be  a  powerful  approach  to  the  problem  despite  the  lack  of  success 
found  in  using  it.  The  inadequacy  of  this  approach  was  by  no  means 
proven  so  it  seems  worthwhile  to  present  the  technique  for  the  record* 

The  use  of  the  edge  source  on  the  sub -r eflector  for  obtaining  the  dif¬ 
fracted  field  is  an  excellent  choice  for  several  reasons.  Besides 
providing  an  accurate  result  for  the  diffracted  field  of  small  sub- 
reflectors  as  was  discussed  in  Chapter  3,  the  edge  diffracted  field  is 
independent  of  the  reflector  cross-section  and  depends  only  upon  the 
boundary  conditions  such  as  were  stated  in  the  synthesis  method  of 
Chapter  1.  These  conditions  are  the  position  of  the  edge  and  the  slope 

of  the  edge  which  are  determined  by  X,  ,  X0  ,  OC  (Q,  V,  and  8(6,  ). 
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Hence  the  diffracted  field  is  a  fixed  field  distribution  which  illumina¬ 
tes  only  the  main  reflector.  Since  it  does  not  depend  upon  Y^(X^)  as 
a  Kirchhoff  integration  would,  it  lends  itself  to  an  extension  of  the 
synthesis  method  of  Chapter  1. 


It  is  not  at  all  clear,  however,  whether  a  solution  exists  to  this 
problem,  Unfortunately,  this  problem  will  remain  unresolved  in 
this  work,  although  it  is  the  feeling  of  the  author  that  the  results  found 
in  this  chapter  imply,  by  no  means  prove,  that  a  solution  does  not 
exist  under  any  conditions  except  that  J0  =  0. 


It  is  well  known*  that  an  arbitrary  phase  and  amplitude  cannot  be 

specified  in  the  aperture  of  a  single  reflector  system  for  an  arbitrary 

source  illumination.  Referring  to  Figure  1,  we  see  that  and 

must  be  controlled  by  shaping  both  Y^(X^)  and  Y^(X^)  in  order  that 

upon  complex  addition  to  and  C^,  a  specified  amplitude  and  phase, 

1 2  and  C,  will  result  in  the  aperture.  From  an  iterative  point  of 

view,  whenever  I_  and  C  are  altered  to  compensate  for  I0  and  C  » 
ip  p  ^  3e  e 

the  change  in  Y^(X^)  and  Y^(X^)  also  alters  1^  and  C^.  It  is  not 
clear  then  whether  an  iteration  would  converge.  The  results  of  Chapter 
.5  indicate  that  such  an  iteration  would  not  converge,  although  by  no 
means  conclusively.  If  a  solution  does  not  exist  of  course  such  an 
iteration  would  not  converge.  The  inverse  is  not  necessarily  true, 
however,  since  the  problem,  mathematically,  is  not  a  linear  one.  if 
were  very  large,  it  is  clear  that  a  solution  would  not  exist.  The 
dominance  of  I^e  over  1^  would  preclude  the  possibility  of  correction 
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of  I^e  by  adjustment  of  .  On  the  other  hand,  even  when  I^e  is  very 

much  less  than  I ^  ,  it  is  not  clear  whether  or  not  a  solution  exists. 

3p 


The  above  method  of  reasoning  was  inconclusive  so  that  a  more  formal 

or  mathematical  approach  to  the  problem  is  necessary.  Before  actually 

attempting  the  solution  mathematically,  let  us  see  from  Figure  2  how 

the  various  laws  and  principles  relevant  to  the  problem  fit  together  . 

A  normal  congruence  for  the  rays  and  wavefronts  is  defined  for  the 

primary  and  edge  sources  when  C^o  and  Ceo  are  specified.  We  know 

that  this  normal  congruence  is  maintained  upon  successive  reflections 

if  the  appropriate  optical  principles  are  applied.  Application  of 

Snell’s  law  at  reflector  #1  insures  a  normal  congruence  for  the 

waves  incident  upon  reflector  #2.  Now  the  theorum  of  Malus  applied 

for  C  and  C  in  the  aperture  insures  a  normal  congruence  foi  the 
p  e 

waves  and  rays  reflected  from  the  main  reflector.  Hence  it  would 

appear  that  Snell's  law  applied  at  the  main  reflector  would  be  redundant. 

However,  the  two  energy  conservation  equations  taken  together  with 

those  principles  already  stated  do  not  define  a  unique  solution  to  the 

problem  for  Y^(X^)  and  Y^X^).  Since  and  are  not  known  (and 

neither  are  C'  =  dC  /dXf  nor  C  '  =  dC  /dXl),  then  it  may  be  true 

p  p  c  e  e  £ 

that  inclusion  of  Snell's  law  at  the  main  reflector  is  not  redundant. 
Inclusion  of  both  these  principles  as  independent  conditions  would  mean 
that  the  problem  is  overstated;  that  is,  there  is  no  solution.  On  the 


other  hand,  since  C  and  I_.  are  given  functions,  it  may  be  that  one  of 
the  two  Snell's  law  applications  is  necessary,  whereas  the  other  is 


redundant.  This  reasoning  is  certainly  inconclusive,  but  unfortunately 
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that  is  the  state  of  this  attempt  to  an  exact  solution  to  this  problem.  A 
reasonably  satisfactory  approximate  solution  to  this  problem  is  obtained 
in  Chapter  5.  The  remainder  of  this  chapter  is  devoted  to  explicitely 
formulating  the  above  ideas  mathematically. 


II.  Equations  for  the  Attempted  Exact  Solution 

The  mathemtaic al  formulation  of  the  problem  described  in  Section  I 
above  will  be  performed  briefly  in  this  section.  We  will  attempt  to 
arrive  at  a  set  of  six  simultaneous  first  order  differential  equations 
with  the  dependences  indicated  below: 


'i- 

£  yz (  K  Cp'  Ce’  X2>| 

(i) 

/From  Energy  Conservation 

•;  ■ 

£&1  (0r  Cp’  Ce’  X2J  j 

(2) 

O 

T3  - 

II 

(Cp«V  9r  Vj,  Y2;X^  ' 

(3) 

1 

From  Theorum  of  Malus 

1 

c  = 

e 

fCe  |Ce'  <S2;  X2> 

1 

(4) 

v'  = 

1  1 

f Y i  (0i«  Ce’  Cp’  Yl*  Y2;  X2^  \ 

(5) 

1 

y  From  Snell's  Law 

Y'  = 
2 

£Yz  (0,,  Ce,  Cp.  Yj.  Y 

2’  X2> 

(6) 

/ 


These  equations,  as  shown,  form  a  complete  interdependent  set. 

i 

The  primed  variables  denote  total  differentiation  with  respect  to 
except  in  the  case  of  X^.  The  total  derivative  of  X^  with  respect  to 
will  always  be  denoted  as  dX2/dX^.  Rather  arbitrarily,  all  the 
variables  may  be  considered  dependent  variables  with  X^  the  independent 
variable.  The  remaining  variables  of  significance  in  the  problem  can 


169 


be  found  from  those  in  (1)  through  (6)  with  the  following  dependence: 


X3  =  X3(Ce,7f2;  X*)  (7) 

Y 3  =  Y3(Ce')fz-'  Xz>  (8> 

X2  =  X2(Cp,01,.YrY2;  X-)  (9) 

Xr  =  Xjiej.Yj)  (10) 


As  we  derive  the  explicate  forms  of  the  above  equations,  we  will 
point  out  the  deficiencies  in  the  formulation. 


The  energy  conservation  equations  leading  to  (1)  and  (Z)  are 

^2  =  ^e/^e  ^2^  j  [X2////^in  ^2 

and 


The  quantities  I^e  and  1^  can  be  expressed  in  forms  of  the  specified 
aperture  distribution  amplitude,  1 ^  ,  and  phase,  C,  as  will  be  done 
shortly. 


From  expression  of  the  theorum  of  Malus  at  the  aperture  we  obtain 


and 


C ' 
e 


sin  2f3 


(X* 


(14) 


which  will  lead  to  equations  (3)  and  (4)  . 


Snell's  law  for  each  reflector  and  each  ray  system  is  expressed  as 


tan 


r  -| 

1 

lxz  -  X1  1 

h  -92 

X]  =  tan  J 

91  ■  atcta"(imp?1j 

’ 

x;  (is) 

1  r 


Y2  =  "tan  7 


e2  -  63 


dX 


dX!, 


=  -tan 


X. 


X, 


arc tan 


v+ Y2-  y. 


/  X2  -  X2 
arc  tan  | - — 


dX2 

dX' 


(16) 


and 


1 

- 

1 

/  X3  \ 

=  -tan  ^ 

2f  2  '  &  3 

X3  —  ~tan 

arctan  hr+a- 

\  3  e/ 

X2  -  X3\ 


Y3  ! 


(17) 


Equations  (15)  and  (16)  lead  to  equations  (5)  and  (6).  It  is  here  that  the 
difficulties  or  unresolved  portions  of  the  problem  exist.  First  of 
all,  as  will  be  seen  shortly,  it  turns  out  that  equation  (17)  is  not 
necessary  to  arrive  at  the  formulation  expressed  by  equations  (1) 
through  (10).  This  may  be  explained  by  stating  that  equation  (4)  or 
(14)  for  the  Theorum  of  Malus  makes  (17)  redundant.  In  this  case, 
however,  we  would  expect  that  (16)  would  also  be  redundant  with  (3) 


or  (13).  However,  it  is  not  clear  that  (1)  through  (6)  are  not  a 
complete  set  of  differential  equations  with  a  unique  solution  unless 
(6)  or  (16)  is  included.  Excluding  (6)  leads  to  the  conclusion  that 


^2^2^  arbitrary*  It  is  conceivable  that  this  is  the  case,  and  that 
the  Cauchy -Lipschitz.  conditions  (see  Appendix  1)  are  not  satisfied 
except  for  certain  Y^X^)  and  in  some  cases  for  no  function  Y^(X^). 
This  last  conclusion  is  the  only  one  which  appears  consistent  with 
all  the  results  obtained.  It  is  however  very  unsatisfactory  from  an 
engineering  point  of  view  since  appropriate  choices  of  Y  ?(X'^)  are 
too  difficult  to  find  and  their  existence  is  not  certain  in  any  specific 
case.  The  general  complexity  of  the  equations  precluded  the 
possibility  of  applying  the  Cauchy-Lipschitz  conditions  in  any  useful 
way.  In  addition,  as  the  results  in  Chapter  3  have  shown,  we 
frequently  find  I^^j;  I^p  for  some  X^  and  may  expect  that  very  often 
a  solution  will  not  exist  in  practical  cases.  The  formulation  is 
however  uncertain  and  the  possibility  remains  that  a  solution  does 
not  exist  even  when  . 


Continuing  the  formulation,  we  find  for  the  unknown  path  lengths  or 
phases  at  the  aperture  for  each  source,  the  primary  and  edge  sources, 
that 


C 

P 


(/-Yjl/cos  Oj  4  y(Xz  -  XL)Z+  (V+Y^Y^ 

+V'X l  -x/  +y|  '  +  cpo<c,> 


where 
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Xj  =  (^-YjJtanOj,  (19) 

and 

ce  =  (Y3  +  de(y2 ))/cos«'2  +V(X2  '  X3>+^+  Ceo(y2)  (20) 


where 


X3  =  (Y 3  +  de  (tf2))  tan^. 


(21) 


The  inclusion  of  C  and  C  in  these  expressions  defines  the  primary 
po  eo  r 

and  edge  source  rays  and  wavefronts  to  be  a  normally  congruent 
system . 


If  we  let^(Xi,)  =  I3a>^3p  -  I3  *  and^^  =  *3e  '  then  we  can  add  the 
two  fields  in  the  aperture  by  complex  addition, 


J|exp(kc)  =i)3p  exp(kCp)  +  ^3eexp(kCe),.  and  find  that 


bp  =  ba(X2>  -T7 

*  cin 


sin  k(C(X2)  -  C  ) 


(22) 


sin  k(C  -  C  ) 
'  p  e' 


and 


be  =  ha<X2> 


sin2k(C(X^  -  Cp) 
sin2k(Cp  -  Ce) 


(23) 
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We  will  now  outline  the  substitutions  that  are  necessary  to  arrive  at 
equations  (1)  through  (10)  from  (11)  through  (2  3)  and  in  addition 
point  out  some  c£  the  difficulties.  In  the  following  we  read  "(11) 


+  (23)*^  "as"  equations 

(11)  and  (2  3)  imply  .  11 

The  plus  sign  does 

not  necessarily  indicate 

summation.  The  reductions  follow: 

^(11)  +  (23)  (1) 

(24) 

[(12)  +(22)  (2) 

(25) 

[(18)  +  (19)  ^(9) 

(26) 

[(2  0)  +  (21)  -*>  (7) 

(27) 

[(20)  +  (21)  (8) 

(28) 

[(26)  +  (13)  (3) 

(29) 

[(27)  +  (28)  +  (14)-* 

(4)] 

(30) 

[(25)  +  (19)  ==>  XJ 

.  "DXl  f  ,  ?>Xl 

T,91  01  +  £Y2 

d 

(31) 

[  15)  +  (19)  +  (26)  +  (31)  -=>(5)1 

(32) 

1(26)  +  (29)  +  (25)  +  (32)  dXZ/dxl  =^2  C 

l  X  2  0Cp  p 

3X2  VXl  ?>X2  2x>  1 

+  Wi  1  +  QY i  Yi  +  3X1  *-W~z  Y2 ] 

(33) 

and 


(16)  +  (19)  +  (26)  +  (33)  j-|^(6) 


(34) 


With  the  exception  of  (33)  and  (34),  no  difficulty  with  the  above 
substitutions  occurs  and  they  lead  directly  to  explicite  forms  of  (1) 
through  (5).  Upon  substituting  (33)  in  the  expression  for  Yi,  ,  equation 
(16)  ,  and  solving  for  Y-,  as  directed  in  (34)  we  find  that  Y^  is  identi¬ 
cally  zero.  Such  a  result  does  not  make  physical  sense  and  no 
alternative  or  method  for  finding  a  usable  explicite  form  of  (6) 
has  been  found.  It  is  possible  that  this  result  is  due  to  a  redundancy 
of  Snell's  law  and  the  Theory  of  Malus  as  discussed  earlier.  In 
this  case  we  are  left  with  the  indefinite  results  and  conclusions  drawn 
earlier.  We  will  therefore  conclude  the  discussion  of  this  attempt 
at  an  exact  solution  by  outlining  more  clearly  the  derivations  indicated 
in  (34). 


By  performing  the  operations  indicated  in  (34)  we  obtain 

t>x2 


Y i  =  -tan 


9  2  '  93 

- 2 - 


1  +  tan 


$X2 

*X2 

+  eiW7 

I °Z  -e3 

\  ?xz  ' 

l  ^ 

1  ®V2J 

+  Y 


Z>*2 

i  W, 


(35) 


We  then  obtain  equation  (9)  from  (18)  and  (19)  as  directed  by  (26). 
Differentiating  (9)  then  gives  the  partial  derivatives  necessary  for 
for  evaluating  (35). 
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Consider  the  denominator  of  (35)  first.  Evaluating  we  get 

0Y2 


£X2  (co  s  0^  +  cos  0^) 

£Y2  ( sin  ©2  -  sin  9^  ) 


(36) 


l»z  -  e3 


We  find  therefore,  upon  expanding  tan!  ^  I  >  that 


jl  +  tan  | 


92  -  93  V  ^X2 


=  1  - 


sin(92  -  9^)  \  I cos©2  +  cos93 


1  +  cos(92  -  9^) 


\  lcosQ 2  + 

rNyrr 


>In9 


3 

(37) 


and  that  therefore  the  denominator  of  (35)  goes  to  zero  if  either  ©2  =  0 
or  ©3  =  0.  We  will  find  however,  that  the  bracketed  terms  in  the 
numerator  vanish  for  all  ©2  and  9^. 

From  equation  (13)  we  have  that 

C  1  =  sin  9.,  . 

P  3 

Differentiating  (9)  gives 


3X2 

*7 


(sin  ©2 


sin  9^) 


(38) 


and 


^X2  -sin  9 3  S' 

QX'Z  {sin  92  -  sin  9^), 


(39) 


17 


so  that 


t>X2 


C' 

p 


2>x2 
+  7*? 


5  o. 


This  leaves  (35)  as 


With  considerable  algebraic  effort  for  the  case  wherein 


the  case  of  a  fixed  center  of  phase  for  the  primary  source,  we  can 
show  that 


t>*2 


+  Y  \ 


Wz 

^r 


5  0 


for  all  values  of  the  variables. 


dC 


and. 


TG" 


po 


are  not 


independent  if  a  normally  congruent  ray  and  wavefront  system  is 
defined  for  the  primary  source.  It  is  presumed,  although  not  shown, 


that  the  same  result  would  be  obtained  with 


d/& 


^  o  and 


dC 


po  ^ 


dQ, 


0. 


The  case  when  these  quantities  are  zero  is,  however,  by  far  the  most 
important  practical  case.  Differentiation  of  (9)  gives 


1?7 


and 


ax2 

wr 


sec  0 


1 


sinO^tanQ^  +  cosG^ 


sinO^) .. 


(42) 


These  equations  together  with  (40)  may  be  used  to  verify  the 
above  assertions.  In  order  to  complete  the  result,  the  explicite  form 
of  and  Yj  must  first  be  obtained  by  carrying  out  the  substitutions 
indicated  in  (2  5)  and  (32). 


If  the  difficulties  indicated  in  the  above  derivations  could  be  overcome 
this  method  for  obtaining  an  exact  solution  would  be  very  powerful 
since,  for  one  reason,  the  six  simultaneous  first  order  differential 
equations  (1)  through  (6)  are  readily  solved  by  the  high  speed  computers 
available  today.  However,  since  in  this  work  the  problem  of  the  edge 
source  diffracted  field  has  been  defined  and  a  partial  correction 
obtained  (see  Chapter  5)  no  further  effort  for  an  exact  solution  will 


be  made. 
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CHAPTER  5 


CORRECTING  THE  EDGE  DIFFRACTED  FIELD 
BY  APPLICATION  OF  THE  SYNTHESIS  METHOD 

Introduction 

The  effect  of  the  aperture  edge  field  and  its  significance  was  dis¬ 
cussed  in  detail  in  Chapter  3.  As  was  evident  from  the  results  of 
Chapter  3,  the  shadow  boundary  lobe  could  be  reduced  to  a  negligible 
width  by  increasing  the  size  of  the  subreflector  in  wavelengths.  This 
approach,  however,  merely  substitutes  aperture  blockage  for  the 
aperture  edge  field  effect.  To  date  this  has  been  the  standard  approach 
to  this  problem.  In  this  chapter,  we  shall  show  that  much  of  the  deleter 
ious  effect  of  the  shadow  boundary  lobe  can  be  reduced  by  application 
of  the  synthesis  method  developed  in  Chapter  1. 

It  was  asserted  in  Chapter  3  that  the  most  significant  aspect  of  the 
shadow  boundary  lobe  effect  was  the  fact  that  it  causes  the  total  field 
in  the  aperture  in  the  region  of  the  lobe  to  be  180°  out  of  phase  with  the 
rest  of  the  approximately  uniformly  phased  aperture.  Computations 
will  be  presented  later  to  verify  this  assertion.  The  results  obtained 
in  this  chapter  by  repeated  application  of  the  synthesis  method  is  that 
the  shadow  boundary  lobe  phase  can  be  made  very  nearly  uniform  with 
the  phase  of  the  geometrical  optics  field  over  the  entire  aperture. 

Hence,  even  though  essentially  no  correction  of  the  amplitude  of  this 
lobe  was  obtainable,  it  is  possible  to  substantially  increase  the  gain 
and  efficiency  of  the  antenna  by  this  method.  As  mentioned  in  Chapter  3 
this  correction  could  potentially  increase  the  effective  aperture  area 
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by  the  order  of  30%  for  a  2  A  diameter  subreflector  and  on  the 
order  of  10%  for  a  6A  diameter  subreflector. 

The  method  used  for  correcting  the  shadow  boundary  lobe  is  by  no 
means  limited  to  the  degree  of  success  so  far  obtained  and  discussed 
above.  While  a  100%  correction  of  this  effect  does  not  seem  possible, 
substantial  correction  of  both  the  amplitude  and  phase  does  appear  pos¬ 
sible  by  continued  judicious  application  of  the  technique. 


The  method  to  be  used  for  correction  is  essentially  an  iterative 

approach  to  the  solution  with  only  one  iteration  taken.  It  will  be  clear 

from  the  computed  results  presented  later  that  repeated  application 

or  iteration  of  this  method  shows  no  promise  of  uniform  convergence, 

or  at  least  rapid  convergence.  This  is  to  be  expected,  since  it 

appears  very  likely  that  an  exact  solution  does  not  exist  in  any  case. 

The  phasor  diagram  (a)  illustrates  how  the  subreflector  edge  field 

adds  to  the  geometrical  optics  field  to  give  the  total  field  amplitude 

and  phase  C^.  (C^  will  be  denoted  as  the  total  phase  of  the 

resultant  field  and  normalized  to  a  length  approximately  equal  to  that 

of  C  .  We  designate  C„  as  the  unnormalized  phase.  Hence  we 
p  °  R 

get  CR/k  <  2tt  and  tankCj  =  tankCR).  This  phasor  diagram 
applies  to  a  specific  value  of  x!,  in  the  aperture.  If  we  presume 
that  and  C  will  remain  essentially  unchanged  upon  a  small 

deformation  of  the  reflector  shapes,  then  we  may  apply  the  synthesis 
method  of  Chapter  1  to  obtain  an  aperture  distribution  given  by  J)  and 
C  as  indicated  in  the  phasor  diagram  (b).  If  4L  and  Ce  actually 
did  remain  unchanged  for  this  aperture  point,  x'^,  then  we  would 
expect  to  obtain  the  original  geometric  optics  field  X)3p  and  C 


^  W  or  h  Cft 
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at  this  point.  We  will  define  the  original  geometric  optics  field  at 
the  "desired"  distribution  we  would  like  to  obtain  finally. 


Upon  synthesizing  for  A)( x ' )  and  C  ( x )  we  change  the  reflector 
shapes  from  y^(x^)  and  y^x?)  to  the  new  shapes  y  ]n(xpn)  and 
y^ n^X2n^  anc*-  aS  a  cons  ecIuence  >  the  new  edge  diffracted  field  in  the 
aperture,  and  C  ,  is  not  equal  to  the  former  edge  field.  As 

might  be  expected,  we  find  in  general  that  the  correction  method  is 
most  sensitive  to  the  changes  in  C  to  C  .  This  will  be  evident  in 
the  computed  results  presented  later.  However,  it  is  precisely  for  this 
reason  that  the  correction  obtained  with  the  total  resultant  phase  distri¬ 
bution,  is  very  satisfactory  even  when  the  resultant  amplitude  dis¬ 

tribution,  lDC  ,  is  not  substantially  changed.  This  may  be  seen  from 

Kb_L  A  fl 

a  phasor  diagram  wherein  is  approximately  equal  to  i  but 


wherein  C  and  C  may  be  greatly  different.  F or  a  very  large 
e  en  b 

range  of  ^en  we  will  find  a  good  phase  correction  even  when  a  poor 
amplitude  correction  results.  This  is  illustrated  in  the  phasor  dia¬ 


gram  (c)  wherein  the  phasor 


connected  to 


<Q  is  al 


lowed  to 


swing  over  a  wide  arc  of  a  circle.  The  resultant  phasor,  regardless 
of  the  amplitude,  ends  up  with  a  substantially  corrected  phase.  The 
computed  results  presented  later  will  illustrate  this  effect., 


There  is  no  reason  why  we  must  start  the  synthesis  for 
and  y  (x  )  with  the  same  boundary  conditions  from  which  y^(x^) 
was  computed.  It  is,  in  fact,  desirable  to  adjust  the  boundary  condi¬ 
tions  so  that  C  ^  C  near  the  shadow  boundary  lobe,  namely  at 
e  en 

x1  s  .  An  adjustment  of  x.  bv  a  fraction  of  a  wavelength  makes 

2  2m  J  An) 

it  possible  to  obtain  C  and  C  equal  to  the  extent  that 


V/' 


y  \kt 

(«> 


COMPLEX  PLAME 
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tan  kC  «.  tan  k  C  (i.  e , ,  k  C  +  mi  «kC  ),  Since  x% 

e  p  '  e  p'  Zm 

is  always  very  large  in  wavelengths,  this  alteration  of  x.,  can  have 

Ct  XTa 

no  significant  effect  on  the  general  properties  of  the  antenna  aside 

from  the  fact  that  since  k  C  +  mr  ^k  C  we  will  find  that 

e  p 

tan  k  C  tan  k  C  c^rtan  k  C  .  Hence,  in  the  vicinity  of  the 
en  e  p  3 

shadow  boundary  lobe,  we  can  insure  that  a  point,  specifically  the 
boundary  point,  will  have  C  «  C  (or  just  as  well,  k  C  C  -+•  tt 

if  desired).  This  particular  adjustment  does  not,  unfor tunately,.  lead 
to  any  substantial  correction  of  the  shadow  boundary  lobe  in  amplitude. 
It  merely  causes  the  shadow  boundary  lobe  to  remain  in  essentially 
the  same  position  in  the  aperture  before  and  after  correction.  In  other 
words,  this  change  in  boundary  condition  causes  the  new  main  reflector 
shape  to  be  less  altered  in  shape  in  the  vicinity  of  the  shadow  boundary 
lobe.  This  adjustment  of  the  boundary  does  lead,,  on  the  other  hand, 
to  a  substantially  improved  phase  correction,  as  will  be  seen  later. 

For  the  purposes  of  correction  in  this  chapter,  we  will  choose  the 
parabola-hyperbola  system  used  in  Chapter  3  for  analysis  purposes. 
Hence  y^x^)  and  y^(x^)  will  be  respectively  a  parabola  and  hyper¬ 
bola.  We  will,  however,  alter  somewhat  what  we  consider  to  be  stand¬ 
ard  input  parameters  or  boundary  conditions.  The  change  will  be  in 
the  value  of  x£m  chosen.  As  Standard  Set  No.  2  input  parameters, 
we  will  choose  the  following: 


is  any  integer,  usually  positive,  since  usually  C 


C 


e 


n 


188 


STANDARD  SET  NO,  2  PARAMETERS 

Same  as  Standard  Set  No,  1  except  that: 

=  Maximum  value  indicated  on  Graphs. 

2  m  L 

OPTION  =  -1  is  the  standard  option  indicating  that 

for  the  new  reflectors  is  computed  so  that 

k  C  ^  k  C  at  xi  3  x  *  If  x-  is  not 
e  en  2  2  m  2m 

chosen  this  way,  then  an  OPTION  =  +l  will  be 
indicated  where  appropriate,, 

As  indicated  in  Chapter  3  also,  the  above  Standard  Set  No.  2 
parameters  may  be  assumed  in  any  discussion  or  graphical  illus¬ 
tration  unless  some  modification  is  specifically  mentioned. 


As  was  mentioned  earlier,  an  improved  correction  of  the  shadow 
boundary  lobe  may  be  possible  if  the  amplitude  and  phase,  Js)  and  C  , 
which  are  synthesized  by  the  method  c*7  Chapter  1,  are  altered  somewhat 
in  anticipation  of  the  fact  that  the  altered  reflector  shapes  v/ill  make 


and  C  different  than  Q0,  and  C  .  Such  a  change  in 
C7  3en  en  3e  e 

and  C  should  be  made  judiciously  after  a  careful  investigation.  This 
change  is  not  attempted  here.  Inadvertantly,  however,  in  the  process 
of  making  machine  computations,  the  value  of  C  was  changed,  in 
other  words  this  phasor  was  rotated  in  a  manner  described  later  in 
this  chapter..  The  results  of  these  computations  are  presented,  since 
it  is  interesting  that  although  the  specific  values  of  the  computed  varia¬ 
bles  is,  of  course,  different  than  the  case  with  ^  JkC)1  not  rotated, 
the  results  are  in  all  qualitative  aspects  the  same  as  before  rotation. 
Computed  results  that  are  illustrated  in  the  figures  of  this  chapter  will 
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be  denoted  by  ROT  on  graphs  if  such  a  rotation  was  made. 


Before  presenting  the  computed  results,  we  will  develop  the  mathe¬ 
matics  of  the  correction  procedure.  Several  approximations  which 
afford  a  simplification  in  numerical  computations  are  indicated  in 
this  development. 

II.  Application  of  the  Synthesis  Method 

Before  we  actually  apply  the  synthesis  method  to  correct  the  edge 
field,  we  will  compute  the  total  field  amplitude,  phase,  and  phase 
derivative  in  the  aperture. 


The  relative  phase,  ,  of  the  total  uncorrected  aperture  distribution 
is  easily  found  as 


'R 


—  arctan 
k 


sin  k  C 


/  3e  \ 

+  <9  3p~  sin  k  c< 


d 


cos  k  C_  +  ( )  cos  k  C( 


(1) 


It  is  convenient  for  purposes  of  comparison  and  computation  to  add  an 
integral  number  of  half  wavelengths  to  so  that  is  approxi¬ 

mately  equal  to 


C 

pm 


C  (x'  ) 

p  2m' 


C 

P 


(2) 


We  define  this  "normalized"  phase  as  Gj  and  find  its  value  by 


=  C 


R 


pm 


-  1/k  arctan  ^:an  (k  C„m)  ) 


(3) 


The  last  two  terms  on  the  right  hand  side  of  (3)  add  up  to  an  integral 
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number  of  half  wavelengths  and  thus  provide  the  "normalization" 
that  is  sought. 


Now  the  rotation  of  ,  which  was  discussed  in  the  introduction, 
and  is  indicated  in  the  figures  by  the  designation  ROT,  was  made 
by  choosing,  inadvertently, 


C,  =  CD  +  C  1/k  arctan  (k  C  )  .  (3a) 

1  R  pm  '  pm  '  1 

The  amplitude  of  the  total  field  is  found  as 


:eR  *  bp  +  2  J  bp  '  be’  cos  (k  tcp-Ce])  +  IJe  .  (4) 

These  formulas  are  easily  derived  from  the  phasor  diagrams  pre¬ 
sented  earlier.  Note  that  we  continue  the  convention  that 


1eR  =  4r  ’  XRSL  =  ^RSL*  *3p  =  ^3p  *  ^e  “  ^3e  ’  etc>  1 


In  order  to  obtain  C  and  <0-  >  the  phase  and  amplitude  for  which  we 
will  synthesize  new  reflector  shapes,  we  either  let  C  — +  Wk 
-  4.  -  -  4.  ,  with  the  result  that  this  negative  edge  field  added 
to  the  geometrical  optics  field  gives  for  the  phase 


kc„-4e> 


C  =  1/k  arctan  |+  C 


pm 


Lcos  k  Cp'(^)  cos  k  CeJ 

-  1/k  arctan  (tan  k  Cpm) , 


(5) 


and  for  the  amplitude 


I  =  I3p  -  2  |  I3p  *  l3e'  cos  (k  [cp  -  ceJ)  +  l3e  • 


l3p 


(6) 


Note  that  the  value  of  C  must  be  normalized  by  the  addition  of 


pm 


1/k  arctan  (tan  k  Cpm) 


in  order  that  the  synthesis  yield  a  usable  result.  Note  also  that 


when  tan  k  C  =  tan  k  C'  then  sin  k  C  =  ±  sin  k 


and 


coskCp  =  ±  cos  k  Ce  and  equation  (5)  reduces  to 


C  =  C 


pm  * 


(7) 


This  holds  true  for  the  standard  OPTION  =  -1  when  tan  k  C  = 

P 

tan  k  C  at  x->  = 

e  2  2m 


In  order  to  actually  carry  out  the  synthesis  we  still  need  the 
derivative  function 

dC  _  _dC_ 

’  dx2 

where  the  value  in  the  above  derivative  is  that  point  in  the 
aperture  which  is  designated  in  the  uniform  phase  case  as  =  x^. 

A  different  value  of  x^  exists  for  the  generally  non-uniform  phase 
distribution  C  but  is  not  needed  in  the  derviations  we  shall  consider. 
This  may  be  made  somewhat  clearer  by  the  diagram  below  (d) 
wherein  the  ray  passing  normally  through  the  aperture  belongs  to  C 
for  which  x  =  x'  and  the  ray  passing  obliquely  through  the  apertuie 
belongs  to  C  and  has  the  same  xi>  value  but  is  incident  upon  the  main 
reflector  at  a  point  different  in  value  than  xi,. 

d  G3 

As  an  approximation  in  obtaining  we  will  assume  that 

<&3p  ~  d|3e 

d  x^  dx^  * 


C-ray 

'  Cp-  ray 
— - ►*- 

X £  for  botr?  rays 

for  Cp  -  raj 

for  C~ra y  onftj 
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From  the  results  illustrated  in  chapter  three  it  is  evident  that 
(or  ^p)  is  a  very  smooth  function  in  comparison  with-S3g 

<or  i3e). 


This  is  especially  true  when^^e  is  large  or  on  the  order  of  ^p* 

When  A 

3e  ^  s3p  the  approximation  is  not  as  good  but  it  also  is 
not  necessary  here  since  no  correction  is  needed,  Hence  we  may 
state  that 

d  T 


d  /^3e_; 
P 


dxzV  $3 


3 


1_  dfi3e 
3p  ~3TT 


1 


3p7 


“3e 


2/^3e_N' 


(8) 


In  practice,,  one  of  the  most  important  cases  to  which  this  correction 
may  be  applied  is,  incidently,  the  case  wherein  =  constant  and 
/  d  =  0.  In  the  present  case  we  have,  identically)  that 


dC 


-P  =  0 


(9) 


We  now  must  evaluate 


d  ^"3  c  d  C 

^ —  explicitly  in  order  to  obtain  g-— 


Referring  to  equation  (90)  of  chapter  three,  we  obtain 


dl- 

3e 

cTxT 


=  I 


where 


"d  x. 

^r. 


ctnP. 


te) 


(10) 


is  defined  by  equation  (87)  in  chapter  three.  From 


equation  (39)  of  chapter  three  we  obtain  by  differentiation 
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dJ ' 

t 

tt: 


±  [■ 


G2  I.  G 
e  lm 


2 1 

lm  J 


E 


Y !  *  Y j  c o s 2^  +  Y2  Y 'z  sin2£  +Y2  Y2  sin$  cosf-f 


(11) 


The  remaining  terms  in  (il)  are  defined  in  chapter  three., 


d  x 

In  order  to  obtain  2  2 

drz  d  *2 


we  differentiate 


Q) 


and  get 


d^  x2 

^2  d  X2 


'  .2 

d  X3 

^2  dX2 


cf  x  d 


d  x. 


^7 


fa  \/d‘  x. 


tan 

dr  ^ 


r3  w-  a3  \  /dy_\  /dx  \  ,  .  . 

V^TvUr,  a?-,)  tan2r3  +(dij)(dr“)sec  ^3  (a^“ 


/dV 

\^r. 


"djr3\  2  d2r3 

yW7J  sec  * 3  +  y3 


sec  r 


'd/ 

•3  v^r 


3  V  2 


+  ^y,  l  TTpr—  )  sec "jr\  tan  df 


2  / 


d^3 
3  V  dx. 


(12) 


Most  of  the  quantities  in  equation  (12)  are  defined  already  in  chapter 
three.  That  is,  in  chapter  three  we  have  y  2  defined  by  (85), 
by  (80),  %  by  (81),  d*3  by  (88),  and  dx3  by  (89).  The 

d*3  d_7T  dr2 

remaining  new  differentials  are  readily  found  in  terms  of  those  already 
defined.  We  have 


d r. 


dir. 


dx. 


<1x1 


w. 


dr. 


(13) 


'dy.t 


directly.  By  differentiating  (  g— — \  we  get 
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where  f^  is  the  main  reflector  focal  length  derived  in  Chapter  3. 
Differentiating  (dy^/dx.^  we  get 


dZ  i  ,/dx3^ 

dx3 d p2  '  ’  TTr  v3^ ) 


where 


F  rom  (16)  we  get 


dy3 


(15) 


(16) 


(17) 


The  final  remaining  derivative  is  obtained  by  differentiating 
dx3  /  ar2)  with  the  result  that 


dZ  x3 
dr2  dx2 


:tn  r  . 


•  2n 
sm 


/dL 


dx,, 


-  1 


(18) 


With  the  above  derivatives  defining  (  "3e  )  completely  and 

\  dx”  / 

(dCe  /dx.,)  already  defined  in  Chapter  3  by  equation  (94)  we  can 
obtain  ^dC/dx^  ,  For  notational  convenience  we  define  first 


(19) 


kQ  =  kC  -  kC 

pm 


+  arctan  (tan  (kCpm)  ), 


D  =  cos  k  C 
nc  p 


cos  k  C  , 
e 


and 


from  equation  (8)  above.  We  find  then  for  (dC/dx^): 


(20) 


(21) 


We  are  now  in  a  position  to  use  the  results  of  Chapter  1  and  synthesize 
new  reflector  cross-sections  y.,  (x0  )  and  y,  (x,  ).  The  variables 

in  the  new  dual  reflector  system  will  be  subscripted  by  the  letter  n 
with  the  exception  of  xi,  .  This  variable  remains  the  same,  by 
definition,  in  both  systems.  Hence  a  reference  to  a  point  xJ,  in  the 
aperture  applies  to  both  reflector  systems.  The  diagram  below  (e) 
illustrates  that  the  primary  source  rays  for  each  reflector  system 
have  been  chosen  to  cross  in  the  aperture.  This  is  convenient  for 
comparison  of  the  results  when  they  are  tabulated  and  in  programming 
for  machine  computations. 
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The  boundary  conditions  for  the  new  two  first  order  differential 
equations,  0jn  and  y1^  ,  will  remain  unchanged,  that  is: 

0 ,  _  =  (initial  0  ,)  =  arctan  (x.  /  R  )  (23) 

1  max  '  1'  '  1  max  J~J '  '  J> 

and 

^  1  max  =  (initial  yj  =  °-  <24) 

Tiiese  are  the  conditions  which  would  have  been  used  if  the  parabola- 
hyperbola  system  (y^  (x^)  and  y^  (x^)  )  had  been  derived  by  means 
of  the  synthesis  method  of  Chapter  1.  The  boundary  conditions  are 
consistent  with  Standard  Set  2  parameters. 

The  two  simultaneous  first  order  differential  equations  for  the  new 
reflector 


and 

yin 


where  the  pertinent  variables  for  the  system  are  defined  as  set  forth 
in  Chapter  1: 

Xln=  </S-  Vln>  l“  ,  •  <27> 


system  are 


’  2 
sin  0 


ln> 


(25) 


2  .  .  (*ln  - 


ij!$m  Yln)  sec  Qln  •  tan  \ 


In 


1  +  tan  6  j  .  tan(— ^2 - — 


) 


(26) 


B  =  C  -  (/  -  yln)  sec  0^ 


(28) 
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x 


2n  ” 


(29) 


*2* 


■  -CO  ( 

(xtn  -*i)z  -n„>2  ~b2J 

( xt  -  xl)  +  B 
(V- v.  p  ^dx2'  ln  2 

'  lli^  - - — - - 

/  &  J 

0 


2n  =  arctan 


x 


2n 


In 

-  y 


(30) 


(31) 


0O  =  arc sin 
3n 


(*C  \ 

\  x2  J 


(32) 


and 


^2n 


d  x 


2n 


The  value  of 


is  not  necessary  for  integrating  (25)  and 


(26)  but  will  be  used  to  find  and  subsequently 


(33) 


The  values  of  C  ,  I0  ,  and  I„  OT  will  not  be  computed  at  x' 
en  3en  RSL  ^ 

but  at  a  different  point  in  the  aperture  defined  as  x.^.  •  This  is  done 

for  computational  convenience.  Although  in  principle  we  could  compute 

these  variables  at  xJL  and  not  define  x'  ,  this  would  involve  a 

2  2n 

considerably  more  complicated  computer  program.  Furthermore 
we  find,  upon  making  the  computations,  that  x^  ^  x^n  and  in  the 
region  of  the  shadow  boundary  lobe  the  difference  between  xi,  and 
is  negligible.  The  value  of  x£  is  defined  by 


20( 


x" 


2n 


•  x2n+  y2n  tan  (  ^3n' 


(34) 


where 


n  /dy 

^3n  *  "zn  +  2  arctan  l  35 


3n 


3n 


(35) 


P2n  =  arctan  (x^  /  (y^  +  V  )  ), 


3n 


3n 


(36) 


and 


X3n  ~  X2n 


y3n  y2n 


dy 


3  n 


dy 


Zn 


dx 


3n 


dx 


Zn 


(37) 


The  edge  source  ray  making  an  angle  P^n  with  the  and  y^  axes 
is  thus  seen  to  intersect  y^n  (x^  )  at  P°*nt  (x2n>  ^2n^‘  ^*nce  ^ 
©2  >  then  the  reflected  ray  intersects  the  aperture  not  xi,  but 
at  the  newly  defined  point  x^  .  These  rays  are  illustrated  in  the 
diagram  below  (f) . 

In  order  to  compute  Cen  at  x^n  we  first  compute  y^4n  andy^n  to 
that  point  and  find 


2n 


/^4n  =  /*3^  +  (v+  y^)z 


r  3n; 


(38) 


and 


ft*  =  y3n  /  cos  . 


(39) 


We  find  R^,  the  reference  phase  (or  CeQn)  by  using  P  from 
(36)  and  the  results  of  Chapter  3  (equation  (41)  in  that  chapter). 
With  these  values  we  obtain 

C  *  /^n  +^5n  -  R 
en  •  I  n 

or  alternatively  (since  a8  exP*ained  Chapter  3) 

C  *  n  f  ^5n  +  C 

cn  '  '  eon. 


From  equation  (39)  of  Chapter  3  we  obtain  the  value  of  J^n  (P ^  ). 
In  order  to  obtain  1^  from  the  energy  conservation  equation 


i.  « r 

3en  .en 


we  must  first  find 
numerically  from 


The  value  of  this  derivative  is  found 


dx2n 


2nLS 


2n 


*T2 


2n 


2nLS 


2n 


■) 


where  x'2nLS  311,3  ^2nLS  are  points  adjacent  to  x'^n  and  ^n'  FoT 
the  purpqses  of  the  computations  made  in  this  report,  the  interval 

xi  -  xl  .  o  was  almost  a  constant  0.  05  (note  that  xi  and  not  x' 

2n  2nLo  ^  4,11 

is  the  independent  vafiable;  the  equivalent  interval  was  precisely 
0.  05  for  x').  Since  P~.  was  an  almost  linearly  varying  function  of 
xi  also,  the  app«roximatioii  for  (d  xl  / d  P  )  and  consequently  I3en 
was  very  g«od. 
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If,  as  explained  earlier,  we  assume  that  x^  -^x^  >  we  can  calculate 

^RSL  directly-  Actually  this  assumption  is  not  necessary  for  finding 

*RSL  a*  X2  ^  v-eiW  convenient  for  machine  computations.  Since 

this  approximation  was  made  in  the  computational  results  that  will 

be  presented  shortly,  we  make  the  approximation  in  these  derivations 

also.  The  approximation 'is  excellent  in  the  vicinity  of  the  shadow 

boundary  lobe  which  is  most  important.  With  this  approximation  we 

assume  that  I,  I,  ,  C.  and  C  are  all  computed  at  the  same  point 
3en  en 

instead  of  at  the  different  points  x^  and  x^.  We  find  therefore,  for  I^g^ 
that 


W  ?  1  +  2  /*  •  '3«„  c°»  k  <C  -  Ce„>  +  hen. 


(43) 


With  the  same  assumptions  made  in  the  last  paragraph  we  can 
compute  the  normalized  phase  of  I^g^,  ,  which  is  found  as 

§3—  -  N 

(0 


C2  * 


arctan 


sin  k  C  + 


>en 


sin  k  C 


en 


cos  k  C 


"X  COS  k  Cen 


(44) 


+  C  -  arctan  (tan  (k  C  )  ). 

pm  v  k  /  pm' 


Equation  (44)  completes  the  derivations  necessary  for  an  evaluation 
of  this  method  for  correcting  the  aperture  edge  field  effect.  It  might 
be  pointed  out,  in  addition,  that  the  syntheses  indicated  above  and  presented 
in  the  next  section  are  more  general  illustrations  of  the  synthesis  method 
than  was  presented  in  Chapter  1.  Neither  the  phase,  C,  nor  the  amplitude 
,  are  linear  for  these  cases. 
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IIL  Computer  Edge  Field  Correction  Results 

Extensive  computations  were  made  to  test  the  aperture  edge  field 
correction  technique  discussed  and  derived  in  Sections  I  and  II.  The 
results  of  these  computations  are  illustrated  in  Figures  1  through  47 
of  this  chapter .  For  the  most  part  the  results  depicted  in  the  figures 
are  self  explanatory  when  taken  together  with  the  discussion  in 
Section  I.  Hence  this  section  will  be  devoted  to  merely  pointing  out 
some  of  the  more  significant  aspects  of  the  results  and  just  how  the 
figures  are  organized. 

Figures  1  through  23  use  the  value  for  Cj  given  by  equation  (3), 

In  other  words,  they  do  not  have  9  exp  (kC)  rotated  in  the  complex 
plane.,  Figures  24  through  47  depict  the  "rotated11  results  and  have 
the  ROT  designation  noted  in  the  figures  .  The  amount  of  rotation, 
incidentally,  will  be  different  for  each  case  of  different  k. 

Figures  1  through  23  consist  of  four  complete  correction  attempts  for  the 
three  frequencies  k  =  2  tt ,  4tt?  and  6tt.  These  lower  frequencies  are 
the  most  significant  from  the  point  of  view  that  an  edge  field  correction 
is  necessary  for  these  small  sub-reflector  diameters  if  they  are  to  be 
of  practical  value  Note  that  a  frequency  of  k  =  utt  implies  a  sub¬ 
reflector  diameter  of  H  wavelengths  since  1*  0.  For  the  frequency 

k  =  2tt?  a  case  with  the  option  value  + 1  in  addition  to  the  case  with  the 
standard  option  value  -1  (Standard  Set  2  parameters;  see  Section  I  of 
this  chapter)  is  illustrated.  The  differences  between  the  results 
obtained  with  each  option  value  is  very  much  the  same  for  all  k.  Figures 
22  arid  23  illustrate  the  shape  of  the  s ub- reflector  before  and  after 
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correction,  (x^)  and  (x^),  f°r  rather  large  k  values, 

k  =  16tt  and  32tn  The  convergence  of  y^  (x^)  y^  (x^)  should  be 

noted. 

Figures  24  through  47  consist  of  four  complete  cases  with  k  equal  to  3Zirr 
2t r,  6tt,  and  16tTo  For  k  =  6tt  there  are  two  cases  illustrated  with 
option  =  +1  in  one  case  and  option  =  -  1  in  the  other  case.  Figure  47, 
for  frequency  k  =  32tt,  illustrates  y^  (x^)  anc^  Y  \n  (xjn)  ^or  this  high 
frequency.  The  results  for  the  rotated  cases  (ROT  =  designation  on 
Figures  24  through  47)  are  qualitatively  the  same  as  those  for  the 
non-rotated  cases  This  result  is  of  some  significance  and;  is  the  primary 
reason  for  presenting  these  figures  in  this  work. 


In  each  case  the  following  curves  are  illustrated: 


C  and  C  versus  x*  ; 
e  en  2 


S3p'  S3e'  93en  versus  x2  ; 


C',-  C'  ,  C  versus  xi  ; 

X  lL  p  Lj 

S3p'  Ser'  \sL  versM  x2  1 


and  y2  (x2),  y2n  yj  (x,),  y]n  (*ln). 


In  each  case  the  independent  variable  xi,  varies  from  x?  -  x0  to 
~  6.  This  range  covers  the  significant  portion  of  the  edge  field 
contribution  to  the  aperture.  The  value  of  is  approximately  equal 

to  1-0 y,  but  of  course  varies  with  k  for  the  standard  option  value  -1, 
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The  polarization  case  0  was  the  only  case  treated  herein.  No  change 
of  any  qualitative  significance  is  expected  for  the  ^  case. 

However  it  may  be  expected  that  the  correction  in  reflector  shapes 
will  be  different  for  this  polarization  if  the  same  correction  technique 
is  applied.  When  the  primary  pattern  has  vector  symmetry  (polarization 
angle  with  the  sub-reflector  edge  independent  of  rotation  angle)  then 
this  presents  no  problem.  The  very  significant  case  of  circular 
polarization  is  such  a  case  with  vector  symmetry.  Otherwise  some 
average  correction  or  possibly  a  simultaneous  correction  must  be 
searched  for. 

In  all  cases  the  degree  of  equality  occurring  simultaneously  between 

C  and  C  and  also  between  <90  and  is  a  measure  of  the 

e  en  3e  3en 

correction  obtainable.  Although  general  agreement  between  anc^ 

£  ^  was  obtained  the  correlation  between  C  and  C  was  poor.  This 
3en  e  en 

was  expected  since  the  amplitude  is  more  sensitive  to  changes  in  the 

second  derivative  of  the  reflector  contour  and  the  phases,  G  and  C  , 

£  e  en 

are  sensitive  to  both  position  and  slope  of  the  reflectors. 

The  correlation  between  and  was  actually  very  poor  for  the 

3e  3en 

option  =  +1  as  is  clearly  evident  from  Figure  2,  When  the  value  of 

x.  is  adjusted  so  that  tan  k  C  ^tan  k  C  at  x'  =  x  ,  the  condition 
2m  J  p  e  Z  Zm 

when  option  =  -1,  then  very  good  correlation  between  <9^  an<^  (9^en 
is  obtained.  The  improvement  is  evident  from  a  comparison  of  Figure  7 
results  with  those  of  Figure  Z, 
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The  change  in  option  choice  does  not  affect  the  agreement  between  C 

e 

and  Cen  except  near  w^ere»  as  discussed  in  Section  I,  a 

close  agreement  is  expected.  This  change  does  however  result  in  a 
considerable  improvement  of  the  total  phase  distribution  C^,  as  shown 
in  Figure  8,  The  value  of  obtained  when  option  =  +1,  Figure  3, 
is  not  nearly  as  uniform  as  that  shown  in  Figure  8  when  option  =  -1. 


The  improvement  in  total  phase  for  all  value  of  k  as  illustrated  in 
Figures  8,  13,  and  19  is  actually  excellent.  The  resulting  phase,  C^* 
for  the  new  dual  reflector  system  is  essentially  urjiforrn  for  most 
practical  considerations.  The  phase,  Cj,  before  correction  is  generally 
of  the  order  of  180°  out  of  phase  with  the  remainder  of  the  aperture  in 


the  vicinity  of  the  shadow  boundary  lobe.  This  is  true  for  all  value  of 

k  considered.  As  evidenced  by  the  amplitudes  and  shown  in 

er  3p 

Figures  4,  9,  14,  and  20,  a  substantial  amount  of  energy  is  radiated 


out  of  phase  with  the  rest  of  the  aperture  from  the  shadow  boundary 
lobe  region.  Although  the  final  amplitude,  *s  not  substantially 

different  than  the  uncorrected  amplitude,  c9^,  the  phase  correction  is 
sufficient  to  alter  the  efficiency  of  the  antenna  substantially  as  was 


discussed  earlier.  Actually,  in  the  particular  cases  stub  ted  here,  the 

distribution  ^  ^  is  considerably  more  uniform  than  either  9  or 
RSL  er 

and  would  result  in  improved  gain  for  this  antenna.  The  reason 


that  an  improvement  in  phase,  to  C^,  was  obtained  while  no  improvement 
in  amplitude,  to  was  found  was  explained  qualitatively  in 


Section  I. 
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The  reflector  shapes  before  and  after  correction  are  illustrated  in 
Figures  5,  10,  15,  16,  21,  and  22.  Of  particular  interest  in  the 
cases  of  the  main  reflector  correction  is  the  small  change  near  xl,  ^ 
x?  for  the  option  =  -1  and  the  large  change  when  option  =  +1.  Compare 
Figures  5  and  10  for  this  effect. 

Two  factors  are  worthy  of  notice  in  the  sub-reflector  correction. 

For  smaller  k  it  will  be  noticed  that  the  corrected  sub- reflector, 

Yi  (x ,  V,  is  contoured  so  that  a  great  deal  of  reflected  energy  is 
J  In  in 

directed  towards  the  upper  portion  of  the  main  reflector.  This  is  to 
help  correct  the  large  shadow  boundary  lobe  in  this  region  of  the  main 
reflector.  Notice  that  a  much  larger  segment  of  the  sub-reflector 
ilium: nates  the  same  portion  of  the  main  reflector  after  correction. 

This  effect  decreases  uniformly  as  k  increases,  as  is  evident  from  the 
k  ~  16tt  and  32tt  cases.  Figures  21  and  22, 

In  Figure  5  a  primary  ray  trajectory  was  traced  as  it  was  reflected 
from  the  sub-reflector  and  main  reflector.  The  ray  trajectory  was 
constructed  from  numerical  data  obtained  with  the  ma.chine  computed 
results  and  it  is  seen  that  Snell's  law  is  satisfied  at  each  new  reflector. 

All  the  cases  discussed  above  were  also  analyzed  for  the  case  wherein 
Q  is  rotated  (designated  by  ROT  on  the  graphs)  and  these  results  are 
similarly  illustrated  in  Figures  24  through  47.  The  above  discussion 
applies  equally  to  these  results  in  every  qualitative  respect.  While 
these  results  may  discourage  somewhat  the  hope  that  a  judiciou.: 
variation  of  the  parameters  will  lead  to  both  an  amplitude  and  phase 


correction,  such  a  parameter  study  does  seem  worthwhile  since  the 
improvement  in  antenna  performance  could  be  substantial  if  the  study 
is  successful. 
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Various  discontinuities  and  other  minor  irregularities  on  the  graphs, 
are  explained  by  accompanying  notes  for  the.  most  part.  It  is  worth 
noting  that  these  computations  taxed  almost  the  complete  storage 
capacity  of  the  IBM  7090  at  Berkeley  and  were  very  time  consuming 
and  hence  costly  to  obtain  Computation  of  all  the  variables  in 
Section  II  and  those  of  Chapter  3  where  necessary  for  one  value  of 
consumed  approximately  fifteen  seconds  of  computer  time.  Hence 
relatively  minor  inaccuracies  may  be  expected  due  to  the  .  omputer  in 
addition  to  the  approximations  made  in  the  course  of  the  derivations. 


The  above  inaccuracies  may  give  rise  to  the  result  that  <POCT  does 

R.O.L 

not  appear  to  converge  to  for  very  small^  ,  whereas  X)  does 
1  3p  7  3en  ^er 

converge  (comparatively)  to  f or  small  t  This  may  be  due  to 

the  accuracy  obtained  or  to  the  possibility  that  no  exact  solution  exists. 
If  an  exact  solution  does  not  exist,  then  an  iterative  solution  of  this  type 
will  never  converge  precisely. 


2  320-911 


FIGURE  24 
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APPENDIX  I 

EXISTENCE,  UNIQUENESS  AND  BOUNDS 


With  the  functions 


r 


dy  j 
dx. 


Yj) 


\ 


J  \ 

and  '  / 

dy2  ,  ,  ,  (  02  -  03  \  ( 

SZl  '  f2(Vy2>  =-tanV - T—j 

V  > 

determined  for  each  reflector,  it  is  important  to  know  whether  a 
solution  to  these  equations  subject  to  the  appropriate  boundary  con¬ 
ditions  exists,  and  if  so,  what  the  maximum  possible  value  of  J  y| 
is  for  each  solution,  parti  cularly  the  solution  for  y^x^). 


(I) 


The  Cauchy  Lipschitz  method  for  integrating  (1)  provides  a  useful 
way  to  determine  bounds  to  the  solution.  Briefly,,  the  method  approxi¬ 
mately  integrates  (1)  by  the  Riemann  sum 


n-1 

1 

y  =  >  f ( x . ,  y.)  Ax.  , 

7n  /  v  i  7  r  i 

i=0 

where  nAx.  =  x.  and  Ax^  is  <  0  when  the  solution  is  developed  from 
x  ,  and  also  provides  the  necessary  and  sufficient  conditions  for 
the  limit  (n-^oo)  to  exist.  If  the  limit  exists,  then  the  solution  exists. 

The  conditions  for  existence  also  provide  a  means  for  determining 
an  upper  bound  for  |  y|  . 

1.  H.  T.  Davis,  "Introduction  to  Non-Linear  Differential  and  Integral 

Equations",  pp.  88-93,  U.S.  Atomic  Energy  Commission;  September  I960 
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The  Cauchy -Lipschitz  conditions  state  that: 

I f(x,y)  must  be  uniformly  continuous.  That  is 

|f(x,  y)  -  f(x  +  h,  y)|  <  £  when  h  <  8  for  all  S>0. 

2.  f(x,  y)  must  satisfy  the  Lipschitz  condition.  That  is, 
there  exists  a  number  K  such  that 

|f(x,y)  -  f(x,y  +  h)|  <  K  h  for  all  h  -=»  0. 

3.  |  f(x,  y)  ]  must  have  a  maximum  M,  in  a  rectangular 
region  in  the  (x,  y)  plane  defined  by  x2max  "  x2min; 
+b;  -b;  such  that 


M  x 


max 


^  b 


(2) 


In  (2)  and  what  follows,  the  value  of  x  will  be  taken  as  zero  for 

'  '  mm 

convenience. 


The  first  two  conditions  will  usually  be  satisfied  since  f(x,y)  will, 
in  general,  be  single -valued,  continuous  and  have  continuous  deriva¬ 
tives  in  practical  problems.  This  is  assured  by  making  appropriate 
choices  of  any  branches  of  f(x,  y)  if  such  exist.  The  distributions 
1^(0^)  and  ^3p(x2)  maY  be  discontinuous  even  though  equation  (1) 
remains  analytic.  In  other  words,  if  discontinuities  in  are 


reflected  in  ^3p(x2^  then  the  reflectors  may  be  continuous, 


Satisfaction  of  the  last  condition  also  provides  an  upper  bound  for 
I  y  I  which  is  equal  to  b.  In  this  case,  the  solution  is  cons  trained 
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to  lie  in  the  shaded  triangular  region  depicted  in  Figure  1. 


If  f(x,y)  can  be  found  in  explicit  form  and  can  be  readily  manipulated, 


then  the  value  of  M  and  a  bound  for 


can  be  found  directly  from 


f(x,y).  Furthermore,  a  lower  bound  for  the  radius  of  curvature  of  the 

df  61 . 


reflectors  can  be  found  directly  by  maximizing 


dx 


in 


If  f(x,y)  cannot  be  found  explicitly  or  in  sufficiently  simple  form,  then 


absolute  bounds 


1 


and 


Y: 


for 


y 


and  y^  can  be  found 


from  equation  (1)  and  the  boundary  conditions.  Although 


y  i 


and 


may  be  considerably  higher  than 


Vl 


max 


and 


V: 


max 


respectively,  their  finite  values  are  proof  of  the  existance  of  a  solution. 


Absolute  bounds  for 
boundary  conditions: 


and 


yz 


will  be  found  for  two  cases  of  the 


a.  x 


2max 
type) 


>  0  (leads  to  solutions  similar  to  the  Cas  segrainian 


and 


x7  <  0  (a  negative  maximum  is  implied  and  leads  to  solutions 

La  m  ax 


similar  to  the  Gregorian  type). 


For  convenience,  it  will  be  assumed  that  0O  .  =  0,  .  -  0 . 

3  mm  I  mm 


Since  0 


3  max 


and  0 


1  max 


are  usually  given  boundary  conditions, 


Since  0 ~  -  arc  sin  (dCp/dx*)  and  Cp(x^)  is  given,  then  OJx.1,)  is 

»J  La  Li  u 

known  and  so  are  0,  .  and  0, 

.>  xxiiii  j  max. 
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it  is  clear  from  equation  (1)  that  it  is  sufficient  to  determine  an  absolute 
upper  bound  for 

dy2 


: 

0~  ,  0.,  ,  in  order  to  determine 

Z  Zb 


dx~ 


dyi 


dx  i 


and 


It  is  also  assumed  that  qC  and  ^3  are  given  boundary  conditions 
(note  that  they  may  be  negative  or  positive). 

The  appropriate  geometry  for  case  (a)  is  depicted  in  Figure  2.  Only 
the  end  points  of  the  reflectors  are  depicted.  An  extreme  ray  trajec¬ 
tory  is  depicted  from  the  phase  center  of  1^(0  to  the  point 

(y_  =  0,  x7,  =  x7  ).  The  value  for  07K  is  found  as 

c  Ci*  m  ax  c  d 


l  x2  max  \ 

0-,,  =  arc  tan  l  - — —  I  for  x 


2b 


c/2 


2  max 


(3) 


which  leads  to  a  value  for  y  as 

lb 


r 


=  x 


X 


1  max  tan  i'arc  tan  l.  ('  so  long  as 


0-,,  =  0t  and  with  xT  >  0  * 

Zb  >  1  max  Z  max 


(4) 
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If  0-.,  is  less  than  0.  then  0,  may  be  substituted  for 

2b  1  max  1  max  ' 

02b  in  (4).  The  value  for  |  y2|b  is  found  as  either 


hi 


.  1  )  (  2  max  ,  w  -  -  _ 

,  -  x0  tan  —  <  arc  tan  l — - - /  /for  6  v  0O 

b  2  max  2  I  V  OU  /  2b  >  3  max 


with  x~  >  0 

2  max 


(5) 


or 


y->  tan  -  J  0Q  V  for  0-  =  Q 

;2  b  2 max  2  v  3  max  [  3  max  >  2b 


with  x0  >0 

2  max 


(6) 


There  are  two  possible  geometries  for  case  (b)  which  are  depicted 
in  Figures  3a  and  3b.  The  geometry  in  Figure  3a  leads  to  a  value 
for  0^  as 


0_,  =  arc  tan 

2  b 


2  max  I 

a. 


with  x^  <  0, 

2  max 


(?) 


and  the  geometry  in  Figure  3b  leads  to 


Ix  -f~  x 

2  max  1  max  .  ... 

°2b  ~  arc  tan  1  - JTTft - I  "'lth  *2max«-0  (8) 


The  appropriate  value  of  aS  S*ven  ky  (?)  or  t®)»  use(^  to 

find  |  and  |  Y 2 1 f  Since  0^  and  0^  always  have  opposite  signs 

in  this  case,  the  value  for  I  y^l  ^  is  found  as 


yj  ,  =  x,  tan 

7  II  b  1  max 


0OU  +  9 1  \ 

2b  1  max  ] 

2 


/  with  x2max*“  0  <9> 
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Since  the  signs  of  0^  an d  6^  remain  the  same  for  this  case,,  as 
for  case  (a),  the  value  of  |  y^j  ^  *s  found  as  either 


2  max 


for 


0 


2b 


> 


6^ 

3  max 


and  with  x0  /  0 
2  ma  x 


(10) 


or 


x 


2  max 


tan 


L 

e„ 

3  max 

\ 

_/ 

V  2  / 

f  o  r 


0 

o  max 


>  9 


2b 


with  x..  >  0 

2  max  * 


(ii) 


A  somewhat  improved  estimate  of  an  upper  bound  for  j  y|  can  be 

2  _ 

obtained  if  the  assumption  is  made  that  the  sign  of  —4 ^ —  does 

dx 

not  change  over  the  reflector.  This  assumption  is  usually  valid  for 
reasonably  smooth  l3p(x2)  anc^  For  example,  it  is  true  for 

all  Cassegrain  parabola-hyperbola  and  Gregorian  parabola-ellipse 
designs.  It  is  also  true  in  the  computed  designs  shown  in  Chapter  1. 

d2y 

If  the  assumption  is  not  valid,  that  is  - does  reverse  s’  u  for  a 

dx^ 

small  interval  of  x  -  x  .  .  then  the  proposed  estimate  for  an 

max  mm  L  * 


upper  bound  to  |  y  |  may  be  in  error,  but  may  nevertheless  be  a  closer 

estimate  to  the  true  value  I  y  thaxi  the  absolute  bounds  previously 

•  "  max  • 

con  side  red. 


If  the  second  derivatives 

1 


d2yi 


dx 


1 


d9 j  -  de2 


26< 


and 


do  not  change  sign  then  this  implies,  for  most  cases,  that  (d0^  -  d0 ^ ) 
and  (dG^  -  d9^)  do  not  change  sign.  In  this  circumstance,,  the  differences 
of  the  angles, 


<0,  -  ) 

1  max  2ba 

and 

*02ba  "  03max* 


will,  in  general,  be  the  maximum  differences, 
an  improved  estimate  of  an  upper  bound  for  |  0 
for  cases  (a)  and  (b). 


The  value  9^ba 
J  which  will  be  found 


For  case  (a),  x2max^  an  aPProx^ma^e  uPPer  bound  for  j  0^j  ,  Q^ba* 
can  be  obtained  by  using  the  value  depicted  in  Figure  4.  This  value  is 
given  by 


X  “X 

~  .  ,  2  max  1  rna  x  \  r  x  ^ 

9  -  arc  tan  [  - - - —  )  for  x0  >  0 

2ba  \  x  A  /  2  max 


<c  +  (6 


(12) 


The  corresponding  approximate  upper  bounds  for  j  ^lmaxl  an<^ 
I  ^2  maxi’  Inlba  and  ly2lba’  are  then  given  by 


YtI 


l’ba  1  max 


tan 


e'>u  - 

Zba  1  max 


for  x 


2  max 


>  o  (13) 
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and 


ly2| 


ba  x2  max 


tan 


0,,  -  9-1 

2ba  3  max 


■j 


for  x 


2  max 


>  0  (14) 


For  this  case,  (a),  the  assumption  has  been  made  that 


d2Yl 


dx 


-£-  <  0  for  all  x1 


For  case  (b);  <  03  the  assumption  is  that 

u  max 


,2 

d  y 


dx 


1  >  0  for 

7 


1 


all  x^.  The  approximate  upper  bound  for  |  G^j  ,  ®2ba'  '’'S  ^oun0 

from  Figure  3b,  irrespective  of  the  position  of  the  focal  point  F, 

The  value  of  0^^  *s  same  as  that  given  for  0  by  equation  (9) 

since  9^  and  0^  have  opposite  signs  in  this  case.  It  should  be  added 

that  the  geometry  of  Figure  3b  occurs  less  frequently  than  that  of 

Figure  3a  in  practice  so  that  the  improved  bounds  for  y.l  are 

1  1  l*max 

usually  different  than  the  absolute  bounds.  The  signs  of  9-  and  0, 

both  change  in  this  case,  (b),  so  that  the  improved  bound  for  I  y  J 

r  J  Z  \  max 

is  given  by 


21  ba 


x 


2  max 


tan 


2  ba 


0 


3max 


(15) 


for  x~,  <  0, 

2  max 


Some  computed  bounds  which  are  applicable  to  the  numerical  solutions 
shown  in  Chapter  1  (Figures  5  and  6)  are  presented  in  the  table  below. 
The  values  of  these  bounds  would  justify  a  computation  and  are  valuable 
for  that  purpose. 


|K^)I<  m  (R. 
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APPENDIX  II 

ON  THE  AXIAE  CAUSTIC  OF  THE  EDGE  FIELD 
FOR  SURFACES  OF  REVOLUTION 

In  this  appendix  we  shall  justify  the  use  of  the  factor  G^e  which  is  the 
Gaussian  curvature  of  the  field  wavefront  in  the  axial  case.  We  will  also 
justify  the  fact  that  a  very  accurate  approximation  of  the  field  near  the  caus¬ 
tic  is  not  necessary,  since  only  a  very  small  solid  angle  is  subtended  in 
this  region. 

1.  The  Factor  G 

In  equation  (6)  of  Chapter  3,  the  amplitude  of  the  field  Ue  varies 
as  1  f  Tier"!  This  is  typical  for  the  far  zone  of  any  cylindrical 
field,  when  r  is  the  distance  from  the  source  to  the  field  point. 
Energy  is  thus  conserved  in  any  angular  sector  per  unit  length  along 
the  cylindrical  axis.  For  the  edge  field  the  caustic  or  source  point  is 
the  edge  itself,  and  as  seen  in  the  diagram  (a),  the  normal  area  sub¬ 
tended  by  a  tube  of  rays  from  the  caustic  depends  on  the  value  of  kr. 
Alternatively,  the  ratio  of  areas  is  given  by  kr^/kr^  for  areas  No.  1 
and  No.  2  in  the  diagram.  The  curvature  of  the  wavefronts  is  propor¬ 
tional  to  kr  .  We  thus  define  G^e  for  the  cylindrical  case  as 

GAe  *  '/^  • 

When  the  edge  is  part  of  a  circular  curve  as  in  the  axially  symmetric 
case,  the  ratio  of  two  areas  subtended  by  a  tube  of  rays  through  which 
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a  constant  rate  of  energy  flow  passes  is  proportional  to  in  the 

far  zone  except  near  the  caustic  axis.  Consider  the  diagram  .(b)  for 
this  case.  The  ratio  of  areas  now  obtained  is 


(k 


r2> 


kr2  sin  £  +  kxlml/(krtHkr 


sin  L  I  kx 


lm 


In  this  case  the  curvature  of  the  wavefronts  is  defined  as  the  Gaussian 
curvature  which  is 


GAe  =  1/4kr!kr  sin  PZ  +  kxlm^ 


For  the  bottom  edge  when  is  sti.ll  in  the  same  direction,  we  have 


Ae 


yJkrlkr  sin  P2  -  k*]mr  ( 


We  assume,  therefore,  that  for  the  curved  edges  of  the  axially  sym¬ 
metric  case,  the  angular  dependence  of  the  edge  field  is  identical  to 
that  of  the  straight  edge  case  (cylindrical  case)  except  as  modified  by 
the  Gaussian  curvature  given  above.  This  assumption  is  the  same  as 
that  employed  by  Keller  for  diffraction  from  a  flat  aperture  with  excel¬ 
lent  results  obtained  in  that  case. 


At  4>  =  tt  -  &L*  we  have  defined  G*i  ( 


Ae 


p 


2m 


)  = 


k ^  sin  |  0 


p2m 


2m 


for  the  axial  case,  where 


t-,  01  .We  should  note  that 

2m  2m  2  m 

is  only  proximately  to  w*thin  the  accuracy  of  the  far  field 

approximation.  This  is  explained  as  follows.  We  define  not 

as  the  angle  when  the  edge  ray  intercepts  the  edge  of  the  large 

reflector,  but  rather  the  value  of  /?,  when  intercepts  the  edge 


Axis  oi  Sy 

Cb) 
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of  the  large  reflector.  Hence,  H,  =  when  ^  ^ 

being  the  angle  the  normal  to  the  total  edge  diffracted  field  wavefront 
makes  with  the  axis),  While  this  distinction  is  really  not  very 

important  to  the  theory  of  the  analysis  and  negligible  when  the  far 
field  approximations  are  excellent,  it  nevertheless  is  important  when 
numerical  computations  are  made,,  especially  by  machine  computors. 
The  machine  will  not  o\erlook  a  small  difference  between  the  true  far 

'-S 

field  and  approximate  far  field  and  may  not  complete  the  calculation. 

Hence,  the  difference  between  the  defined  P>  and  must  be 

Cm  cm 

c  onsidered  carefully. 

2,  The  Caustic  Axis  of  the  Axially  Symmetric  Case 

In  general,  for  the  fax'  field  we  will  have  r  xlm  ^or  ec^ge 

fields.  Hence,  GA  reduces  to  the  same  expression  for  the  top 
and  bottom  edges 

GAe  =  kr  si*1/2  |  P2  | 

providing  that  is  not  small.  This  assumption  is  valid  so  long  as 


When  is  small,  we  take  as  an  approximation  for  the  above  inequality 

the  value  for  r  given  by  rc  =  •  This  is  a  rough  approxima¬ 

tion  for  the  distance  from  the  edge  to  reflector  No.  2,  when  Q  is 
small.  A  more  accurate  approximation  will  not  be  of  great  value  us  we 

More  suitable  choices  for  rc  can  be  made  for  specific  problems. 
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shall  show  shortly.  On  this  basis  we  take  for  the  value  of  G^e: 


When  >  0  and  near  zero,  we  will  obtain  an  infinity  from  G^e  due 
to  the  bottom  edge.  To  eliminate  this,  we  make  the  further  approx¬ 
imation  that 

G  kl  * 

for  both  edges . 

At  =  0,  this  value  for  G^e  is  correct  for  both  edges.  For 


this  value  is  again  correct  for  both  edges.  Hence,  it  is  only  for  a 
small  angular  sector  in  f ^  and  a  very  small  solid  angle  that  the 
last  given  value  for  G^e  is  in  error.  In  the  case  of  vector  symmetry 
for  the  axial  system,  we  have  =  0  (equation(13)  )  on  axis 

(  =  0)  and  the  total  energy  in  this  small  solid  angle  will  be  very 

small  as  compared  with  the  surface  reflected  primary  field  in  this 
same  region.  Even  when  A/  4  0  as  when  we  "average"  an  incident 
linearly  polarized  primary  field,  the  net  energy  in  this  solid  angle  is 
still  small  and  the  value  for  GAe  at  the  center  (  Q  =  0)  is  fairly 
accurate. 

As  an  example,  let  us  define  as  that  angle  at  which 
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The  total  solid  angle  s  ubtended  by  the  range  (0  j  2.  J  rgp]  )  i1 
very  small  for  j  x lm/rc  j^l.  The  percentage  solid  angle  in  this 
range  of  the  total  4tt  ste radians  is 


(  -1--  cos  )  x  100% 

2 

Hence  if  ^2R  were  as  large  as  Z0  degrees,  the  percentage  solid 
angle  would  be  only  6%,  The  percentage  solid  angle  of  that  subtended 
by  reflector  No,  Z  is  much  larger,  however,  but  flo  =  20°  is 
very  large  also.  For  ■=»  5°,  the  subtended  solid  angle  is  only 

0.  2%. 


The  above  numbers  indicate  that  a  precise  correction  for  the  fields  on 
the  caustic  is  not  necessary  and  will  not  alter  the  final  results  in  any 
significant  way.  Hence  a  reasonable  approximation  from  a  computa¬ 
tional  point  of  view  was  chosen. 


